SEARCH · Engineering Papers

Results for “monotonicity”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Correlation between Mechanical Behavior and Actuator-type Performance of Ni-Ti-Pd High-temperature Shape Memory Alloys

High-temperature shape memory alloys in the NiTiPd system are being investigated as lower cost alternatives to NiTiPt alloys for use in compact solid-state actuators for the aerospace, automotive, and power generation industries. A range of ternary NiTiPd alloys containing 15 to 46 at.% Pd has been processed and actuator mimicking tests (thermal cycling under load) were used to measure transformation temperatures, work behavior, and dimensional stability. With increasing Pd content, the work output of the material decreased, while the amount of permanent strain resulting from each load-biased thermal cycle increased. Monotonic isothermal tension testing of the high-temperature austenite and low temperature martensite phases was used to partially explain these behaviors, where a mismatch in yield strength between the austenite and martensite phases was observed at high Pd levels. Moreover, to further understand the source of the permanent strain at lower Pd levels, strain recovery tests were conducted to determine the onset of plastic deformation in the martensite phase. Consequently, the work behavior and dimensional stability during thermal cycling under load of the various NiTiPd alloys is discussed in relation to the deformation behavior of the materials as revealed by the strain recovery and monotonic tension tests.

Bigelow, Glen S.↗

NASA NTRS2006

Dependence of Microelastic-plastic Nonlinearity of Martensitic Stainless Steel on Fatigue Damage Accumulation

Self-organized substructural arrangements of dislocations formed in wavy slip metals during cyclic stress-induced fatigue produce substantial changes in the material microelastic-plastic nonlinearity, a quantitative measure of which is the nonlinearity parameter Beta extracted from acoustic harmonic generation measurements. The contributions to Beta from the substructural evolution of dislocations and crack growth for fatigued martensitic 410Cb stainless steel are calculated from the Cantrell model as a function of percent full fatigue life to fracture. A wave interaction factor f(sub WI) is introduced into the model to account experimentally for the relative volume of material fatigue damage included in the volume of material swept out by an interrogating acoustic wave. For cyclic stress-controlled loading at 551 MPa and f(sub WI) = 0.013 the model predicts a monotonic increase in Beta from dislocation substructures of almost 100 percent from the virgin state to roughly 95 percent full life. Negligible contributions from cracks are predicted in this range of fatigue life. However, over the last five percent of fatigue life the model predicts a rapid monotonic increase of Beta by several thousand percent that is dominated by crack growth. The theoretical predictions are in good agreement with experimental measurements of 410Cb stainless steel samples fatigued in uniaxial, stress-controlled cyclic loading at 551 MPa from zero to full tensile load with a measured f(sub WI) of 0.013.

Cantrell, John H.↗

NASA NTRS2010

A Singular Perturbation Approach for Time-Domain Assessment of Phase Margin

This paper considers the problem of time-domain assessment of the Phase Margin (PM) of a Single Input Single Output (SISO) Linear Time-Invariant (LTI) system using a singular perturbation approach, where a SISO LTI fast loop system, whose phase lag increases monotonically with frequency, is introduced into the loop as a singular perturbation with a singular perturbation (time-scale separation) parameter Epsilon. First, a bijective relationship between the Singular Perturbation Margin (SPM) max and the PM of the nominal (slow) system is established with an approximation error on the order of Epsilon(exp 2). In proving this result, relationships between the singular perturbation parameter Epsilon, PM of the perturbed system, PM and SPM of the nominal system, and the (monotonically increasing) phase of the fast system are also revealed. These results make it possible to assess the PM of the nominal system in the time-domain for SISO LTI systems using the SPM with a standardized testing system called "PM-gauge," as demonstrated by examples. PM is a widely used stability margin for LTI control system design and certification. Unfortunately, it is not applicable to Linear Time-Varying (LTV) and Nonlinear Time-Varying (NLTV) systems. The approach developed here can be used to establish a theoretical as well as practical metric of stability margin for LTV and NLTV systems using a standardized SPM that is backward compatible with PM.

