Jumper is not working satisfactorily in ANDE

[PandaRP01]

Hi Everyone
I am working in ANDES for my dynamic simulation. I want to merge two buses through Jumper. Simulation is running fine. However, when I try to deactivate the Jumper using "ssPV.Jumper.alter("u", 1, 0)" . However, Jumper power is not zero after using this command.
I am using ANDES1.8.3 and also ANDES 1.7.6.

Could you suggest why Jumper is not working for me?

Is there any alternative to switch ON/OFF the line? As making u=1/0, cause simulation stop (converges to zero).

Thanks
PandaRP01

[cuihantao]

What's the error message?

[PandaRP01]

After running TDS for 10 second.. I ran below command and error massage is attached below
ssPV.Line.alter("u", "Line_5", 0)
ssPV.Line.alter("u", "Line_6", 1)
ssPV.TDS.config.tf =25
ssPV.TDS.run()

Error massage:

Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=10.3300 s.
Simulation to t=25.00 sec completed in 0.0660 seconds.

Also kindly suggest, the why Jumper is not deactivate after executing ssPV.Jumper.alter("u", 1, 0). It is not creating any error. However, there is no change in jumper power too.

[cuihantao]

Thanks. If you can test the develop branch and let me know. There was a recent fix on the time stepping algorithm.

[PandaRP01]

Hi Cui

Thank you so much for your kind response. I definitely report the issues while doing on develop branch.

Could you suggest, what is the issues with Jumper. Or, provide some more documents to on Jumper

[PandaRP01]

Hi Cui
Thank you for fixing the time stepping error while line switching.

I am using Jumper to merge two buses in ANDES simulation. When I ran below attached commands, it is not showing any error, however, the jumper power is not changed to zero also.

ssPV.Jumper.alter("u", 1, 0)
ssPV.TDS.config.tf =25
ssPV.TDS.run()

As I am changing the u of Jumper from 1 to zero. I am expecting that, Power through jumper must be zero.
Could you suggest, what are the issues with Jumper or did I make some mistakes while considering Jumper.

[cuihantao]

I pushed a change to my develop and should fix it. Please test and let me know.

[PandaRP01]

Hello Cui
Thank you for your kind response. I have installed ANDES 1.8.4 and repeat the method to disconnect the line. The commands are attached below

ssPV.Line.alter("u", "Line_5", 0)
ssPV.Line.alter("u", "Line_6", 1)
ssPV.TDS.config.tf =25
ssPV.TDS.run()

However, error massage is remain same. Could you kindly suggest the mistakes I am doing.
Error massage is attached below.

Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=10.0000 s.
Simulation to t=25.00 sec completed in 0.0500 seconds.

Also, after disconnection of this line, No Islanding in the system. system connectivity is Okay.
Kindly suggest the solution for this problem.

[cuihantao]

I cannot replicate the issue. Below is the code I used:

import andes

ss = andes.run(andes.get_case("ieee14/ieee14_jumper.xlsx"))

ss.TDS.run()

ss.Jumper.alter("u", 1, 0)
ss.TDS.config.tf =25
ss.TDS.run()



ss.Jumper.as_df()

ss.Jumper.algebs

ss.Jumper.u

ss.TDS.plt.plot(ss.Bus.v, a=(8, 9))

ss.TDS.plt.plot(ss.Jumper.p)

ss.Line.alter("u", "Line_5", 0)
ss.Line.alter("u", "Line_6", 1)
ss.TDS.config.tf =40
ss.TDS.run()


ss.TDS.plt.plot(ss.Bus.a)

Please post the full minimal working example that I can run.

[PandaRP01]

Hi Cui

Thank you so much for your effort.
After discussing with my colleagues, I get back to you with minimal working examples to replicate our issues.

Regards
Rakesh Panda

[PandaRP01]

Hello Mr. Cui

Thank you for resolving the Jumper issue in the ANDES 1.8.5. release. Now, Jumper is working fine. However I am facing same time step error as previously mentioned. For demonstration of my time step issues, I took kundur_full.xlsx system
I am trying to simulate kundur_full.xlsx file with toggler. I am using ANDES 1.8.5 with ss.TDS.config.tstep = 0.33, ss.TDS.config.shrinkt = 0, ss.TDS.config.check_conn= 1, ss.TDS.config.tf =5. It is showing same time step error. The error is attached below.

Time Domain Simulation Summary:
Sparse Solver: KLU
Simulation time: 0.0-5 s.
Fixed step size: h=330 ms.
Initialization for dynamics completed in 0.1330 seconds.
Initialization was successful.

Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=2.0001 s.
Simulation to t=5.00 sec completed in 0.0980 seconds.

