Affiliations
AffiliationItem in Clipboard
Suppressing cascades of load in interdependent networksCharles D Brummitt et al. Proc Natl Acad Sci U S A. 2012.
. 2012 Mar 20;109(12):E680-9. doi: 10.1073/pnas.1110586109. Epub 2012 Feb 21. AffiliationItem in Clipboard
AbstractUnderstanding how interdependence among systems affects cascading behaviors is increasingly important across many fields of science and engineering. Inspired by cascades of load shedding in coupled electric grids and other infrastructure, we study the Bak-Tang-Wiesenfeld sandpile model on modular random graphs and on graphs based on actual, interdependent power grids. Starting from two isolated networks, adding some connectivity between them is beneficial, for it suppresses the largest cascades in each system. Too much interconnectivity, however, becomes detrimental for two reasons. First, interconnections open pathways for neighboring networks to inflict large cascades. Second, as in real infrastructure, new interconnections increase capacity and total possible load, which fuels even larger cascades. Using a multitype branching process and simulations we show these effects and estimate the optimal level of interconnectivity that balances their trade-offs. Such equilibria could allow, for example, power grid owners to minimize the largest cascades in their grid. We also show that asymmetric capacity among interdependent networks affects the optimal connectivity that each prefers and may lead to an arms race for greater capacity. Our multitype branching process framework provides building blocks for better prediction of cascading processes on modular random graphs and on multitype networks in general.
Conflict of interest statementThe authors declare no conflict of interest.
FiguresFig. 1.
The power grid of the…
Fig. 1.
The power grid of the continental United States, illustrating the three main regions…
Fig. 1.The power grid of the continental United States, illustrating the three main regions or “interconnects”—Western, Eastern, and Texas—and new lines (in red) proposed by American Electric Power to transport wind power (11, 13).
Fig. 2.
A random three- and four-regular…
Fig. 2.
A random three- and four-regular graph connected by Bernoulli-distributed coupling with interconnectivity parameter…
Fig. 2.A random three- and four-regular graph connected by Bernoulli-distributed coupling with interconnectivity parameter p = 0.1 [R(3)-B(0.1)-R(4)]. We also illustrate the shedding branch distribution qab(rba,bb), the chance that an ab shedding causes rba, rbb many ba, bb sheddings at the next time step. Note these random graphs are small and become tree-like when they are large (⪆1,000 nodes).
Fig. 3.
The multitype branching process (red…
Fig. 3.
The multitype branching process (red curves) approximates simulations of sandpile cascades on the…
Fig. 3.The multitype branching process (red curves) approximates simulations of sandpile cascades on the power grids
c,
d(blue curves) surprisingly well, given that power grids are among the most difficult network topologies on which to predict other dynamics (54). We plot the four marginalized avalanche size distributions (using the labels
a,
bfor grids
c,
d). For example, Fig. 3
A, plotting
versus
ta, is the chance an avalanche begun in
atopples
tanodes in
a. Vertical axes’ labels are the expressions in the plots. Here dissipation
f= 0.1.
Fig. 4.
Interconnectivity is locally stabilizing, but…
Fig. 4.
Interconnectivity is locally stabilizing, but only up to a critical point. The main…
Fig. 4.Interconnectivity is locally stabilizing, but only up to a critical point. The main plot, the results of simulations on
R(3)-
B(
p)-
R(3) (2 × 10
6grains;
f= 0.01; 2 × 10
3nodes/network), shows that large local cascades decrease and then increase with
p, whereas large inflicted cascades only become more likely. Their average (gold curve) has a stable minimum at
p∗≈ 0.075 ± 0.01. This curve and its critical point,
p∗, are stable for cutoffs
in the regime
. (
Insets) Rank-size plots on log–log scales of the largest cascades in network
a(
Left) for
p= 10
-3, 10
-2, 10
-1and in power grid
dconnected to
cby 0, 8, or 16 edges.
Fig. 5.
Increasing the interconnectivity p between…
Fig. 5.
Increasing the interconnectivity p between two random three-regular graphs extends the tail of…
Fig. 5.Increasing the interconnectivity p between two random three-regular graphs extends the tail of the total avalanche size distribution s(t), which does not distinguish whether toppled nodes are in network a or b. The inset shows a rank-size plot on log–log scales of the number of topplings t in the largest 104 avalanches (with 2 × 106 grains of sand dropped), showing that adding more edges between random three-regular graphs enlarges the largest global cascades by an amount on the order of the additional number of interconnections. As expected theoretically (39), when a and b nodes are viewed as one network, s(t) ∼ t-3/2 for large t (green line).
Fig. 6.
( A ) Networks mitigating…
Fig. 6.
( A ) Networks mitigating the smallest cascades 1 ≤ T a ≤…
Fig. 6.(A) Networks mitigating the smallest cascades 1 ≤ Ta ≤ 51 seek isolation p = 0, whereas (B) networks suppressing intermediate cascades 100 ≤ Ta ≤ 150 seek isolation p = 0 or strong coupling p = 1, depending on the initial interconnectivity p in relation to the unstable critical point p∗ ≈ 0.05. But networks like power grids mitigating large cascades (C and D) would seek interconnectivity at the stable equilibrium p∗ ≈ 0.075 ± 0.01. The bottom figures and the location of p∗ are robust to changes in the window ℓ ≤ Ta ≤ ℓ + 50 for all 400 ≤ ℓ ≤ 1,500.
Fig. 7.
The critical load disparity at…
Fig. 7.
The critical load disparity at which inflicted cascades in d become equally large…
Fig. 7.The critical load disparity at which inflicted cascades in d become equally large as its local cascades is r∗ ≈ 15 (for 16 interconnections, “bridges”). (For eight interconnections, r∗ > 20; see Fig. S9 ). Here we show a rank-size plot in log–log scales of the largest 104 avalanches in power grid d, distinguishing whether they begin in c (inflicted cascades) or in d (local cascades), for 8 and 16 interconnections (solid, dashed curves), for r = 1 (Main), r = 15 (Inset), in simulations with dissipation of sand f = 0.05, 106 grains dropped after 105 transients.
Fig. P1.
The chance that a network…
Fig. P1.
The chance that a network a coupled to another network b suffers a…
Fig. P1.The chance that a network a coupled to another network b suffers a cascade spanning more than half the network (gold curve) has a stable minimum at a critical amount of interconnectivity p∗. Networks seeking to mitigate their largest cascades would prefer to build or demolish interconnections to operate at this critical point p∗. The blue (red) curve is the chance that a cascade that begins in a (b) topples at least 1,000 nodes in a. Increasing interconnectivity only exacerbates the cascades inflicted from b to a (red), but interestingly, it initially suppresses the local cascades in a. (From simulations on coupled random three-regular graphs; the inset depicts a small example with 30 nodes per network and p = 0.1).
Similar articlesTuralska M, Swami A. Turalska M, et al. Sci Rep. 2021 Feb 8;11(1):3276. doi: 10.1038/s41598-021-82843-8. Sci Rep. 2021. PMID: 33558578 Free PMC article.
Turalska M, Burghardt K, Rohden M, Swami A, D'Souza RM. Turalska M, et al. Phys Rev E. 2019 Mar;99(3-1):032308. doi: 10.1103/PhysRevE.99.032308. Phys Rev E. 2019. PMID: 30999482
Mikaberidze G, D'Souza RM. Mikaberidze G, et al. Chaos. 2022 May;32(5):053121. doi: 10.1063/5.0095094. Chaos. 2022. PMID: 35649989
Tan F, Xia Y, Wei Z. Tan F, et al. Phys Rev E Stat Nonlin Soft Matter Phys. 2015 May;91(5):052809. doi: 10.1103/PhysRevE.91.052809. Epub 2015 May 18. Phys Rev E Stat Nonlin Soft Matter Phys. 2015. PMID: 26066214
Peng H, Liu C, Zhao D, Hu Z, Han J. Peng H, et al. Sensors (Basel). 2020 Oct 28;20(21):6123. doi: 10.3390/s20216123. Sensors (Basel). 2020. PMID: 33126431 Free PMC article.
Boccaletti S, Bianconi G, Criado R, Del Genio CI, Gómez-Gardeñes J, Romance M, Sendiña-Nadal I, Wang Z, Zanin M. Boccaletti S, et al. Phys Rep. 2014 Nov 1;544(1):1-122. doi: 10.1016/j.physrep.2014.07.001. Epub 2014 Jul 10. Phys Rep. 2014. PMID: 32834429 Free PMC article. Review.
Karsai M, Iñiguez G, Kikas R, Kaski K, Kertész J. Karsai M, et al. Sci Rep. 2016 Jun 7;6:27178. doi: 10.1038/srep27178. Sci Rep. 2016. PMID: 27272744 Free PMC article.
Min B, Zheng M. Min B, et al. PLoS One. 2018 Apr 18;13(4):e0195539. doi: 10.1371/journal.pone.0195539. eCollection 2018. PLoS One. 2018. PMID: 29668730 Free PMC article.
Parra-Morales AM, Rudas J, Vargas JA, Gómez F, Enciso-Olivera CO, Trujillo-Rodriguez D, Martínez D, Hernandez J, Ordóñez-Rubiano EG, Marín-Muñoz JH. Parra-Morales AM, et al. Medicine (Baltimore). 2019 May;98(19):e15620. doi: 10.1097/MD.0000000000015620. Medicine (Baltimore). 2019. PMID: 31083258 Free PMC article.
Cho JH, Gao J. Cho JH, et al. PLoS One. 2016 Feb 9;11(2):e0148674. doi: 10.1371/journal.pone.0148674. eCollection 2016. PLoS One. 2016. PMID: 26859840 Free PMC article.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4