SEARCH · Engineering Papers

Results for “reactive flow”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Seasonal mixing from intermittent flow drives concentration-discharge behaviour in a stream affected by coal mine drainage

Abandoned mining operations continue to severely degrade many ecosystems worldwide by releasing acidic water and/or heavy metals into surface and groundwater. Contaminant concentrations in affected streams vary with discharge in patterns that reflect both geochemical reactions and variable mixing of contaminated and non-contaminated waters. However, controls on concentration-discharge (C-Q) patterns remain unclear, particularly for constituents that experience changing solubility across redox and pH gradients. Understanding the C-Q behaviour of contaminants aids in predicting both downstream transport and effects on aquatic life under variable flow. Here, we examined the C-Q behaviours of non-reactive (Na, K, Ca, Mg, Cl - ) and reactive (Fe, Mn, Al, H + , SO 4 2- ) solutes in a stream contaminated with acid mine drainage in northeastern Ohio, USA. Concentration-discharge patterns at the watershed outlet primarily reflected mixing of contaminated baseflow with intermittent inputs of high pH water draining from a passive limestone treatment system into the stream. The treatment system acted as an ephemeral tributary that mitigated contamination in the stream by diluting solutes, raising pH, and driving metal precipitation, but only when flow was present during wet seasons. Consequently, AMD-derived reactive solutes (H + , Fe, Mn, Al) decreased with increasing stream discharge while relatively conservative solutes (e.g., Ca, Mg, K, Na) decreased only slightly or were chemostatic. This study highlights both the unique C-Q patterns of reactive solutes when compared to those of non-reactive solutes and the potential for intermittent streams to control C-Q behaviour in headwater catchments.

58 GEOSCIENCES↗

DOE OSTI2023

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Visual Tool for Assessing Stability of DER Configurations on Three-Phase Radial Networks

