Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HELICAL 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.

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

At least 307 records · Page 17

Large-Scale Flow Structure in Turbulent Nonpremixed Flames under Normal- And Low-Gravity Conditions

It is well known that buoyancy has a major influence on the flow structure of turbulent nonpremixed jet flames. Buoyancy acts by inducing baroclinic torques, which generate large-scale vortical structures that can significantly modify the flow field. Furthermore, some suggest that buoyancy can substantially influence the large-scale structure of even nominally momentum-dominated flames, since the low velocity flow outside of the flame will be more susceptible to buoyancy effects. Even subtle buoyancy effects may be important because changes in the large-scale structure affects the local entrainment and fluctuating strain rate, and hence the structure of the flame. Previous studies that have compared the structure of normal- and micro-gravity nonpremixed jet flames note that flames in microgravity are longer and wider than in normal-gravity. This trend was observed for jet flames ranging from laminar to turbulent regimes. Furthermore, imaging of the flames has shown possible evidence of helical instabilities and disturbances starting from the base of the flame in microgravity. In contrast, these characteristics were not observed in normal-gravity. The objective of the present study is to further advance our knowledge of the effects of weak levels of buoyancy on the structure of transitional and turbulent nonpremixed jet flames. In later studies we will utilize the drop tower facilities at NASA Glenn Research Center (GRC), but the preliminary work described in this paper was conducted using the 1.25-second drop tower located at the University of Texas at Austin. A more detailed description of these experiments can be found in Idicheria et al.

Clemens, N. T.↗

Extended magnetohydrodynamics simulations of thin-foil Z-pinch implosions with comparison to experiments

Cylindrical foil liners, with foil thicknesses on the order of 400 nm, are often used in university-scale Z-pinch experiments (~1 MA in 100 ns) to study physics relevant to inertial confinement fusion efforts on larger-scale facilities (e.g., the MagLIF effort on the 25-MA Z facility at Sandia National Laboratories). The use of ultrathin foil liners typically requires a central support rod to maintain the structural integrity of the liner target assembly prior to implosion. The radius of this support rod sets a limit on the maximum convergence ratio achievable for the implosion. In recent experiments with a support rod and a preimposed axial magnetic field, helical instability structures in the imploding foil plasma were found to persist as the foil plasma stagnated on the rod and subsequently expanded away from the rod. We have now used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) to study these experiments. The results suggest that it is the support rod which is responsible for the helical structures persisting beyond stagnation. Furthermore, we find that as the radius of the support rod decreases (i.e., as the convergence ratio increases), the integrity and persistence of the helical modes diminish. In the limit with no support rod, we find that the structure of the final stagnation column is governed by the structure of the central precursor plasma column. Furthermore, these simulation results and their comparisons to experiment are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent boundary layer on a rotating helical blade.

Investigation of the boundary-layer characteristics on a helical blade of large chord length, enclosed in an annulus and rotating in a fluid otherwise at rest. The three-dimensional form of momentum integral equations is derived, and is used to predict the boundary-layer growth and limiting streamline angles on the blade surface. The measurements are in general agreement with the predictions. The wall shear stress correlation, which includes both Reynolds number and rotation parameters, valid for a rotating blade operating at zero pressure gradient, is derived. Radial and tangential velocity profiles, the tangential component of turbulence intensity, and blade static pressures are measured at several locations on the blade surface. The nature of flow near the blade tip is discussed. An expression for the radial velocity profile, valid in the outer region of the boundary layer, is derived theoretically.

Lakshminarayana, B.↗

An Evaluation of Linear Instability Waves as Sources of Sound in a Supersonic Turbulent Jet

Mach wave radiation from supersonic jets is revisited to better quantify the extent to which linearized equations represent the details of the actual mechanism. To this end, we solve the linearized Navier-Stokes equations (LNS) with precisely the same mean flow and inflow disturbances as a previous direct numerical simulation (DNS) of a perfectly expanded turbulent M = 1.92 jet. We restrict our attention to the first two azimuthal modes, n = 0 and n = 1, which constitute most of the acoustic field. The direction of peak radiation and the peak Strouhal number matches the DNS reasonably well, which is in accord with previous experimental justification of the linear theory. However, it is found that the sound pressure level predicted by LNS is significantly lower than that from DNS. In order to investigate the discrepancy, individual frequency components of the solution are examined. These confirm that near the peak Strouhal number, particularly for the first helical mode n = 1, the amplification of disturbances in the LNS closely matches the DNS. However, away from the peak frequency (and generally for the azimuthal mode n = 0), modes in the LNS are damped while those in the DNS grow at rates comparable to those at the peak Strouhal number.