Also when I made toggler off and run the system and made line_8, u=0 at 2 sec using ss.Line.alter("u", "Line_8", 0). The same error has been occurred. Could you please suggest why this error is occurring in my system?

[cuihantao]

Since you turned off adaptive time stepping (shrinkt=0), the program will terminate if a non-convergence is detected.

What you can do is to turn on shrinkt or set tstep to a smaller value like 0.033.

[PandaRP01]

Hello Mr. Cui

Thank you for highlighting my mistakes. However, I turned ON the adaptive time stepping (shrinkt=1). However, the error is still persist
Case-1: ss.TDS.config.shrinkt = 1, ss.TDS.config.tstep = 0.33. The errors are

Time Domain Simulation Summary:
Sparse Solver: KLU
Simulation time: 0.0-5 s.
Fixed step size: h=330 ms. Shrink if not converged.
Initialization for dynamics completed in 0.0980 seconds.
Initialization was successful. 

Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=2.0111 s.
Simulation to t=5.00 sec completed in 3.7604 seconds.

Case-2 : ss.TDS.config.shrinkt = 1, ss.TDS.config.tstep = 0.033. The errors are

-> Time Domain Simulation Summary:
Sparse Solver: KLU
Simulation time: 0.0-5 s.
Fixed step size: h=33 ms. Shrink if not converged.
Initialization for dynamics completed in 0.0990 seconds.
Initialization was successful.
Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=2.0111 s.
Simulation to t=5.00 sec completed in 3.0229 seconds.

Case-3: ss.TDS.config.shrinkt = 0, ss.TDS.config.tstep = 0.033. The errors are

-> Time Domain Simulation Summary:
Sparse Solver: KLU
Simulation time: 0.0-5 s.
Fixed step size: h=33 ms.
Initialization for dynamics completed in 0.1130 seconds.
Initialization was successful.
Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=2.0001 s.
Simulation to t=5.00 sec completed in 0.2601 seconds.

Kindly suggest the solution for the issues

Regards
Rakesh Panda

[cuihantao]

You may run the second case with andes -v 10 run ... and check the debugging output. You might paste the part for time-domain simulation here as well. It may be possible that the system does not have a power flow solution after the line disconnection.

[PandaRP01]

Hello Mr. Cui

I also run the simulation in command window with command
andes -v 10 run D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx

the power flow is also converging and the massage is

(ANDES185) C:\Users\rakeshpa>andes -v 10 run D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx
Namespace(add_book=None, addfile=None, command='run', config_option='', convert='', convert_all='', dime_address=None, filename=['D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx'], flat=False, from_csv=None, init=False, input_path='', kqrt=None, ncpu=2, no_output=False, no_pbar=False, no_preamble=False, output_path='', pert='', pool=False, profile=False, qrt=False, routine=('pflow',), shell=False, state_matrix=False, tf=None, verbose=10)

    _           _         | Version 1.8.5
   /_\  _ _  __| |___ ___ | Python 3.8.15 on Windows, 01/11/2023 12:55:03 AM
  / _ \| ' \/ _` / -_|_-< |
 /_/ \_\_||_\__,_\___/__/ | This program comes with ABSOLUTELY NO WARRANTY.

Logging to file "C:\Users\rakeshpa\AppData\Local\Temp\2\andes-gknw5wv2\andes.log"
Working directory: "C:\Users\rakeshpa"
Found files: ['D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx']
> Loaded config from file "andes.rc"
> Loaded generated Python code in "C:\Users\rakeshpa\.andes\pycode".
Generated code for <Fault> is stale.
Numerical code generation (rapid incremental mode) started...
Generating code for 1 models on 2 processes.
Saved generated pycode to "C:\Users\rakeshpa\.andes\pycode"
> Reloaded generated Python code of module "pycode".
Collecting Model.calls into System.
Generated numerical code for 1 models in 2.9535 seconds.
Input format guessed as xlsx.
Parsing input file "D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx"...
Input file parsed in 1.9690 seconds.
Setting internal address for Bus
Setting internal address for PQ
Setting internal address for PV
Setting internal address for Slack
Setting internal address for Line
Setting internal address for Area
System internal structure set up in 0.1093 seconds.
Entering connectivity check.
-> System connectivity check results:
  No islanded bus detected.
  System is interconnected.
  Bus indices in interconnected system (0-based): [[0, 1, 2, 3, 4, 5, 6, 7, 8, 9]]
  Each island has a slack bus correctly defined and enabled.

-> Power flow calculation
           Numba: Off
   Sparse solver: KLU
 Solution method: NR method