Here, we present a method and tool for evaluating the placement of Distributed Energy Resources (DER) on distribution circuits in order to control voltages and power flows. Our previous work described Phasor-Based Control (PBC), a novel control framework where DERs inject real and reactive power to track voltage magnitude and phase angle targets. Here, we employ linearized power flow equations and integral controllers to develop a linear state space model for PBC acting on a three-phase unbalanced network. We use this model to evaluate whether a given inverter-based DER configuration admits a stable set of controller gains, which cannot be done by analyzing controllability nor by using the Lyapunov equation. Instead, we sample over a parameter space to identify a stable set of controller gains. Our stability analysis requires only a line impedance model and does not entail simulating the system or solving an optimization problem. We incorporate this assessment into a publicly available visualization tool and demonstrate three processes for evaluating many control configurations on the IEEE 123-node test feeder (123NF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Decentralized Control of Cascaded H-Bridge Inverters for Medium-Voltage Grid Integration: Preprint

This paper introduces a fully decentralized control strategy that allows for individualized active and reactive power control of H-bridge inverters connected in series. After formulating a model of the system power flows, we outline a communications-free approach where each inverter can act on independent commands. Inverter-level feedback is accomplished locally via a measurement of current flowing through the stack. A key feature is the lack of any centralized control routines or phase-locked loop. In this framework, the reactive power control loops achieve voltage synchronization among the converters and the grid, and the active power control loops modulate the terminal voltage amplitudes that influence power delivery. Furthermore, pole placement via state feedback is used to ensure stability of voltage synchronization while also enabling tracking of bidirectional power commands. Simulations are presented to validate the proposed approach.

communication-free modular converter control↗

DOE OSTI2023

Linearized Distribution Optimal Power Flow for OEDI SI

This research is to meant to demonstrate the OEDI SI use case for distributed optimal power flow (DOPF). The goal was to formulate the optimal power flow problem in the distribution system for active and reactive power setpoints of PV systems using topology information and voltage measurements. The co-simulation runs every 15 minutes as outlined within the scenario file for the given feeder configuration. The linked GitHub repository includes five federates to achieve DOPF for the small, medium, large, and IEEE 123 feeder scenarios. We are using the OEDI SI framework, as well as the example feeder, sensor, recorder, and estimator federates provided in the example repository for OEDI SI. We also provide a runner script for switching between scenarios.

algorithm↗

NASA NTRS1976

High performance N2O4/amine elements

An analytical and experimental investigation was conducted to develop an understanding of the mechanisms that cause reactive stream separation, commonly called blowapart, for hypergolic propellants. The investigation was limited to a N2O4/MMH propellant combination and to a range of engine-operating conditions applicable to the space tug and space shuttle attitude control and orbital maneuvering engines. Primary test variables were: chamber pressure (1 to 20 atm), fuel injection temperature (283 to 400 K)m and propellant injection velocity (9 to 50 m/s). The injector configuration studied was the unlike doublet. The reactive stream separation experiments were conducted using special combustors designed to permit photography of the near-injector spray combustion flow field. Analysis of color motion pictures provided the means of determining the occurrence of reactive stream separation.

Falk, A. Y.↗

NASA NTRS1976

High Performance N2O4/Amine Elements

An analytical and experimental investigation was conducted to develop an understanding of the mechanisms that cause reactive stream separation, commonly called blowapart, for hypergolic propellants. The investigation was limited to the N2O4/MMH propellant combination and to a range of engine-operating conditions applicable to the space tug and space shuttle attitude control and orbital maneuvering engines. Primary test variables were: chamber pressure (1 to 20 atm), fuel injection temperature (283 to 400 K)m and propellant injection velocity (9 to 50 m/s). The injector configuration studied was the unlike doublet. The reactive stream separation experiments were conducted using special combustors designed to permit photography of the near-injector spray combustion flow field. Analysis of color motion pictures provided the means of determining the occurrence of reactive stream separation.

Falk, A. Y.↗

DOE OSTI2022

Pele: An Exascale-Ready Suite of Combustion Codes

High fidelity simulations of realistic combustion devices are extremely demanding computationally because of the requirements to capture complex fuel chemical decomposition, its intricate interactions with turbulent, often multiphase, flows, and the wide separation of space and time scales between the thin flame and the device boundaries. Software required to carry out such computations tends to be extremely complex, particularly when designed to exploit hardware accelerators, and can be difficult to port and maintain. We present Pele, a performance portable suite of tools for the simulation of combustion systems, including codes to evolve reactive multiphase configurations in the low Mach number and compressible flow regimes, along with a set of inter-compatible post processing and in situ analysis tools. The Pele suite of tools is built on top of the AMReX framework for block-structured adaptive mesh refinement, which provides efficient data structures and algorithms that enable the development of a wide variety of efficient mesh and particle based PDE integration schemes. A hierarchical MPI+X parallelism scheme supports CPU-only and accelerated architectures, where X can be OpenMP, CUDA, and HIP based approaches for intra-node computational work distribution. The algorithms and data structures underlying the Pele simulation and analysis tools are highly scalable and performant across a wide variety of high-performance computing platforms, including DOEs newest exascale-class machines, Frontier and Aurora. The simulation and analysis tools are fully documented and freely distributed as open source via GitHub. We present key algorithmic and software challenges, solution strategies, performance and resulting set of capabilities.

AMReX↗

DOE OSTI2020

Etching of Si 3 N 4 induced by electron beam plasma from hollow cathode plasma in a downstream reactive environment

An etching system based on the interaction of electrons extracted from a direct current (DC) hollow-cathode (HC) Ar plasma and injected towards a Si3N4 covered silicon substrate located in the downstream reactive environment created by an additional remote CF 4 /O 2 plasma source was developed and evaluated. By controlling the properties of the injected beam electrons, this approach allows to deliver energy to a surface functionalized by exposure to reactive species, and initiate surface etching. The energy of the primary beam electrons is controlled by the acceleration voltage relative to the HC discharge. Ar atoms flow from the high-pressure HC discharge into the low pressure downstream reactive environment in the process chamber. For an acceleration voltage greater than the ionization potential of Ar and/or process gas species, the energetic primary beam electrons produce a secondary plasma in the process chamber, and can also cause additional dissociation. We have characterized the properties of the secondary plasma and also surface etching of Si 3 N 4 as a function of process parameters, including acceleration voltage (0-80 V), discharge current of the HC discharge (1-2 A), pressure (3.5-20 mTorr), source to substrate distance (1.5 to 5 cm), and feed gas composition (20% and 80% O2 in CF 4 /O 2 ). The electron energy probability function (EEPF) measured with a Langmuir probe about 2.5 cm below the extraction ring suggests several major groups of electrons for this situation, including high energy primary beam electrons with an energy that varies as the acceleration voltage is changed and low energy electrons produced by beam electron-induced ionization of the Ar gas in the process chamber. When a remote CF 4 /O 2 plasma is additionally coupled to the process chamber, Si 3 N 4 surfaces can be functionalized, and by varying the energy of the beam electrons, Si 3 N 4 etching can be induced by electron-neutral synergy effect with plasma-surface interaction. For conditions without beam electron injection, the remote plasma etching rate (ER) of Si 3 N 4 depends strongly on the O 2 concentration in the CF 4 /O 2 processing gas mixture and can be suppressed for O 2 -rich process conditions by the formation of a SiONF passivation layer on the Si 3 N 4 surface. The combination of the HC electron beam source with the remote plasma source (HCEB-RP) makes it possible to induce Si 3 N 4 etching for O 2 -rich remote plasma conditions where remote plasma by itself produces negligible Si 3 N 4 etching. The electron enhanced etching of Si 3 N 4 depends strongly on O 2 /CF 4 mixing ratio reflecting changing arrival rates of O and F species at the surface. Optical emission spectroscopy was used to estimate the ratio of gas phase F and O densities and found to be controlled by the gas mixing ratio and independent of HC EB operating conditions. At this time the detailed sequence of events operative in the etching mechanism is unclear. While the increase of the electron energy is ultimately responsible for initiating surface etching, at present we cannot rule out a role of ions from the simultaneously produced secondary plasma in plasma surface interaction mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Particle-based high-temperature thermochemical energy storage reactors

Solar and other renewable energy driven gas-solid thermochemical energy storage (TCES) technology is a promising solution for the next generation energy storage systems due to its high operating temperature, efficient energy conversion, ultra-long storage duration, and potential high energy density. Experimental and theoretical studies suggest that the respective gravimetric and volumetric TCES energy storage densities vary from 200 to 3000 kJ kg –1 and 1–3 GJ m –3 . Solar radiation or heat generated from electric furnaces powered by renewable electricity can be stored in the form of chemical energy through endothermic reactions, while the stored chemical energy can be converted to thermal energy via an exothermic reaction when needed. The design of highly effective reactors requires a deep understanding of materials, thermodynamics, chemical kinetics, and transport phenomena. At time of writing, TCES reactors are yet to be deployed at commercially relevant scales, leaving a substantial gap between development efforts and commercial feasibility. Therefore, this review aims to examine the state-of-the-art design and performance of particle-based TCES reactors with different reactive materials. Fundamentals related to TCES reactive materials, reaction conditions, thermodynamics and kinetics, and transport phenomena are reviewed in detail to provide a comprehensive understanding of the reactor design and operation. Five major types of TCES reactors have been comprehensively reviewed and compared, including fixed, moving, rotary, fluidized, and entrained bed reactors. Most reported prototype reactors in the literature operate at lab scale with thermal inputs below 40 kW, and scaled TCES reactors (e.g., at megawatt level) are yet to be demonstrated. The nominal reactor operating temperatures range from 300 to 1500 °C, depending on the selected chemistry, reactive material, and heat sources. To evaluate their designs, the reactors are assessed in aspects of performance, cost, and durability. Discrepancies in performance indicators of energy storage density, extent of reaction, and various energy efficiencies are highlighted. The scale-up of reactors and power block integration, which hold the key to the successful commercialization of TCES systems, are critically analyzed. Furthermore, advanced materials (both reactive materials and ceramic reactor housing materials), effective particle flow control, advanced modeling tools, and novel system design may bring significant improvement to the energy efficiency, storage density and cost competitiveness of particle-based TCES reactors.

25 ENERGY STORAGE↗

DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2001

Sedimentation along the Eastern Chenier Plain Coast: Down Drift Impact of a Delta Complex Shift

The Mississippi River Chenier Plain is a shore parallel landform (down-drift from the Atchafalaya distributary of the Mississippi River) consisting of an alternating series of transgressive sand-shell ridges and regressive, progradational mudflats. The late 1940s shift of 1/3 of the flow of the Mississippi to the newly developing Atchafalaya delta complex to the west has resulted in injection of the river waters and suspended sediment into the westward flowing currents of the coastal current system. This has reactivated the dormant processes of mud accumulation along this coast. These environmental circumstances have provided the opportunity to: (1) investigate the depositional processes of the prograding, fine grained, mud flat facies of the open Chenier main coast and (2) to test the hypothesis that the impacts of the frequent cold front passages of fall, winter and spring exceed those of the occasional and more localized hurricane in shaping the coast and powering the dominant sedimentary processes. We conducted field investigations with the benefit of multi - scale, time series environmental surveillance by remote sensing systems, including airborne and satellite sensors. These systems provided invaluable new information on areal geomorphic patterns and the behavior of the coastal waters. This is a classic case of weather impacting inner shelf waters and sediments and causing the development of a new landform. It is clear that mud flats of the eastern chenier plain are prograding seaward, as well as progressively growing in a westerly direction.

Huh, Oscar K.↗

NASA NTRSDate not supplied

Effect of Injection Dynamics on Wave Propagation in a Linear Detonation Combustor

Injector response to self-excited and sustained detonation wave propagation in an optically-accessible nonpremixed natural gas-oxygen linear detonation combustor is investigated. 100 kHz measurements of planar laser induced fluorescence (PLIF) of trace quantities of acetone injected into the fuel supply and simultaneous OH∗ chemiluminescence measurements were used to correlate the effect on injection dynamics on detonation wave structure and propagation. Measurements over a range of equivalence ratios (0.61-1.48) revealed wave propagation through regions of high fuel stratification in the axial and transverse direction of reactant injection. Consequently, the propagation of the detonation wave along the channel is through a flow-field with significant variation in fuel mass fraction, reactivity, and sound speeds. At lean conditions, wave propagation frequency is generally lower (∼6 kHz) resulting in longer time between wave passages for mixing and pre-heat resulting in steeper detonation wave-fronts anchored near the injector exit plane. Conversely, at rich equivalence ratios the wave frequency is higher (∼10 kHz) resulting in shorter intra-cycle period for reactant fill and mixing resulting waves lifted from the injection plane that are more compact. Higher chemiluminescence intensity in the post-wave region with deflagrative combustion is observed in these cases. Phase averaged representation of the acetone-PLIF measurements indicate that the fraction of time between wave passages available for reactant fill, mixing and pre-heat is consistent across all equivalence ratios.

Propulsion↗

NASA NTRSDate not supplied

Effect of Injection Dynamics on Detonation Wave Propagation in a Linear Detonation Combustor

Propulsion↗

DOE OSTI2020

A Combined Day-ahead and Real-time approach for Scheduling Real and Reactive power dispatch of Battery Energy Storage System

A combined day-ahead (look-ahead) and real-time optimization approach is proposed to schedule BESS real and reactive power dispatch. The look ahead scheduling agent maximizes profit and essentially performs the function of energy arbitrage by looking at forecasted prices, network constraints and BESS operating conditions to generate real power dispatch schedule. While this schedule is followed, the real time optimization agent computes an optimum reactive power set-point (for a particular objective function) using optimal power flow framework. The concepts are implemented by co-simulating legacy tools and demonstrated on a large real world distribution system. In addition to the profits from energy arbitrage, optimum reactive power set point computation results in achieving desirable network attributes like flat voltage profiles and minimizing voltage deviations from desired value.

Bhatti, Bilal A.↗