Mohseni, Kamran↗

Shear-Sensitive Liquid Crystal Coating Method: Surface-Inclination Effects on Shear Vector Measurements

The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Reda, Daniel C.↗

Fast-Running Aeroelastic Code Based on Unsteady Linearized Aerodynamic Solver Developed

The NASA Glenn Research Center has been developing aeroelastic analyses for turbomachines for use by NASA and industry. An aeroelastic analysis consists of a structural dynamic model, an unsteady aerodynamic model, and a procedure to couple the two models. The structural models are well developed. Hence, most of the development for the aeroelastic analysis of turbomachines has involved adapting and using unsteady aerodynamic models. Two methods are used in developing unsteady aerodynamic analysis procedures for the flutter and forced response of turbomachines: (1) the time domain method and (2) the frequency domain method. Codes based on time domain methods require considerable computational time and, hence, cannot be used during the design process. Frequency domain methods eliminate the time dependence by assuming harmonic motion and, hence, require less computational time. Early frequency domain analyses methods neglected the important physics of steady loading on the analyses for simplicity. A fast-running unsteady aerodynamic code, LINFLUX, which includes steady loading and is based on the frequency domain method, has been modified for flutter and response calculations. LINFLUX, solves unsteady linearized Euler equations for calculating the unsteady aerodynamic forces on the blades, starting from a steady nonlinear aerodynamic solution. First, we obtained a steady aerodynamic solution for a given flow condition using the nonlinear unsteady aerodynamic code TURBO. A blade vibration analysis was done to determine the frequencies and mode shapes of the vibrating blades, and an interface code was used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor was used to interpolate the mode shapes from the structural dynamic mesh onto the computational dynamics mesh. Then, we used LINFLUX to calculate the unsteady aerodynamic forces for a given mode, frequency, and phase angle. A postprocessor read these unsteady pressures and calculated the generalized aerodynamic forces, eigenvalues, and response amplitudes. The eigenvalues determine the flutter frequency and damping. As a test case, the flutter of a helical fan was calculated with LINFLUX and compared with calculations from TURBO-AE, a nonlinear time domain code, and from ASTROP2, a code based on linear unsteady aerodynamics.

Reddy, T. S. R.↗

A Model for Energy Buildup and Eruption Onset in Coronal Mass Ejections

Coronal mass ejections (CMEs) and eruptive flares (EFs) are the most energetic explosions in the solar system. Their underlying origin is the free energy that builds up slowly in the sheared magnetic field of a filament channel. We report the first end-to-end numerical simulation of a CME/EF, from zero-free-energy initial state through filament channel formation to violent eruption, driven solely by the magnetic-helicity condensation process. Helicity is the topological measure of linkages between magnetic flux systems, and is conserved in the corona, building up inexorably until it is ejected into interplanetary space. Numerous investigations have demonstrated that helicity injected by small-scale vortical motions, such as those observed in the photosphere, undergoes an inverse cascade from small scales to large, “condensing” at magnetic-polarity boundaries. Our new results verify that this process forms a filament channel within a compact bipolar region embedded in a background dipole field, and show for the first time that a fast CME eventually occurs via the magnetic-breakout mechanism. We further show that the trigger for explosive eruption is reconnection onset in the flare current sheet that develops above the polarity inversion line: this reconnection forms flare loops below the sheet and a CME flux rope above, and initiates high-speed outward flow of the CME. Our findings have important implications for magnetic self organization and explosive behavior in solar and other astrophysical plasmas, as well as for understanding and predicting explosive solar activity.

filaments↗

Unsteady wing surface pressures in the wake of a propeller