Zhu, J. Jim↗

NASA NTRS2010

Efficient Web Services Policy Combination

Large-scale Web security systems usually involve cooperation between domains with non-identical policies. The network management and Web communication software used by the different organizations presents a stumbling block. Many of the tools used by the various divisions do not have the ability to communicate network management data with each other. At best, this means that manual human intervention into the communication protocols used at various network routers and endpoints is required. Developing practical, sound, and automated ways to compose policies to bridge these differences is a long-standing problem. One of the key subtleties is the need to deal with inconsistencies and defaults where one organization proposes a rule on a particular feature, and another has a different rule or expresses no rule. A general approach is to assign priorities to rules and observe the rules with the highest priorities when there are conflicts. The present methods have inherent inefficiency, which heavily restrict their practical applications. A new, efficient algorithm combines policies utilized for Web services. The method is based on an algorithm that allows an automatic and scalable composition of security policies between multiple organizations. It is based on defeasible policy composition, a promising approach for finding conflicts and resolving priorities between rules. In the general case, policy negotiation is an intractable problem. A promising method, suggested in the literature, is when policies are represented in defeasible logic, and composition is based on rules for non-monotonic inference. In this system, policy writers construct metapolicies describing both the policy that they wish to enforce and annotations describing their composition preferences. These annotations can indicate whether certain policy assertions are required by the policy writer or, if not, under what circumstances the policy writer is willing to compromise and allow other assertions to take precedence. Meta-policies are specified in defeasible logic, a computationally efficient non-monotonic logic developed to model human reasoning. One drawback of this method is that at one point the algorithm starts an exhaustive search of all subsets of the set of conclusions of a defeasible theory. Although the propositional defeasible logic has linear complexity, the set of conclusions here may be large, especially in real-life practical cases. This phenomenon leads to an inefficient exponential explosion of complexity. The current process of getting a Web security policy from combination of two meta-policies consists of two steps. The first is generating a new meta-policy that is a composition of the input meta-policies, and the second is mapping the meta-policy onto a security policy. The new algorithm avoids the exhaustive search in the current algorithm, and provides a security policy that matches all requirements of the involved metapolicies.

Vatan, Farrokh↗

NASA NTRS2011

A New Inversion Routine to Produce Vertical Electron-Density Profiles from Ionospheric Topside-Sounder Data

Two software applications have been produced specifically for the analysis of some million digital topside ionograms produced by a recent analog-to-digital conversion effort of selected analog telemetry tapes from the Alouette-2, ISIS-1 and ISIS-2 satellites. One, TOPIST (TOPside Ionogram Scalar with True-height algorithm) from the University of Massachusetts Lowell, is designed for the automatic identification of the topside-ionogram ionospheric-reflection traces and their inversion into vertical electron-density profiles Ne(h). TOPIST also has the capability of manual intervention. The other application, from the Goddard Space Flight Center based on the FORTRAN code of John E. Jackson from the 1960s, is designed as an IDL-based interactive program for the scaling of selected digital topside-sounder ionograms. The Jackson code has also been modified, with some effort, so as to run on modern computers. This modification was motivated by the need to scale selected ionograms from the millions of Alouette/ISIS topside-sounder ionograms that only exist on 35-mm film. During this modification, it became evident that it would be more efficient to design a new code, based on the capabilities of present-day computers, than to continue to modify the old code. Such a new code has been produced and here we will describe its capabilities and compare Ne(h) profiles produced from it with those produced by the Jackson code. The concept of the new code is to assume an initial Ne(h) and derive a final Ne(h) through an iteration process that makes the resulting apparent-height profile fir the scaled values within a certain error range. The new code can be used on the X-, O-, and Z-mode traces. It does not assume any predefined profile shape between two contiguous points, like the exponential rule used in Jackson s program. Instead, Monotone Piecewise Cubic Interpolation is applied in the global profile to keep the monotone nature of the profile, which also ensures better smoothness in the final profile than in Jackson s program. The new code uses the complete refractive index expression for a cold collisionless plasma and can accommodate the IGRF, T96, and other geomagnetic field models.

