Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stabilization”

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

Polyethylene Glycol Crowder’s Effect on Enzyme Aggregation, Thermal Stability, and Residual Catalytic Activity

Protein stability and performance in various natural and artificial systems incorporating many other macromolecules for therapeutic, diagnostic, sensor, and biotechnological applications attract increasing interest with the expansion of these technologies. Here we address the catalytic activity of lysozyme protein (LYZ) in the presence of a polyethylene glycol (PEG) crowder in a broad range of concentrations and temperatures in aqueous solutions of two different molecular mass PEG samples (Mw = 3350 and 10000 g/mol). The phase behavior of PEG–protein solutions is examined by using dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS), while the enzyme denaturing is monitored by using an activity assay (AS) and circular dichroism (CD) spectroscopy. Molecular dynamic (MD) simulations are used to illustrate the effect of PEG concentration on protein stability at high temperatures. The results demonstrate that LYZ residual activity after 1 h incubation at 80 °C is improved from 15% up to 55% with the addition of PEG. The improvement is attributed to two underlying mechanisms. (i) Primarily, the stabilizing effect is due to the suppression of the enzyme aggregation because of the stronger PEG–protein interactions caused by the increased hydrophobicity of PEG and lysozyme at elevated temperatures. (ii) The MD simulations showed that the addition of PEG to some degree stabilizes the secondary structures of the enzyme by delaying unfolding at elevated temperatures. The more pronounced effect is observed with an increase in PEG concentration. This trend is consistent with CD and AS experimental results, where the thermal stability is strengthened with increasing of PEG concentration and molecular mass. Overall, the results show that the highest stabilizing effect is approached at the critical overlap concentration of PEG.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

A Polymer-in-Salt Electrolyte with Enhanced Oxidative Stability for Lithium Metal Polymer Batteries

Lithium (Li) metal polymer batteries (LMPBs) are a promising candidate of solid-state batteries with high safety. However, rare progress has been demonstrated so far in high voltage stability of polyethylene oxide (PEO) based polymer electrolytes for LMPBs. Herein, we revived the polymer-in-salt electrolyte (PISE) strategy based on the PEO-LiFSI system with EO/Li = 8 through a dry process to avoid the contamination of the residual solvent. The obtained PISEs exhibit quite different morphologies and coordination structures which greatly enhance the oxidative stability of polymer electrolytes. P(EO)1LiFSI has a low melting temperature, a high ionic conductivity at 60 ?C and an oxidative stability of ~4.5 V vs. Li/Li+ (overcoming the inherent low oxidative stability of PEO). With an effective interphase rich in inorganic species, in combination with good stability of the hybrid polymer electrolyte towards Li metal, the cycling stability of Li||LiNi1/3Co1/3Mn1/3O2 is greatly improved. The LMPB can retain 74.4% of capacity after 186 cycles at 60 ?C under the charge cutoff voltage of 4.3 V. The findings point out a promising strategy to develop high-voltage stable polymer electrolytes for high energy-density and safe LMPBs.

Polymer-in-salt electrolyte, lithium metal battery↗
DOE OSTI2026

Composition, Activity, and Stability of IrO x Oxygen Evolution Reaction Electrocatalysts