The unsteady nature of the propeller slipstream interacting with a wing has been studied by flow visualization and unsteady wing surface pressure measurements. Flow visualization was performed by marking the propeller tip vortex with smoke. Unsteady wing surface pressures were measured by traversing a wing instrumented with a chordwise array of 16 microphones in a spanwise direction through the propeller wake. This work yielded information on the motion of the propeller wake as it passes over the wing. As the propeller wake passed over the wing: the propeller tip vortex experienced an inviscid interaction at the leading edge; viscous action at the leading edge severed the propeller tip vortex; the propeller tip vortex experienced significant spanwise and chordwise displacements and then deformed in order to reconnect at the trailing edge; axial velocity in the vortex core caused the helical vortex to thicken or stretch near the wing surface; and, the magnitude of the pressure fluctuations decreased in magnitude with distance traveled along the chord.

Johnston, R. T.↗

High Power Light Gas Helicon Plasma Source For VASMIR

The VASIMR space propulsion development effort relies on a high power (greater than 10kW) helicon source to produce a dense flowing plasma (H, D and He) target for ion cyclotron resonance (ICR) acceleration of the ions. Subsequent expansion in an expanding magnetic field (magnetic nozzle) converts ion lunetic energy to directed momentum. This plasma source must have critical features to enable an effective propulsion device. First, it must ionize most of the input neutral flux of gas, thus producing a plasma stream with a high degree of ionization for application of ICR power. This avoids propellant waste and potential power losses due to charge exchange. Next, the plasma stream must flow into a region of high magnetic field (approximately 0.5 T) for efficient ICR acceleration. Third, the ratio of input power to plasma flux must be low, providing an energy per ion-electron pair approaching 100 eV. Lastly, the source must be robust and capable of very long life-times (years). In our helicon experiment (VX-10) we have measured a ratio of input gas to plasma flux near 100%. The plasma flows from the helicon region (B approximately 0.1 T) into a region with a peak magnetic field of 0.8 T. The energy input per ion-electron pair has been measured at 300 plus or minus 100 eV. Recent results at Oak Ridge National Laboratory (ORNL) show an enhanced efficiency mode of operation with a high power density, over 5 kW in a 5 cm diameter tube. Our helicon is presently 9 cm in diameter and operates up to 3.5 kW of input power. An upgrade to a power level of 10 kW is underway. Much of our recent work has been with a Boswell double-saddle antenna design. We are also converting the antenna design to a helical type. With these modifications, we anticipate an improvement in the ionization efficiency. This paper presents the results from scaling the helicon in the VX-10 device from 3.5 to 10 kW. We also compare the operation with a double-saddle to a helical antenna design. Finally, we discuss modeling of these configurations using ORNL's EMIR code.

Squire, J. P.↗

Creep and creep rupture of strongly reinforced metallic composites

A creep and creep damage theory is presented for metallic composites with strong fibers. Application is to reinforced structures in which the fiber orientation may vary throughout but a distinct fiber direction can be identified locally (local transverse isotropy). The creep deformation model follows earlier work and is based on a flow potential function that depends on invariants reflecting stress and the material symmetry. As the focus is on the interaction of creep and damage, primary creep is ignored. The creep rupture model is an extension of continuum damage mechanics and includes an isochronous damage function that depends on invariants specifying the local maximum transverse tension and the maximum longitudinal shear stress. It is posited that at high temperature and low stress, appropriate to engineering practice, these stress components damage the fiber/matrix interface through diffusion controlled void growth, eventually causing creep rupture. Experiments are outlined for characterizing a composite through creep rupture tests under transverse tension and longitudinal shear. Application is made to a thin-walled pressure vessel with reinforcing fibers at an arbitrary helical angle. The results illustrate the usefulness of the model as a means of achieving optimal designs of composite structures where creep and creep rupture are life limiting.

Robinson, D. N.↗

Optimization-Inspired Pin-Fin Array for Supercritical Carbon Dioxide Recuperator