Wang, Yongli↗

NASA NTRS2012

GX 3+1: The Stability of Spectral Index as a Function of Mass Accretion Rate

We present an analysis of the spectral and timing properties observed in X-rays from neutron star (NS) binary GX 3+1 (4U 1744-26) during long-term transitions between the faint and bright phases superimposed on short-term transitions between lower banana (LB) and upper banana (UB) branches in terms of its color-color diagram, We analyze all observations of this source obtained with the Rossi X-ray Timing Explorer and BeppoSAX satellites, We find that the X-ray broadband energy spectra during these spectral transitions can be adequately reproduced by a composition of a low-temperature blackbody component, a Comptonized component (COMPTB), and Gaussian component We argue that the electron temperature kTe of the Compton cloud monotonically increases from 2.3 keY to 4.5 keY, when GX 3+1 makes a transition from UB to LB. We also detect an evolution of noise components (a very low frequency noise and a high-frequency noise) during these LB-UB transitions. Using a disk seed photon normalization of COMPTB, which is proportional to the mass accretion rate, we find that the photon power-law index Gamma is almost constant (Gamma = 2.00 +/- 0.02) when mass accretion rate changes by factor four. In addition, we find that the emergent spectrum is dominated by the strong Comptonized component We interpret this quasi-stability of the index Gamma and a particular form of the spectrum in the framework of a model in which the energy release in the transition layer located between the accretion disk and NS surface dominates that in the disk. Moreover, this index stability effect now established for GX 3+ I was previously found in the atoll source 4U 1728-34 and suggested for a number of other low-mass X-ray NS binaries. This intrinsic behavior of NSs, in particular for atoll sources, is fundamentally different from that seen in black hole binary sources where the index monotonically increases during spectral transition from the low state to the high state and then finally saturates at high values of mass accretion rate.

Seifana, Elena↗

NASA NTRS2015

Using Single-Scattering Albedo Spectral Curvature to Characterize East Asian Aerosol Mixtures

Spectral dependence of aerosol single-scattering albedo (SSA) has been used to infer aerosol composition. In particular, aerosol mixtures dominated by dust absorption will have monotonically increasing SSA with wavelength while that dominated by black carbon absorption has monotonically decreasing SSA spectra. However, by analyzing SSA measured at four wavelengths, 440, 675, 870, and 1020 nm from the Aerosol Robotic Network data set, we find that the SSA spectra over East Asia are frequently peaked at 675 nm. In these cases, we suggest that SSA spectral curvature, defined as the negative of the second derivative of SSA as a function of wavelength, can provide additional information on the composition of these aerosol mixtures. Aerosol SSA spectral curvatures for East Asia during fall and winter are considerably larger than those found in places primarily dominated by biomass burning or dust aerosols. SSA curvature is found to increase as the SSA magnitude decreases. The curvature increases with coarse mode fraction (CMF) to a CMF value of about 0.4, then slightly decreases or remains constant at larger CMF. Mie calculations further verify that the strongest SSA curvature occurs at approx. 40% dust fraction, with 10% scattering aerosol fraction. The nonmonotonic SSA spectral dependence is likely associated with enhanced absorption in the shortwave by dust, absorption by black carbon at longer wavelengths, and also the flattened absorption optical depth spectral dependence due to the increased particle size.

curvature↗

NASA NTRS1964

Theoretical Current-Voltage Curve in Low-Pressure Cesium Diode for Electron-Rich Emission

