Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained control”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Enhanced Communication Network Solution for Positive Train Control Implementation

The commuter and freight railroad industry is required to implement Positive Train Control (PTC) by 2015 (2012 for Metrolink), a challenging network communications problem. This paper will discuss present technologies developed by the National Aeronautics and Space Administration (NASA) to overcome comparable communication challenges encountered in deep space mission operations. PTC will be based on a new cellular wireless packet Internet Protocol (IP) network. However, ensuring reliability in such a network is difficult due to the "dead zones" and transient disruptions we commonly experience when we lose calls in commercial cellular networks. These disruptions make it difficult to meet PTC s stringent reliability (99.999%) and safety requirements, deployment deadlines, and budget. This paper proposes innovative solutions based on space-proven technologies that would help meet these challenges: (1) Delay Tolerant Networking (DTN) technology, designed for use in resource-constrained, embedded systems and currently in use on the International Space Station, enables reliable communication over networks in which timely data acknowledgments might not be possible due to transient link outages. (2) Policy-Based Management (PBM) provides dynamic management capabilities, allowing vital data to be exchanged selectively (with priority) by utilizing alternative communication resources. The resulting network may help railroads implement PTC faster, cheaper, and more reliably.

Policy-Based Management (PBM)↗
DOE OSTI2020

Modeling Flexible Generator Operating Regions via Chance- constrained Stochastic Unit Commitment

Here, we introduce a novel chance-constrained stochastic unit commitment model to address uncertainty in renewables' production uncertainty in power systems operation. For most thermal generators,underlying technical constraints that are universally treated as "hard" by deterministic unit commitment models are in fact based on engineering judgments, such that system operators can periodically request operation outside these limits in non-nominal situations, e.g., to ensure reliability. We incorporate this practical consideration into a chance-constrained stochastic unit commitment model, specifically by in-frequently allowing minor deviations from the minimum and maximum thermal generator power output levels. We demonstrate that an extensive form of our model is computationally tractable for medium-sized power systems given modest numbers of scenarios for renewables' production. We show that the model is able to potentially save significant annual production costs by allowing infrequent and controlled violation of the traditionally hard bounds imposed on thermal generator production limits. Finally, we conduct a sensitivity analysis of optimal solutions to our model under two restricted regimes and observe similar qualitative results.

42 ENGINEERING↗
NASA NTRS1996

Non-Isothermal Experimental Study of the Constrained Vapor Bubble Thermosyphon

Experimental and theoretical techniques to study non-isothermal transport processes in the constrained vapor bubble thermosyphon (CVBT) were developed using a pentane/quartz system. The transport processes can be evaluated by measuring the liquid film profile, which gives the pressure field, and the temperature field. The axial variation in the capillary pressure was measured using an image-analyzing interferometer that is based on computer-enhanced video microscopy of the naturally occurring interference fringes. Thermoelectric coolers were used to control the temperature level in the condensation region and, therefore, the length of the approximately 'adiabatic' surface region which is a function of the temperature difference between the CVBT surface and the surroundings. High values for the axial thermal conductance in the 'adiabatic' surface region were demonstrated under certain conditions.

Karthikeyan, Muthu↗
NASA NTRS2001

Establishing the Dependence of [HO2]/[OH] on Temperature, Halogen Loading, O3, and NO(x) Based on in Situ Measurements from the NASA ER-2

In situ observations of OH and HO2 from the Airborne Southern Hemisphere Ozone Experiment/Measurements for Assessing the Effects of Stratospheric Aircraft (ASHOE/MAESA), Stratospheric TRacers of Atmospheric Transport (STRAT), and Polar Ozone Loss in the Arctic Region in Summer (POLARIS) NASA ER-2 field campaigns are used to examine the partitioning of HO(x) in the lower stratosphere (tropopause to approx.21 km) and upper troposphere (approx.10 km to tropopause). These measurements span a latitude range from 70degS to 90degN and a variety of atmospheric conditions as a result of seasonal changes and altitude. The response of the observed [HO2]/[OH] to changes in temperature, [03], [CO], [NO], [CIO], and [BrO] is investigated. The measured ratio is accurately described (approx.+/-10%) by a steady-state model constrained by the measured mixing ratios of O3, CO, NO, CIO, and BrO, where the model is valid for conditions of HO(x) cycling much faster than HO(x) production and loss. The concentration of HO2 depends on [OH], which, to first order, has been observed to be a simple function of the solar zenith angle in the lower stratosphere. The partitioning between OH and HO2 is controlled by the local chemistry between the HO, radicals and O3, CO, NO, CIO, and BrO. The response of [HO(x)] to changes in [NO(x)] and [O3] is demonstrated. Further observations are necessary to illustrate the response of HO(x) to changes in halogen concentrations. A quantitative understanding of [HO2]/[OH] is important, since many of the reactions that control this ratio are directly involved in catalytic removal of O3 in the lower stratosphere and production of O3 in the upper troposphere.