Additively manufactured heat exchangers are one possible route to cost-effective sCO 2 power cycles. In this paper, experimental results are obtained for two helical pin fin tubes that were designed following parametric optimization of the fin array. Neither the numerical optimization nor the experimental testing have been previously reported in the literature. The two tube designs were, (1) the optimization inspired design, and (2) the optimization-inspired design with a fin diameter increased by a factor of 2. To characterize the print, the designs were scanned using X-ray computed tomography to measure feature sizes and heat transfer area. The optimization was conducted in a commercial, computational fluid dynamics code. The code solved the Reynolds Averaged Navier-Stokes (RANS) and energy equations with turbulence closure provided by the shear stress transport (SST) k-ω model. In the experiments, the Nusselt number augmentation was measured using the Wilson plot technique and the friction factor was determined with mass flow and pressure drop measurements. The experimental testing indicated that the optimization-inspired design had a friction factor that was four times less than the baseline tube design at equal Nusselt number. Additionally, the optimization-inspired design had a 14% improvement in Nusselt number, at equal friction factor, relative to the best performing computational fluid dynamics (CFD) trial points. Tube design (2), with the larger diameter pin fins, had similar performance, within experimental error, as the tube with the smaller diameter pin fins (1). Both tubes achieved overall fin array efficiencies near 1. A performance factor, V/V 0 , equal to the volume of the enhanced heat exchanger divided by the volume of the baseline (no-fins) heat exchanger, is recommended to quantify internal cooling performance. In conclusion, the experimental shell and tube heat exchanger, using the additively manufactured tube, is competitive with printed circuit heat exchangers in its pressure drop class and could be further improved by optimizing a shell-and-tube heat exchanger utilizing this heat transfer enhancement feature.

42 ENGINEERING↗

Ring Bubbles of Dolphins

The article discusses how dolphins create and play with three types of air-filled vortices. The underlying physics is discussed. Photographs and sketches illustrating the dolphin's actions and physics are presented. The dolphins engage in this behavior on their own initiative without food reward. These behaviors are done repeatedly and with singleminded effort. The first type is the ejection of bubbles which, after some practice on the part of the dolphin, turn into toroidal vortex ring bubbles by the mechanism of baroclinic torque. These bubbles grow in radius and become thinner as they rise vertically to the surface. One dolphin would blow two in succession and guide them to fuse into one. Physicists call this a vortex reconnection. In the second type, the dolphins first create an invisible vortex ring in the water by swimming on their side and waving their tail fin (also called flukes) vigorously. This vortex ring travels horizontally in the water. The dolphin then turns around, finds the vortex and injects a stream of air into it from its blowhole. The air "fills-out" the core of the vortex ring. Often, the dolphin would knock-off a smaller ring bubble from the larger ring (this also involves vortex reconnection) and steer the smaller ring around the tank. One other dolphin employed a few other techniques for planting air into the fluke vortex. One technique included standing vertically in the water with tail-up, head-down and tail piercing the free surface. As the fluke is waved to create the vortex ring, air is entrained from above the surface. Another technique was gulping air in the mouth, diving down, releasing air bubbles from the mouth and curling them into a ring when they rose to the level of the fluke. In the third type, demonstrated by only one dolphin, the longitudinal vortex created by the dorsal fin on the back is used to produce 10-15 foot long helical bubbles. In one technique she swims in a curved path. This creates a dorsal fin vortex since centrifugal force has to be balanced by a lift-like force. She then re-traces her path and injects air into the vortex from her blowhole. She can even make a ring reconnect from the helix. In the second technique, demonstrated a few times, she again swims in a curved path, releases a cloud or group of bubbles from her blowhole and turns sharply away (Which presumably strengthens the vortex). As the bubbles encounter the vortex, they travel to the center of the vortex, merge and, in a flash, elongate along the core of the vortex. In all the three types, the air-water interface is shiny smooth and stable because the pressure gradient in the vortex flow around the bubble stabilizes it. A lot of the interesting physics still remains to be explored.

Shariff, Karim↗

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of the TAEs’ saturation phase in the Large Helical Device: MHD burst