Although considerable interest has been shown in the space-charge analysis of low-pressure (collisionless case) thermionic diodes, there is a conspicuous lack in the presentation of results in a way that allows direct comparison with experiment. The current-voltage curve of this report was, therefore, computed for a typical case within the realm of experimental interest. The model employed in this computation is shown in Fig. 1 and is defined by the limiting potential distributions [curves (a) and (b)]. Curve (a) represents the potential V as a monotonic function of position with a slope of zero at the anode; curve (b) is similarly monotonic with a slope of zero at the cathode. It is assumed that by a continuous variation of the anode voltage, the potential distributions vary continuously from one limiting form to the other. Although solutions for infinitely spaced electrodes show that spatically oscillatory potential distributions may exist, they have been neglected in this computation.

Coldstein, C. M.↗

NASA NTRS2017

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗

NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗

NASA NTRS2019

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a “Virtual Constellation” was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating “stable” calibration to within 5%—the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)–were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Student’s T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process.

Tuli, Fatima Tuz Zafrin↗

NASA NTRS2011

TPSAS-NF1676L-12354-DND

This work deals with performance properties of a dynamic traffic model, the Air Traffic Monotonic Lagrangian Grid (ATMLG), which can be used to evaluate new control strategies for conflict avoidance, separation assurance, and traffic management. The model is based on an algorithm and data structure called the Monotonic Lagrangian Grid (MLG), originally developed at NRL in the mid 1980s and since then used as an underpinning for various particle dynamics simulations. The MLG stores positions and other data needed to describe N moving objects, where N can be very large. The MLG algorithm involves sorting and ordering objects. A stationary grid is an alternative to the dynamic grid of MLG. Stationary grids can be attractive in that they do not require sorting. We investigate and report on the relative performances of air traffic simulations based on dynamic (MLG) and static (lat-long) grids.

C Kaplan↗

NASA NTRS2011

TPSAS-NF1676L-12301-DND

Carolyn Kaplan↗

NASA NTRS2022

Asymmetric Warming/Cooling Response to CO2 Increase/Decrease Mainly Due to Non-Logarithmic Forcing, not Feedbacks

We explore the CO2 dependence of effective climate sensitivity (SG) with symmetric abrupt and transient CO2 forcing, spanning the range 1/8×, 1/4×, 1/2×, 2×, 4×, and 8×CO2, using two state-of-the-art fully coupled atmosphere-ocean-sea-ice-land models. In both models, under abrupt CO2 forcing, we find an asymmetric response in surface temperature and SG. The surface global warming at 8×CO2 is more than one third larger than the corresponding cooling at 1/8×CO2, and SG is CO2 dependent, increasing non-monotonically from 1/8×CO2 to 8×CO2. We find similar CO2 dependence in the transient runs, forced with -1%yr−1CO2 and +1%yr−1CO2 up to 1/8×CO2 and 8×CO2, respectively. The non-logarithmic radiative forcing – not the changing feedbacks – primarily explains the dependence of SG on CO2, particularly at low CO2 levels. The changing feedbacks, however, explain SG’s non-monotonic behavior.

climate sensitivity↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

numerical methods↗

NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗

DOE OSTI2021

Resilient Operation of Power Distribution Systems Using MPC-Based Critical Service Restoration: Preprint

Power distribution systems are more prone to disruptions and cause most power system outages. We propose a service restoration technique to recover the system service (electricity delivery) following an extreme event-triggered substation outage. The proposed technique considers the problem of controlling distributed energy resources (DERs) of a distribution system with the objective of achieving maximum load pick up while satisfying network flow and voltage constraints. The problem is formulated as a model predictive control (MPC), where a linearized optimal power flow (OPF) model is employed to describe the network. The formulation is augmented with a ramping (up) reserve product for the DERs to ensure an upward monotonic load restoration as time evolves. We perform simulations considering the IEEE 13-bus test feeder integrated with wind, solar, microturbine, and energy storage devices. We demonstrate the efficacy of the devised technique in restoring the system loads monotonically, without shedding previously restored loads. We also show the benefit of co-optimization of power and reserve products for DERs on service restoration. In addition, the capability of the technique in regulating nodal voltages and reducing renewable power curtailment is demonstrated.

61 RADIATION PROTECTION AND DOSIMETRY↗