The oxygen evolution reaction (OER) is integral to several electrochemical energy conversion and storage technologies, including carbon dioxide reduction to value added fuels, nitrogen reduction to ammonia, reversible fuel cells, rechargeable metal−air batteries, and water electrolysis to produce hydrogen. Iridium oxide (IrO x ) is widely recognized as the benchmark OER catalyst for acidic environments. Despite widespread use of IrO x catalysts, most notably in proton-exchange membrane water electrolyzers (PEMWEs), a comprehensive understanding of the physicochemical properties of commercial catalysts and the impact of these properties on both the activity and stability of these catalysts is lacking. Here, we study commercial IrO x catalysts with different physicochemical properties, three nominally considered amorphous and three rutile, to elucidate how structural and compositional variations affect OER activity and stability. Utilizing standardized aqueous electrochemical protocols, time-resolved dissolution quantification using inductively-coupled plasma mass spectrometry, and physicochemical characterization, including multiple synchrotron X-ray techniques, we systematically correlate catalyst properties with OER performance and degradation behavior aided by principal component analysis (PCA). Our results demonstrate the general trend of amorphous IrO x having higher intrinsic activity but limited stability and crystalline rutile IrO 2 having lower activity but enhanced stability against dissolution. The trends within the amorphous and rutile catalyst groups correlate with inherent material properties, including phase composition and structure, crystallinity, particle size, surface area, and surface structure/chemistry. Notably, we identify a rutile catalyst with the largest crystallite/ domain sizes, moderate surface area, a small fraction of hydrous phase, and a favorable pore structure (trimodal distributions of pore sizes ranging from 2−5 nm) that exhibits the best balance between activity and stability among the six catalysts studied here. These findings illustrate a fundamental structure-governed trade-off between activity and stability and highlight the critical role of surface chemistry modification and structure engineering in IrO x catalyst optimization.

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

Insights Into Nature-Based Climate Solutions: Managing Forests for Climate Resilience and Carbon Stability

Successful implementation of forest management as a nature-based climate solution is dependent on the durability of management-induced changes in forest carbon storage and sequestration. As forests face unprecedented stability risks in the face of ongoing climate change, much remains unknown regarding how management will impact forest stability, or how interactions with climate might shift the response of forests to management across spatiotemporal scales. Here, we used a process-based model to simulate multidecadal projections of forest dynamics in response to changes in management and climate. Simulations were conducted across gradients in forest type, edaphic factors, and management intensity under two alternate radiative forcing scenarios (RCP4.5 and RCP8.5). This allowed for the quantification of forest stability shifts in response to climate change, and the role of management in modulating that response, where ecosystem stability is characterized as the resilience and temporal stability of net primary production, aboveground biomass, and soil carbon. Our results indicate that forest structure is primarily shaped by management, but the same management strategy often produced divergent structures over time, due to interactions with regional climate change. We found that management can be used to increase stability and minimize the release of stored carbon by reducing mortality, but also highlight the regional dependency of management-induced changes in resilience to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗
DOE OSTI2024

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗
DOE OSTI2023

A Stability Analysis Tool for Bulk Power Systems Using Black-Box Models of Inverter-Based Resources

This paper presents a small-signal stability analysis tool for large-scale power systems with high penetration of inverter-based resources (IBRs). Firstly, a network transfer function matrix (NTFM), which represents the information of the system topology, transmission lines, loads, IBRs locations, etc., is derived to model the entire power system network. Secondly, small-signal perturbation method is applied to obtain the sequence impedance/admittance models of the block-box IBRs considering the frequency cross-coupling effects. With the obtained NTFM as well as IBR models, a multi-input, multi-output (MIMO) feedback system is constructed, and the generalized Nyquist criterion (GNC)-based stability method is employed to analyze the stability of the entire power system. Furthermore, based on the developed stability analysis method, sensitivity analysis is conducted on an unstable case to identify which parameter has a high impact on the system stability. As a result, different test cases based on a modified IEEE 14- bus system as well as a reduced 240-bus WECC system are studied to verify the proposed stability analysis tool.

42 ENGINEERING↗
DOE OSTI2020

Low-Frequency Stability Analysis of Inverter-Based Islanded Multiple-Bus AC Microgrids Based on Terminal Characteristics