Lanzendorf, E. J.↗
NASA NTRSDate not supplied

Surface Transportation of the Common Habitat from Lander to Habitation Zone

The Common Habitat Architecture is a feasibility study surrounding the use of an SLS core stage liquid oxygen tank as the pressure vessel for a long-duration habitat intended for use in multiple gravity environments. The Common Habitat is used within this study as the primary habitation element in both Moon and Mars surface base camps. The Common Habitat Architecture offloads the Common Habitat from its lander and transports it to a Habitation Zone instead of leaving the habitat integrated with the lander. In this study, the Landing Zone is assumed to be approximately 3.5 kilometers from the Habitation Zone. Given the physical size and estimated 90-ton mass of the Common Habitat, a four-week trade study encompassing Moon and Mars lander identification, offloading, surface transport, and emplacement was conducted in February 2021 to assess whether there are any credible options for landing the Common Habitat on the Moon or Mars and deliver it to its intended point of use. Constrained to use only public data, the study applied subject matter expert opinion to each component of the study. The surface transportation component of the trade study assumes the habitat has been successfully offloaded from its lander and is at a point of handover to the surface transportation system. It is assumed that there are no crew present, and all human operations are performed remotely by Mission Control personnel. Three core cargo handling elements from prior NASA studies were used as the basis from which to derive surface transportation options – the Chariot, the All-Terrain Hex-Limbed Extra-Terrestrial Explorer, and the Lightweight Surface Manipulator System. Variations and hybrid combinations of these elements were used to develop transportation options for the Moon and Mars, given different habitat masses. Ultimately, several potentially feasible solutions were identified, and a solution was recommended that is common to both the Moon and Mars, with the potential to use a dissimilar system as a backup during lunar Common Habitat delivery. Next steps for sizing and additional development and analysis of the recommended surface transportation system are included as forward work.

Surface Mobility↗
DOE OSTI2025

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗
DOE OSTI2023

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

The solubility of NdPO 4 and DyPO 4 and stability of Nd and Dy chloride and hydroxyl complexes as a function of pH and salinity to 450 °C

The mobility of rare earth elements (REE) in geological systems is often controlled by the stability of monazite and xenotime, which are important hosts for the light and heavy REE. These REE phosphates constitute important resources and provide information on ore formation conditions and timing due to their uses as geothermometers and geochronometers. While the thermodynamic properties for these minerals are well-established up to high temperature and pressure, the properties for the aqueous REE complexes are not well constrained with limited information up to 300 °C. In this study, the solubility of synthetic NdPO 4 and DyPO 4 endmembers were measured in hydrothermal sub- to supercritical NaCl-bearing aqueous solutions from 350 at saturated water vapor pressure to 450 °C at 700 bar. The speciation of Nd and Dy was determined in acidic to alkaline solutions (pH 25 °C values from 2 to 10), and indicates that the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 predominate at low salinity (0.01 mol/kg NaCl) with some contribution of REE chloride species REECl 2+ and REECl 2 + in acidic fluids. The REE phosphate solubility measured in these experiments is up to two orders of magnitude higher than predicted using existing thermodynamic properties from literature for aqueous species extrapolated from lower temperature data to supercritical conditions. To address this discrepancy, the thermodynamic properties of the REE aqueous species were optimized using GEMSFITS to derive the formation constants for the REE chloride (β Cl ) and hydroxyl (β OH ) complexes in the studied temperature range. This study highlights the importance of the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 over a wide range of pH and temperature. As a result, the revised thermodynamic properties for Dy and Nd chloride and hydroxyl species more accurately predict the solubility of NdPO 4 and DyPO 4 and provide insight into the speciation and mobility of the light and heavy REE in hydrothermal supercritical fluids in the crust.