Here, the aim of the present study is to analyze the saturation regime of the toroidal Alfven eigenmodes (TAE) in the Large Helical Device plasma, particularly the MHD burst. The linear and nonlinear evolution of the TAEs are simulated by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The linear simulations indicate the overlapping of 1/2–1/1, 2/3–2/4 and 3/5–3/6 TAEs in the inner-middle plasma region and frequency range of 45–75 kHz, triggered by EPs with an energy of T f = 45 keV and EP β = 0.022. The nonlinear simulations show that 2/3–2/4 and 3/4–3/5 TAEs are further destabilized due to the energy transfer from the 1/1–1/2 TAE, leading to broad TAE radial overlapping and triggering of the MHD burst. The energy of the 1/1–1/2 TAE is also nonlinearly transferred to the thermal plasma destabilizing the 0/0 and 0/1 modes, inducing the generation of shear flows and zonal currents, as well as large deformations in the thermal pressure and EP density radial profiles. The nonlinear simulation reproduces the same succession of instabilities and the same frequency range with respect to the experiment. The instability propagates outward during the bursting phase, showing a large decrease of the EP density profile between the middle-outer plasma, indicating the loss of part of the EP population that explains the decrease in the plasma heating efficiency observed during the MHD burst.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liner implosion experiments driven by a dynamic screw pinch

This paper expands upon recent experimental results [Campbell et al., Phys. Rev. Lett. 125, 035001 (2020)], where thin-foil liner implosions were driven by a dynamic screw pinch (DSP) and found to have magneto-Rayleigh–Taylor instability (MRTI) amplitudes up to three times smaller than in implosions driven by a standard z-pinch (SZP). Here, the expanded discussion presented herein includes: (1) a detailed comparison of the MRTI growth measured in the experiment with that calculated from theory; (2) measurements of axial magnetic field injection into the liner interior prior to the implosion, as well as the subsequent compression of this field during the implosion; (3) an in-depth description of how the helical geometry of the DSP can result in earlier implosion and stagnation times relative to the SZP; and (4) particle-in-cell simulations showing different electron drift behavior in the anode–cathode gap of the DSP relative to the SZP, and how this difference may be related to the different current waveforms recorded during the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Evaluation of Peer-to-Peer Distributed Microgrids Coordination for Voltage Regulation

This paper presents the performance evaluation of a peer-to-peer microgrids coordination algorithm for sub-transmission systems. As distributed energy resources (DERs) in distribution system start to show negative impact to the bulk power system, a paradigm shift is needed for transmission planning and operation. Because distribution substations are located far from the sub-transmission system, and it is hard to use traditional centralized control for real-time control and coordination. Thus, distributed control is a natural choice because it requires less communication and central computation. In this paper, each distribution substation is treated as a microgrid, and the peer-to-peer distributed microgrids control is formulated as a real-time optimal power flow problem to reduce the negative impact in sub-transmission systems. A distributed primal-dual optimization algorithm is adopted to solve the problem. Validation of the peer-to-peer algorithm is performed through the simulation of a real-world sub-transmission system composing of many distribution systems with high renewable penetration. Simulation results show that the peer-to-peer algorithm can achieve satisfactory performance (e.g., voltage regulation) in sub-transmission system by coordinating and controlling DERs in distribution systems.

distributed control↗

Performance Evaluation of Peer-to-Peer Distributed Microgrids Coordination for Voltage Regulation: Preprint

This paper presents the performance evaluation of a peer-to-peer microgrids coordination algorithm for sub-transmission systems. As distributed energy resources (DERs) in distribution system start to show negative impact to the bulk power system, a paradigm shift is needed for transmission planning and operation. Because distribution substations are located far from the sub-transmission system, and it is hard to use traditional centralized control for real-time control and coordination. Thus, distributed control is a natural choice because it requires less communication and central computation. In this paper, each distribution substation is treated as a microgrid, and the peer-to-peer distributed microgrids control is formulated as a real-time optimal power flow problem to reduce the negative impact in sub-transmission systems. A distributed primal-dual optimization algorithm is adopted to solve the problem. Validation of the peer-to-peer algorithm is performed through the simulation of a real-world sub-transmission system composing of many distribution systems with high renewable penetration. Simulation results show that the peer-to-peer algorithm can achieve satisfactory performance (e.g., voltage regulation) in sub-transmission system by coordinating and controlling DERs in distribution systems.

distributed control↗