For system planning of three-phase inverter-based islanded ac microgrids, the low frequency instability issue caused by interactions of inverter droop controllers is a major concern. When internal control information of procured commercial inverters is unknown, impedance-based small-signal stability criteria facilitate prediction of resonances in medium and high frequency ranges, but they usually assume the grid fundamental frequency as constant and thus they are incapable of analyzing the low-frequency oscillation of the fundamental frequency in islanded microgrids. Aiming at solving this issue, this paper proposes two stability analysis methods based on terminal characteristics of inverters and passive connection network including the dynamics of the fundamental frequency for analysis of low-frequency stability in islanded multiple-bus microgrids. Based on the Component Connection Method (CCM) to systematically separate inverters from the passive connection network, a general approach is developed to model the microgrid as a multiple-input-multiple-output (MIMO) negative feedback system in the common system d-q reference frame. By applying the generalized Nyquist stability criterion (GNC) to the return-ratio and return-difference matrices of the MIMO system model, the low-frequency stability related to the fundamental frequency can be analyzed using the measured terminal characteristics of inverters. Finally, analysis and simulation of a 37-bus microgrid verify the effectiveness of the proposed stability analysis methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Moisture Stability of Sulfide Solid-State Electrolytes

In this report we detail a comprehensive study on the moisture stability of sulfide solid-state electrolytes in dry room environments. Although sulfide SSEs have many favorable attributes, this class of materials suffers from poor stability with water. Sulfide SSEs react with water to form gaseous H 2 S and a variety of solid byproducts like Li 3 PO 4 and LiOH, which go on to increase the interfacial impedance of solid-state batteries. Lab-scale research typically utilizes gloveboxes with <1 ppm water, however, the large-scale manufacturing of Li-ion batteries occurs in –40°C dewpoint dry rooms with around 126 ppm water. Consequently, the moisture stability of sulfide SSEs must be addressed if the manufacture of solid-state batteries based on sulfide SSEs is to be scaled up. Here, we are the first to characterize the moisture stability of sulfide SSEs according to both H 2 S and the degradation of ionic conductivity at different moisture setpoints ranging from –76°C to –40°C dewpoint. A variety of different SSE compositions are studied; namely, (Li 2 S) 75 (P 2 S 5 ) 25 , (Li 2 S) 70 (P 2 S 5 ) 30 , (Li 2 O) 7 (Li 2 S) 68 (P 2 S 5 ) 25 , (Li 2 O) 7 (Li 2 S) 63 (P 2 S 5 ) 30 , and (Li 2 S) 75 (P 2 S 5 ) 25 + 20 mol% LiI. We find that moisture stability improves with 75 mol% Li2S modifier content and the introduction of a Li 2 O co-modifier. After a 30 min exposure in a –40°C dewpoint dry room environment we found that (Li 2 S) 75 (P 2 S 5 ) 25 + 20 mol% LiI powder generated 0.1 cc/g H 2 S and its ionic conductivity decreased by over 50%. However, when SSE powder was exposed as a slurry in a dodecane carrier the same SSE composition generated 0 cc/g H 2 S and its ionic conductivity only dropped by 14%. Our results show that sulfide SSEs have acceptable moisture stability when appropriately processed in a dry room environment.

25 ENERGY STORAGE↗
DOE OSTI2020

Analyzing Impact of Distributed PV Generation on Integrated Transmission & Distribution System Voltage Stability—A Graph Trace Analysis Based Approach

The use of a Graph Trace Analysis (GTA)-based power flow for analyzing the voltage stability of integrated Transmission and Distribution (T&D) networks is discussed in the context of distributed Photovoltaic (PV) generation. The voltage stability of lines and the load carrying capability of buses is analyzed at various PV penetration levels. It is shown that as the PV generation levels increase, an increase in the steady state voltage stability of the system is observed. Moreover, within certain regions of stability margin changes, changes in voltage stability margins of transmission lines are shown to be linearly related to changes in the loading of the lines. Two case studies are presented, where one case study involves a model with eight voltage levels and 784,000 nodes. In one case study, a voltage-stability heat map is used to demonstrate the identification of weak lines and buses.

14 SOLAR ENERGY↗
NASA NTRS1952