58 GEOSCIENCES↗
DOE OSTI2023

Bark beetle impacts on forest evapotranspiration and its partitioning

Insect outbreaks affect forest structure and function and represent a major category of forest disturbance globally. However, the resulting impacts on evapotranspiration (ET), and especially hydrological partitioning between the abiotic (evaporation) and biotic (transpiration) components of total ET, are not well constrained. As a result, we combined remote sensing, eddy covariance, and hydrological modeling approaches to determine the effects of bark beetle outbreak on ET and its partitioning at multiple scales throughout the Southern Rocky Mountain Ecoregion (SRME), USA. At the eddy covariance measurement scale, 85 % of the forest was affected by beetles, and water year ET as a fraction of precipitation (P) decreased by 30 % relative to a control site, with 31 % greater reductions in growing season transpiration relative to total ET. At the ecoregion scale, satellite remote sensing masked to areas of >80 % tree mortality showed corresponding ET/P reductions of 9–15 % that occurred 6–8 years post-disturbance, and indicated that the majority of the total reduction occurred during the growing season; the Variable Infiltration Capacity hydrological model showed an associated 9–18 % increase in the ecoregion runoff ratio. Long-term (16–18 year) ET and vegetation mortality datasets extend the length of previously published analyses and allowed for clear characterization of the forest recovery period. During that time, transpiration recovery outpaced total ET recovery, which was lagged in part due to persistently reduced winter sublimation, and there was associated evidence of increasing late summer vegetation moisture stress. Overall, comparison of three independent methods and two partitioning approaches demonstrated a net negative impact of bark beetles on ET, and a relatively greater negative impact on transpiration, following bark beetle outbreak in the SRME.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Experimental fire increases soil carbon dioxide efflux in a grassland long‐term multifactor global change experiment

Abstract Numerous studies have demonstrated that soil respiration rates increase under experimental warming, although the long‐term, multiyear dynamics of this feedback are not well constrained. Less is known about the effects of single, punctuated events in combination with other longer‐duration anthropogenic influences on the dynamics of soil carbon (C) loss. In 2012 and 2013, we assessed the effects of decadal‐scale anthropogenic global change – warming, increased nitrogen (N) deposition, elevated carbon dioxide ( CO 2 ), and increased precipitation – on soil respiration rates in an annual‐dominated Mediterranean grassland. We also investigated how controlled fire and an artificial wet‐up event, in combination with exposure to the longer‐duration anthropogenic global change factors, influenced the dynamics of C cycling in this system. Decade‐duration surface soil warming (1–2 °C) had no effect on soil respiration rates, while +N addition and elevated CO 2 concentrations increased growing‐season soil CO 2 efflux rates by increasing annual aboveground net primary production ( NPP ) and belowground fine root production, respectively. Low‐intensity experimental fire significantly elevated soil CO 2 efflux rates in the next growing season. Based on mixed‐effects modeling and structural equation modeling, low‐intensity fire increased growing‐season soil respiration rates through a combination of three mechanisms: large increases in soil temperature (3–5 °C), significant increases in fine root production, and elevated aboveground NPP . Our study shows that in ecosystems where soil respiration has acclimated to moderate warming, further increases in soil temperature can stimulate greater soil CO 2 efflux. We also demonstrate that punctuated short‐duration events such as fire can influence soil C dynamics with implications for both the parameterization of earth system models ( ESM s) and the implementation of climate change mitigation policies that involve land‐sector C accounting.

Strong, Aaron L.↗
NASA NTRS1984

Stagger angle dependence of inertial and elastic coupling in bladed disks