========== Bus has <pflow_init> = True ==========
Initialization sequence:
a -> v
Bus: assignment initialization phase begins
========== PQ has <pflow_init> = True ==========
Initialization sequence:
a -> v
PQ: assignment initialization phase begins
========== PV has <pflow_init> = True ==========
Initialization sequence:
q -> a -> v
PV: assignment initialization phase begins
========== Slack has <pflow_init> = True ==========
Initialization sequence:
q -> p -> a -> v
Slack: assignment initialization phase begins
========== Line has <pflow_init> = True ==========
Initialization sequence:
a1 -> a2 -> v1 -> v2
Line: assignment initialization phase begins
========== Area has <pflow_init> = True ==========
Initialization sequence:

Area: assignment initialization phase begins
Power flow initialized in 0.0312 seconds.
Jacobian updated at t=0.000000.
Max. algeb mismatch 14.9282832 on v Bus 4
0: |F(x)| = 14.9282832
Jacobian updated at t=0.000000.
Max. algeb mismatch 3.608627841 on v Bus 10
1: |F(x)| = 3.608627841
Jacobian updated at t=0.000000.
Max. algeb mismatch 0.1701107882 on v Bus 7
2: |F(x)| = 0.1701107882
Jacobian updated at t=0.000000.
Max. algeb mismatch 0.002038626956 on a Bus 8
3: |F(x)| = 0.002038626956
Jacobian updated at t=0.000000.
Max. algeb mismatch 3.745104006e-07 on a Bus 8
4: |F(x)| = 3.745104006e-07
Converged in 5 iterations in 0.0312 seconds.
Report saved to "kundur_full_out.txt" in 0.0000 seconds.
-> Single process finished in 6.2508 seconds.

The time domain code is attached below

import andes
import matplotlib
matplotlib.use('nbAgg')

andes.__version__

from andes.utils.paths import get_case

andes.config_logger(stream_level=20)


ss = andes.run(andes.get_case("kundur/kundur_full.xlsx"), setup=False, default_config=True)
#ss = andes.load('D:/Simulator_project/SimuFiles/Phase2/kundur_full.xlsx',
#               setup=False, config_option=["ESD1.plim=0", "System.mva=100"],
#                no_output=False)
    
ss.setup()

# use constant power model for PQ
ss.PQ.config.p2p = 1
ss.PQ.config.q2q = 1
ss.PQ.config.p2z = 0
ss.PQ.config.q2z = 0

# turn off under-voltage PQ-to-Z conversion
ss.PQ.pq2z = 0
ss.PFlow.run()

ss.PVD1.Iqcmd.lim.allow_adjust

# First running the simulation for 5s
ss.TDS.config.tstep = 0.033
ss.TDS.config.shrinkt = 1
ss.TDS.config.check_conn= 1
#ssPV.TDS.config.honest=1
#ssPV.TDS.config.qrt = 1
#ssPV.TDS.config.max_iter = 100
ss.TDS.config.tf =5
ss.TDS.run()

The massage is

Working directory: "C:\Users\rakeshpa\simu_pr\Phase_2"
> Loaded generated Python code in "C:\Users\rakeshpa\.andes\pycode".
Parsing input file "C:\Users\rakeshpa\.conda\envs\ANDES185\lib\site-packages\andes\cases\kundur\kundur_full.xlsx"...
Input file parsed in 1.4274 seconds.
System is not set up. Routines cannot continue.
System internal structure set up in 0.1090 seconds.
-> System connectivity check results:
  No islanded bus detected.
  System is interconnected.
  Each island has a slack bus correctly defined and enabled.

-> Power flow calculation
           Numba: Off
   Sparse solver: KLU
 Solution method: NR method
Power flow initialized in 0.0180 seconds.
0: |F(x)| = 14.9282832
1: |F(x)| = 3.608627841
2: |F(x)| = 0.1701107882
3: |F(x)| = 0.002038626956
4: |F(x)| = 3.745104006e-07
Converged in 5 iterations in 0.0250 seconds.
Report saved to "kundur_full_out.txt" in 0.0080 seconds.
-> Single process finished in 2.1598 seconds.

-> Time Domain Simulation Summary:
Sparse Solver: KLU
Simulation time: 0.0-5 s.
Fixed step size: h=33 ms. Shrink if not converged.
Initialization for dynamics completed in 0.1141 seconds.
Initialization was successful.
41%
40.66/100 [00:02<00:00, 197.24%/s]
<Toggle 1>: Line.Line_8 status changed to 0 at t=2.0 sec.
Time step reduced to zero. Convergence is not likely.
Simulation terminated at t=2.0111 s.
Simulation to t=5.00 sec completed in 2.9131 seconds.
Outputs to "kundur_full_out.lst" and "kundur_full_out.npz".
Outputs written in 0.0040 seconds.

[cuihantao]

Hello, in the second part, you did not have verbose on.

If you'd like to share a minimum working example with me, you can send it to cuihantao@gmail.com or post it here. I will help you debug it.