Wind-tunnel Investigation of the Contribution of a Vertical Tail to the Directional Stability of a Fighter-type Airplane

Contributions of vertical tail and fuselage to directional stability are determined from pressure distributions measured over wide ranges of angle of attack and angle of yaw with stabilizer in each of three vertical positions, with stabilizer removed and with both stabilizer and vertical tail removed. Results point out inadequacies of current design methods and show that more accurate methods must treat separately contributions of vertical tail, fuselage area above stabilizer, and fuselage area below stabilizer.

LOADS, STEADY - TAIL↗
NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗
DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Interphase Stabilization of LiNi 0.5 Mn 1.5 O 4 Cathode for 5 V–Class All–Solid–State Batteries

Employing high voltage cobalt-free spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) as a cathode is promising for high energy density and cost-effectiveness, but it has challenges in all-solid-state batteries (ASSBs). Here, it is revealed that the limitation of lithium argyrodite sulfide solid electrolyte (Li 6 PS 5 Cl) with the LNMO cathode is due to the intrinsic chemical incompatibility and poor oxidative stability. Through a careful analysis of the interphase of LNMO, it is elucidated that even the halide solid electrolyte (Li 3 InCl 6 ) with high oxidative stability can be decomposed to form resistive interphase layers with LNMO in ASSBs. Interestingly, with Fe-doping and a Li 3 PO 4 protective layer coating, LNMO with Li 3 InCl 6 displays stable cycle performance with a stabilized interphase at a high voltage (≈4.7 V) in ASSBs. The enhanced interfacial stability with the extended electrochemical stability window through doping and coating enables high electrochemical stability with LNMO in ASSBs. Furthermore, this work provides guidance for employing high-voltage cathodes in ASSBs and highlights the importance of stable interphases to enable stable cycling in ASSBs.

25 ENERGY STORAGE↗
DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗
DOE OSTI2020

Efficient and robust stability analysis in the internal energy, volume, and moles (UVN) space

Hydrodynamic simulation of fluids in vessels is often based on the internal energy (U), volume (V), and the number of moles of all components ( re N = N 1 , N 2 , ..., N n ) as the working variables to determine the thermodynamic state by finding the global maximum of entropy. The procedure of finding the equilibrium state may be divided into two steps: stability analysis and phase-split computations. Stability analysis is performed first, and the phase-split computation proceeds if the system is unstable. Most of the past work on UVN space has been based on phase-split computations, while the stability analysis is barely investigated. Furthermore, previous studies for the stability analysis use only Newton's method to solve the non-linear algebraic equations. We present an efficient and robust approach for stability analysis where UVN is specified. The successive substitution iteration (SSI) is used to provide good initial guesses for Newton's method. The proposed approach results in a reduced number of unknowns and does not require a large number of iterations to achieve convergence in Newton's method. Overall, our proposed formulation is compatible with different equations of state and is applicable to both pure component and multi-components. The robustness and efficiency of the algorithm for stability analysis are demonstrated in various examples.

02 PETROLEUM↗
DOE OSTI2020

Uncertainty evaluation for twist drilling stability model

This paper describes the first uncertainty analysis for drilling stability using a frequency-domain drilling stability model. The stability model inputs include: the modal parameters for the torsional-axial vibration mode from the twist drill-holder-spindle axial frequency response function; and the mechanistic coefficients that relate the torque and thrust force to chip area for the selected drill-workpiece material combination. Furthermore, Monte Carlo simulation is applied to propagate the input uncertainties to output uncertainty in the predicted stability map, which separates stable from unstable (chatter) zones in the spindle speed-chip width parameter space. Additionally, the mean stability boundary and its 95% confidence intervals are determined for five cases: varying all four inputs simultaneously and varying them individually. This enables the individual sensitivities to be compared. Experimental results from drilling tests are included for comparison to the prediction. Additionally, Matlab code is provided to implement the stability model and Monte Carlo uncertainty analysis.

42 ENGINEERING↗