Conditions which necessitate the inclusion of disk and shaft flexibility in the analysis of blade response in rotating blade-disk-shaft systems are derived in terms of nondimensional parameters. A simple semianalytical Rayleigh-Ritz model is derived in which the disk possesses all six rigid body degrees of freedom, which are elastically constrained by the shaft. Inertial coupling by the rigid body motion of the disk on a flexible shaft and out-of-plane elastic coupling due to disk flexure are included. Frequency ratios and mass ratios, which depend on the stagger angle, are determined for three typical rotors: a first stage high-pressure core compressor, a high bypass ratio fan, and an advanced turboprop. The stagger angle controls the degree of coupling in the blade-disk system. In the blade-disk-shaft system, the stagger angle determines whether blade-disk motion couples principally to the out-of-plane or in-plane motion of the disk on the shaft. The Ritz analysis shows excellent agreement with experimental results.

Crawley, E. F.↗
NASA NTRS1991

Solar Pumped Lasers and Their Applications

Since 1980, NASA has been pursuing high power solar lasers as part of the space power beaming program. Materials in liquid, solid, and gas phases have been evaluated against the requirements for solar pumping. Two basic characteristics of solar insolation, namely its diffuse irradiance and 5800 K blackbody-like spectrum, impose rather stringent requirements for laser excitation. However, meeting these requirements is not insurmountable as solar thermal energy technology has progressed today, and taking advantage of solar pumping lasers is becoming increasingly attractive. The high density photons of concentrated solar energy have been used for mainly electric power generation and thermal processing of materials by the DOE Solar Thermal Technologies Program. However, the photons can interact with materials through many other direct kinetic paths, and applications of the concentrated photons could be extended to processes requiring photolysis, photosynthesis, and photoexcitation. The use of solar pumped lasers on Earth seems constrained by economics and sociopolitics. Therefore, prospective applications may be limited to those that require use of quantum effects and coherency of the laser in order to generate extremely high value products and services when conventional and inexpensive means are ineffective or impossible. The new applications already proposed for concentrated solar photons, such as destruction of hazardous waste, production of renewable fuel, production of fertilizer, and air/water pollution controls, may benefit from the use of inexpensive solar pumped laser matched with the photochemical kinetics of these processes.

Lee, Ja H.↗
NASA NTRS2010

Lessons Learned in Launching Ares I-X

The Ares I-X flight test, scheduled for 2009, is the first opportunity for the Constellation Program and Ares Projects to obtain important data on the in-flight loads, first stage recovery, and ground-handling characteristics of the Ares I crew launch vehicle. The flight test vehicle will incorporate a mix of flight and mockup hardware, reflecting an acceptable representation of the mass and outer mold line characteristics of the operational Ares I vehicle. It will be powered by a four-segment solid rocket motor from the Space Shuttle inventory and will include a fifth, inert spacer segment and new forward structures to make it the same shape as the Ares I first stage. The vehicle also includes mass simulators for the upper stage, Orion crew module, and launch abort system to match the outer mold line of Ares I; an active roll control system; and avionics derived from Shuttle and Atlas hardware. This suborbital mission will take the flight test vehicle from stacking in the Vehicle Assembly Building (VAB) at Kennedy Space Center (KSC) through liftoff, maximum dynamic pressure, first stage separation, and recovery. Ares I-X presented NASA with unique project management challenges. This presentation will address those unique challenges, including managing a virtual nationwide team under a constrained timeline; designing and building the launch vehicle and ground systems; working concurrently with Space Shuttle activities; and integrating technical and management functions. It also will identify key data collected and lessons learned from the flight that will be applied to future exploration missions.

Davis, Stephan R.↗
DOE OSTI2012

Sensor placement algorithm development to maximize the efficiency of acid gas removal unit for integrated gasifiction combined sycle (IGCC) power plant with CO2 capture

Future integrated gasification combined cycle (IGCC) power plants with CO{sub 2} capture will face stricter operational and environmental constraints. Accurate values of relevant states/outputs/disturbances are needed to satisfy these constraints and to maximize the operational efficiency. Unfortunately, a number of these process variables cannot be measured while a number of them can be measured, but have low precision, reliability, or signal-to-noise ratio. In this work, a sensor placement (SP) algorithm is developed for optimal selection of sensor location, number, and type that can maximize the plant efficiency and result in a desired precision of the relevant measured/unmeasured states. In this work, an SP algorithm is developed for an selective, dual-stage Selexol-based acid gas removal (AGR) unit for an IGCC plant with pre-combustion CO{sub 2} capture. A comprehensive nonlinear dynamic model of the AGR unit is developed in Aspen Plus Dynamics® (APD) and used to generate a linear state-space model that is used in the SP algorithm. The SP algorithm is developed with the assumption that an optimal Kalman filter will be implemented in the plant for state and disturbance estimation. The algorithm is developed assuming steady-state Kalman filtering and steady-state operation of the plant. The control system is considered to operate based on the estimated states and thereby, captures the effects of the SP algorithm on the overall plant efficiency. The optimization problem is solved by Genetic Algorithm (GA) considering both linear and nonlinear equality and inequality constraints. Due to the very large number of candidate sets available for sensor placement and because of the long time that it takes to solve the constrained optimization problem that includes more than 1000 states, solution of this problem is computationally expensive. For reducing the computation time, parallel computing is performed using the Distributed Computing Server (DCS®) and the Parallel Computing® toolbox from Mathworks®. In this presentation, we will share our experience in setting up parallel computing using GA in the MATLAB® environment and present the overall approach for achieving higher computational efficiency in this framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Sensor placement algorithm development to maximize the efficiency of acid gas removal unit for integrated gasification combined cycle (IGCC) power plant with CO{sub 2} capture

Future integrated gasification combined cycle (IGCC) power plants with CO{sub 2} capture will face stricter operational and environmental constraints. Accurate values of relevant states/outputs/disturbances are needed to satisfy these constraints and to maximize the operational efficiency. Unfortunately, a number of these process variables cannot be measured while a number of them can be measured, but have low precision, reliability, or signal-to-noise ratio. In this work, a sensor placement (SP) algorithm is developed for optimal selection of sensor location, number, and type that can maximize the plant efficiency and result in a desired precision of the relevant measured/unmeasured states. In this work, an SP algorithm is developed for an selective, dual-stage Selexol-based acid gas removal (AGR) unit for an IGCC plant with pre-combustion CO{sub 2} capture. A comprehensive nonlinear dynamic model of the AGR unit is developed in Aspen Plus Dynamics® (APD) and used to generate a linear state-space model that is used in the SP algorithm. The SP algorithm is developed with the assumption that an optimal Kalman filter will be implemented in the plant for state and disturbance estimation. The algorithm is developed assuming steady-state Kalman filtering and steady-state operation of the plant. The control system is considered to operate based on the estimated states and thereby, captures the effects of the SP algorithm on the overall plant efficiency. The optimization problem is solved by Genetic Algorithm (GA) considering both linear and nonlinear equality and inequality constraints. Due to the very large number of candidate sets available for sensor placement and because of the long time that it takes to solve the constrained optimization problem that includes more than 1000 states, solution of this problem is computationally expensive. For reducing the computation time, parallel computing is performed using the Distributed Computing Server (DCS®) and the Parallel Computing® toolbox from Mathworks®. In this presentation, we will share our experience in setting up parallel computing using GA in the MATLAB® environment and present the overall approach for achieving higher computational efficiency in this framework.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2016

Advanced Energy Storage Management in Distribution Network

With increasing penetration of distributed generation (DG) in the distribution networks (DN), the secure and optimal operation of DN has become an important concern. In this paper, an iterative mixed integer quadratic constrained quadratic programming model to optimize the operation of a three phase unbalanced distribution system with high penetration of Photovoltaic (PV) panels, DG and energy storage (ES) is developed. The proposed model minimizes not only the operating cost, including fuel cost and purchasing cost, but also voltage deviations and power loss. The optimization model is based on the linearized sensitivity coefficients between state variables (e.g., node voltages) and control variables (e.g., real and reactive power injections of DG and ES). To avoid slow convergence when close to the optimum, a golden search method is introduced to control the step size and accelerate the convergence. The proposed algorithm is demonstrated on modified IEEE 13 nodes test feeders with multiple PV panels, DG and ES. Numerical simulation results validate the proposed algorithm. Various scenarios of system configuration are studied and some critical findings are concluded.

linearization↗