def setup_simulation(self, scenario='default'):
"""Setup a simulation."""
self.events_list = []
self.grid_list = []
self.DER_model_list = []
self.sim_list = []
self.results_list = []
self.n_instances = self.return_settings(scenario=scenario,
parameter='n_DER',
settings_type='setup')
power_rating = self.return_settings(scenario=scenario,
parameter='power_rating',
settings_type='setup')
SinglePhase = self.return_settings(scenario=scenario,
parameter='SinglePhase',
settings_type='setup')
SteadyState = self.return_settings(scenario=scenario,
parameter='SteadyState',
settings_type='setup')
for i in range(self.n_instances):
self.events_list.append(SimulationEvents())
self.grid_list.append(Grid(events=self.events_list[-1]))
if SinglePhase:
self.DER_model_list.append(
SolarPV_DER_SinglePhase(
grid_model=self.grid_list[-1],
events=self.events_list[-1],
standAlone=True,
Sinverter_rated=power_rating,
STEADY_STATE_INITIALIZATION=SteadyState))
else:
self.DER_model_list.append(
SolarPV_DER_ThreePhase(
grid_model=self.grid_list[-1],
events=self.events_list[-1],
standAlone=True,
Sinverter_rated=power_rating,
STEADY_STATE_INITIALIZATION=SteadyState))
self.sim_list.append(
DynamicSimulation(grid_model=self.grid_list[-1],
PV_model=self.DER_model_list[-1],
events=self.events_list[-1],
LOOP_MODE=False,
COLLECT_SOLUTION=True))
self.sim_list[
-1].jacFlag = False #Provide analytical Jacobian to ODE solver
self.sim_list[
-1].DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results_list.append(
SimulationResults(simulation=self.sim_list[-1], PER_UNIT=True))
def test_run_simulation(self):
"""Test run simulation method."""
events = SimulationEvents()
PVDER = SolarPVDERThreePhase(events=events,
configFile=config_file,
**{
**self.flag_arguments,
**self.ratings_arguments,
**self.voltage_arguments
})
sim = DynamicSimulation(PV_model=PVDER,
events=events,
jacFlag=True,
verbosity='DEBUG',
solverType='odeint')
sim.tStop = 10.0
sim.tInc = 1 / 120.
sim.run_simulation()
self.assertEqual(sim.t[-1], sim.tStop)
self.assertTrue(sim.SOLVER_CONVERGENCE)
def test_init(self):
"""Test Dynamic Simulation initialization."""
events = SimulationEvents()
PVDER = SolarPV_DER_ThreePhase(events = events,
Sinverter_rated = self.power_rating,Vrms_rated = self.Vrms, #175
gridVoltagePhaseA = self.Va,
gridVoltagePhaseB = self.Vb,
gridVoltagePhaseC = self.Vc,
gridFrequency = self.wgrid,
standAlone = False,STEADY_STATE_INITIALIZATION=True,verbosity = 'DEBUG')
sim = DynamicSimulation(PV_model=PVDER,events = events,jacFlag = True,verbosity = 'DEBUG',solver_type='odeint')
self.assertIsInstance(sim, DynamicSimulation)
self.assertTrue(sim.jacFlag)
def create_DER_simulation_instances(self, n_instances):
"""Create specified number of DER instances."""
assert n_instances >= 1, 'Number of simulation instances should be greater than or equal to one!'
self.events_list = []
self.grid_list = []
self.DER_list = []
self.sim_list = []
print('{} DER model instances will be created.'.format(n_instances))
for instance in range(n_instances):
print('Creating Instance:', instance + 1)
self.events_list.append(SimulationEvents())
self.grid_list.append(
Grid(events=self.events_list[instance],
unbalance_ratio_b=1.0,
unbalance_ratio_c=1.0))
if self._SINGLE_PHASE:
self.DER_list.append(
SolarPV_DER_SinglePhase(
grid_model=self.grid_list[instance],
events=self.events_list[instance],
Sinverter_rated=self._rating,
standAlone=True,
STEADY_STATE_INITIALIZATION=True,
verbosity='INFO'))
else:
self.DER_list.append(
SolarPV_DER_ThreePhase(grid_model=self.grid_list[instance],
events=self.events_list[instance],
Sinverter_rated=self._rating,
standAlone=True,
STEADY_STATE_INITIALIZATION=True,
verbosity='INFO'))
self.sim_list.append(
DynamicSimulation(grid_model=self.grid_list[instance],
PV_model=self.DER_list[instance],
events=self.events_list[instance]))
def test_init(self):
"""Test Dynamic Simulation initialization."""
events = SimulationEvents()
PVDER = SolarPVDERThreePhase(events=events,
configFile=config_file,
**{
**self.flag_arguments,
**self.ratings_arguments,
**self.voltage_arguments
})
sim = DynamicSimulation(PV_model=PVDER,
events=events,
jacFlag=True,
verbosity='DEBUG',
solverType='odeint')
self.assertIsInstance(sim, DynamicSimulation)
self.assertTrue(sim.jacFlag)
def setup_PVDER_simulation(self, model_type='model_2'):
"""Setup simulation environment."""
self.max_sim_time = self.max_sim_time_user
events = SimulationEvents(events_spec=self.env_events_spec,
verbosity='INFO')
grid_model = Grid(events=events)
PVDER_model = DERModel(
modelType=self.env_model_spec[model_type]['DERModelType'],
events=events,
configFile=self.env_model_spec[model_type]['configFile'],
gridModel=grid_model,
derId=self.env_model_spec[model_type]['derId'],
standAlone=True,
steadyStateInitialization=True)
PVDER_model.DER_model.LVRT_ENABLE = False #Disconnects PV-DER using ride through settings during voltage anomaly
PVDER_model.DER_model.DO_EXTRA_CALCULATIONS = True
#PV_model.Vdc_EXTERNAL = True
self.sim = DynamicSimulation(gridModel=grid_model,
PV_model=PVDER_model.DER_model,
events=events,
verbosity='INFO',
solverType='odeint',
LOOP_MODE=False) #'odeint','ode-vode-bdf'
self.sim.jacFlag = True #Provide analytical Jacobian to ODE solver
self.sim.DEBUG_SOLVER = False #Give information on solver convergence
self.results = SimulationResults(simulation=self.sim)
self.results.PER_UNIT = False
self.results.font_size = 18
self.generate_simulation_events()
PVDER_model.Qref_EXTERNAL = True #Enable VAR reference manipulation through outside program
self.sim.DEBUG_SOLVER = False #Give information on solver convergence
self.sim.tInc = self.env_sim_spec[
'sim_time_step'] #1/60.0 #self.time_taken_per_step #
def test_run_simulation(self):
"""Test run simulation method."""
events = SimulationEvents()
PVDER = SolarPV_DER_ThreePhase(events = events,
Sinverter_rated = self.power_rating,Vrms_rated = self.Vrms, #175
gridVoltagePhaseA = self.Va,
gridVoltagePhaseB = self.Vb,
gridVoltagePhaseC = self.Vc,
gridFrequency = self.wgrid,
standAlone = False,STEADY_STATE_INITIALIZATION=True,verbosity = 'DEBUG')
sim = DynamicSimulation(PV_model=PVDER,events = events,
jacFlag = True,verbosity = 'DEBUG',solver_type='odeint')
sim.tStop = 10.0
sim.tInc = 1/120.
sim.run_simulation()
self.assertEqual(sim.t[-1],sim.tStop)
self.assertTrue(sim.SOLVER_CONVERGENCE)
def setup(self, nodeid):
try:
VpuInitial = 1.0
SinglePhase = False
if 'myconfig' in OpenDSSData.config and 'DERParameters' in OpenDSSData.config[
'myconfig']:
if 'DERParameters' in OpenDSSData.config['myconfig']:
pvderConfig = OpenDSSData.config['myconfig'][
'DERParameters']
power_rating = OpenDSSData.config['myconfig'][
'DERParameters']['power_rating'] * 1e3
voltage_rating = OpenDSSData.config['myconfig'][
'DERParameters']['voltage_rating']
SteadyState = OpenDSSData.config['myconfig'][
'DERParameters']['SteadyState']
if 'nodenumber' in OpenDSSData.config['myconfig']:
DER_location = str(os.getpid()) + '-' + 'bus_' + str(
OpenDSSData.config['myconfig']
['nodenumber']) + '-' + 'node_' + nodeid
else:
DER_location = str(os.getpid()) + '-' + 'node_' + nodeid
else:
print(
'DERParameters not found in `OpenDSSData` object - using default ratings and parameters!'
)
pvderConfig = None
SteadyState = True
DER_location = 'node_' + nodeid
if SinglePhase:
power_rating = 10.0e3
else:
power_rating = 50.0e3
Va = (.50 + 0j) * Grid.Vbase
Vb = (-.25 - .43301270j) * Grid.Vbase
Vc = (-.25 + .43301270j) * Grid.Vbase
events = SimulationEvents()
if SinglePhase:
self.PV_model = PV_model = SolarPV_DER_SinglePhase(
events=events,
Sinverter_rated=power_rating,
Vrms_rated=voltage_rating,
gridVoltagePhaseA=Va * VpuInitial,
gridVoltagePhaseB=Vb * VpuInitial,
gridVoltagePhaseC=Vc * VpuInitial,
gridFrequency=2 * math.pi * 60.0,
standAlone=False,
STEADY_STATE_INITIALIZATION=SteadyState,
pvderConfig=pvderConfig,
identifier=DER_location)
else:
self.PV_model = PV_model = SolarPV_DER_ThreePhase(
events=events,
Sinverter_rated=power_rating,
Vrms_rated=voltage_rating,
gridVoltagePhaseA=Va * VpuInitial,
gridVoltagePhaseB=Vb * VpuInitial,
gridVoltagePhaseC=Vc * VpuInitial,
gridFrequency=2 * math.pi * 60.0,
standAlone=False,
STEADY_STATE_INITIALIZATION=SteadyState,
pvderConfig=pvderConfig,
identifier=DER_location)
self.PV_model.LVRT_ENABLE = True #Disconnects PV-DER based on ride through settings in case of voltage anomaly
self.sim = DynamicSimulation(PV_model=PV_model,
events=events,
LOOP_MODE=True,
COLLECT_SOLUTION=True)
self.sim.jacFlag = True #Provide analytical Jacobian to ODE solver
self.sim.DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results = SimulationResults(simulation=self.sim,
PER_UNIT=True)
self.lastSol = copy.deepcopy(self.sim.y0) # mutable,make copy
self.lastT = 0
except Exception as e:
OpenDSSData.log("Failed Setup PVDER at node:{}!".format(nodeid))
class PVDERModel:
def __init__(self):
self._pvModel = {}
self._sim = {}
self._lastSol = {}
self._lastT = 0
def setup(self, nodeid):
try:
VpuInitial = 1.0
SinglePhase = False
if 'myconfig' in OpenDSSData.config and 'DERParameters' in OpenDSSData.config[
'myconfig']:
if 'DERParameters' in OpenDSSData.config['myconfig']:
pvderConfig = OpenDSSData.config['myconfig'][
'DERParameters']
power_rating = OpenDSSData.config['myconfig'][
'DERParameters']['power_rating'] * 1e3
voltage_rating = OpenDSSData.config['myconfig'][
'DERParameters']['voltage_rating']
SteadyState = OpenDSSData.config['myconfig'][
'DERParameters']['SteadyState']
if 'nodenumber' in OpenDSSData.config['myconfig']:
DER_location = str(os.getpid()) + '-' + 'bus_' + str(
OpenDSSData.config['myconfig']
['nodenumber']) + '-' + 'node_' + nodeid
else:
DER_location = str(os.getpid()) + '-' + 'node_' + nodeid
else:
print(
'DERParameters not found in `OpenDSSData` object - using default ratings and parameters!'
)
pvderConfig = None
SteadyState = True
DER_location = 'node_' + nodeid
if SinglePhase:
power_rating = 10.0e3
else:
power_rating = 50.0e3
Va = (.50 + 0j) * Grid.Vbase
Vb = (-.25 - .43301270j) * Grid.Vbase
Vc = (-.25 + .43301270j) * Grid.Vbase
events = SimulationEvents()
if SinglePhase:
self.PV_model = PV_model = SolarPV_DER_SinglePhase(
events=events,
Sinverter_rated=power_rating,
Vrms_rated=voltage_rating,
gridVoltagePhaseA=Va * VpuInitial,
gridVoltagePhaseB=Vb * VpuInitial,
gridVoltagePhaseC=Vc * VpuInitial,
gridFrequency=2 * math.pi * 60.0,
standAlone=False,
STEADY_STATE_INITIALIZATION=SteadyState,
pvderConfig=pvderConfig,
identifier=DER_location)
else:
self.PV_model = PV_model = SolarPV_DER_ThreePhase(
events=events,
Sinverter_rated=power_rating,
Vrms_rated=voltage_rating,
gridVoltagePhaseA=Va * VpuInitial,
gridVoltagePhaseB=Vb * VpuInitial,
gridVoltagePhaseC=Vc * VpuInitial,
gridFrequency=2 * math.pi * 60.0,
standAlone=False,
STEADY_STATE_INITIALIZATION=SteadyState,
pvderConfig=pvderConfig,
identifier=DER_location)
self.PV_model.LVRT_ENABLE = True #Disconnects PV-DER based on ride through settings in case of voltage anomaly
self.sim = DynamicSimulation(PV_model=PV_model,
events=events,
LOOP_MODE=True,
COLLECT_SOLUTION=True)
self.sim.jacFlag = True #Provide analytical Jacobian to ODE solver
self.sim.DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results = SimulationResults(simulation=self.sim,
PER_UNIT=True)
self.lastSol = copy.deepcopy(self.sim.y0) # mutable,make copy
self.lastT = 0
except Exception as e:
OpenDSSData.log("Failed Setup PVDER at node:{}!".format(nodeid))
def prerun(self, gridVoltagePhaseA, gridVoltagePhaseB, gridVoltagePhaseC):
"""Prerun will set the required data in pvder model. This method should be run before
running the integrator."""
try:
# set grid voltage. This value will be kept constant during the run
# as opendss will sync will grid_simulation only once during the run call.
# opendss will solve power flow and will set gridVoltagePhase value.
self.PV_model.gridVoltagePhaseA = gridVoltagePhaseA * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseB = gridVoltagePhaseB * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseC = gridVoltagePhaseC * math.sqrt(
2) / Grid.Vbase
self.sim.t = self.lastT + 1 / 120.0
return None
except Exception as e:
OpenDSSData.log("Failed prerun in the pvder model")
def postrun(self, sol, t, KVAbase=50):
"""Postrun will gather results. This method should be run after running the integrator."""
try:
if self.sim.COLLECT_SOLUTION:
self.sim.collect_solution(sol, t)
#np array is mutable, hence make a copy
self.lastSol = copy.deepcopy(sol[-1])
self.lastT = t[-1]
# get S
S = self._getS()
return S * KVAbase
except Exception as e:
OpenDSSData.log("Failed postrun in the pvder model")
def _run(self,
gridVoltagePhaseA,
gridVoltagePhaseB,
gridVoltagePhaseC,
KVAbase=50):
#t = [self.lastT, self.lastT + OpenDSSData.data['DNet']['DER']['DERParameters']['dt']]
t = [self.lastT, self.lastT + 1 / 120.0]
try:
# set grid voltage. This value will be kept constant during the run
# as opendss will sync will grid_simulation only once during the run call.
# opendss will solve power flow and will set gridVoltagePhase value.
self.PV_model.gridVoltagePhaseA = gridVoltagePhaseA * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseB = gridVoltagePhaseB * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseC = gridVoltagePhaseC * math.sqrt(
2) / Grid.Vbase
self.sim.t = t
sol, info, ConvergeFlag = self.sim.call_ODE_solver(
self.sim.ODE_model, self.sim.jac_ODE_model, self.lastSol, t)
if self.sim.COLLECT_SOLUTION:
self.sim.collect_solution(sol, t)
#sol,info=odeint(odeFunc,y0,t,Dfun=jac_odeFunc,full_output=1,printmessg=True,hmax = 1/120.,mxstep=50,atol=1e-4,rtol=1e-4)
#sol,info=odeint(odeFunc,y0,t,full_output=1,printmessg=True,hmax = 1/120.,mxstep=50,atol=1e-4,rtol=1e-4)
#np array is mutable, hence make a copy
self.lastSol = copy.deepcopy(sol[-1])
self.lastT = t[-1]
# get S
S = self._getS()
return S * KVAbase
except Exception as e:
#OpenDSSData.log("Failed run in the pvder model")
OpenDSSData.log(
"{}:ODE solver failed between {:.4f} s and {:.4f} s!".format(
self.PV_model.name, t[0], t[-1]))
S = 0
return S
def _getS(self):
try:
return self.sim.PV_model.S_PCC # value at end of PV-DER simulation
except Exception as e:
OpenDSSData.log("Failed getS in the pvder model")
def setup_simulation(self, scenario='default'):
"""Setup a simulation."""
self.events_list = []
self.grid_list = []
self.DER_model_list = []
self.sim_list = []
self.results_list = []
self.n_instances = self.return_settings(scenario=scenario,
parameter='n_DER',
settings_type='setup')
flag_arguments = {
'standAlone':
True,
'steadyStateInitialization':
self.return_settings(scenario=scenario,
parameter='steadyStateInitialization',
settings_type='setup'),
'verbosity':
'DEBUG'
}
ratings_arguments = {
'powerRating':
self.return_settings(scenario=scenario,
parameter='powerRating',
settings_type='setup')
}
SinglePhase = self.return_settings(scenario=scenario,
parameter='SinglePhase',
settings_type='setup')
for i in range(self.n_instances):
self.events_list.append(SimulationEvents())
self.grid_list.append(Grid(events=self.events_list[-1]))
if SinglePhase:
self.DER_model_list.append(
SolarPVDERSinglePhase(events=self.events_list[-1],
configFile=config_file,
**{
"gridModel": self.grid_list[-1],
"identifier": scenario,
**flag_arguments,
**ratings_arguments
}))
else:
self.DER_model_list.append(
SolarPVDERThreePhase(events=self.events_list[-1],
configFile=config_file,
**{
"gridModel": self.grid_list[-1],
"identifier": scenario,
**flag_arguments,
**ratings_arguments
}))
self.sim_list.append(
DynamicSimulation(gridModel=self.grid_list[-1],
PV_model=self.DER_model_list[-1],
events=self.events_list[-1],
LOOP_MODE=False,
COLLECT_SOLUTION=True))
self.sim_list[
-1].jacFlag = False #Provide analytical Jacobian to ODE solver
self.sim_list[
-1].DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results_list.append(
SimulationResults(simulation=self.sim_list[-1], PER_UNIT=True))
def setup(self, nodeid, V0):
try:
VpuInitial = 1.0
SinglePhase = False
if 'myconfig' in OpenDSSData.config and 'DERParameters' in OpenDSSData.config[
'myconfig']:
#pvderConfig = copy.deepcopy(OpenDSSData.config['myconfig']['DERParameters'])
DERFilePath = OpenDSSData.config['myconfig']['DERFilePath']
DERModelType = OpenDSSData.config['myconfig']['DERModelType']
DERSetting = OpenDSSData.config['myconfig']['DERSetting']
DERParameters = OpenDSSData.config['myconfig']['DERParameters']
if DERSetting == 'default':
pvderConfig = self.get_derconfig(DERParameters['default'])
DERArguments = self.get_derarguments(
DERParameters['default'])
elif DERSetting == 'PVPlacement': # for manual feeder config based on PVPlacement
if nodeid in DERParameters['PVPlacement']:
pvderConfig = self.get_derconfig(
DERParameters['PVPlacement'][nodeid])
DERArguments = self.get_derarguments(
DERParameters['PVPlacement'][nodeid])
else:
raise ValueError(
'Distribution node {} not found in config file!'.
format(nodeid))
else:
raise ValueError(
'{} is not a valid DER setting in config file!'.format(
DERSetting))
if 'nodenumber' in OpenDSSData.config['myconfig']:
DERLocation = str(os.getpid()) + '-' + 'bus_' + str(
OpenDSSData.config['myconfig']
['nodenumber']) + '-' + 'node_' + nodeid
else:
DERLocation = str(os.getpid()) + '-' + 'node_' + nodeid
else:
print(
'DERParameters not found in `OpenDSSData` object - using default ratings and parameters!'
)
DERArguments = {}
pvderConfig = {}
SteadyState = True
if SinglePhase:
powerRating = 10.0e3
else:
powerRating = 50.0e3
DERLocation = 'node_' + nodeid
DERArguments.update({'pvderConfig': pvderConfig})
DERArguments.update({'powerRating': powerRating})
DERArguments.update({'SteadyState': SteadyState})
#Va = cmath.rect(DERArguments['VrmsRating']*math.sqrt(2),0.0)
#Vb = utility_functions.Ub_calc(Va)
#Vc = utility_functions.Uc_calc(Va)
a = utility_functions.Urms_calc(
V0['a'], V0['b'], V0['c']) / DERArguments['VrmsRating']
Va = (V0['a'] / a) #Convert node voltage at HV side to LV
Vb = (V0['b'] / a)
Vc = (V0['c'] / a)
DERArguments.update({'identifier': DERLocation})
DERArguments.update({'derConfig': pvderConfig})
DERArguments.update({'standAlone': False})
DERArguments.update({'gridFrequency': 2 * math.pi * 60.0})
DERArguments.update({'gridVoltagePhaseA': Va})
DERArguments.update({'gridVoltagePhaseB': Vb})
DERArguments.update({'gridVoltagePhaseC': Vc})
logging.debug(
'Creating DER instance of {} model for {} node.'.format(
DERModelType, DERLocation))
print('Creating DER instance of {} model for {} node.'.format(
DERModelType, DERLocation))
events = SimulationEvents()
PVDER_model = DERModel(modelType=DERModelType,
events=events,
configFile=DERFilePath,
**DERArguments)
self.PV_model = PVDER_model.DER_model
self.PV_model.LVRT_ENABLE = True #Disconnects PV-DER based on ride through settings in case of voltage anomaly
self.sim = DynamicSimulation(PV_model=self.PV_model,
events=events,
LOOP_MODE=True,
COLLECT_SOLUTION=True)
if DERModelType in self.sim.jac_list:
self.sim.jacFlag = True #Provide analytical Jacobian to ODE solver
else:
self.sim.jacFlag = False
self.sim.DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results = SimulationResults(simulation=self.sim,
PER_UNIT=True)
self.lastSol = copy.deepcopy(self.sim.y0) # mutable,make copy
self.lastT = 0
except Exception as e:
OpenDSSData.log("Failed Setup PVDER at node:{}!".format(nodeid))
class PVDERModel:
def __init__(self):
self._pvModel = {}
self._sim = {}
self._lastSol = {}
self._lastT = 0
def setup(self, nodeid, V0):
try:
VpuInitial = 1.0
SinglePhase = False
if 'myconfig' in OpenDSSData.config and 'DERParameters' in OpenDSSData.config[
'myconfig']:
#pvderConfig = copy.deepcopy(OpenDSSData.config['myconfig']['DERParameters'])
DERFilePath = OpenDSSData.config['myconfig']['DERFilePath']
DERModelType = OpenDSSData.config['myconfig']['DERModelType']
DERSetting = OpenDSSData.config['myconfig']['DERSetting']
DERParameters = OpenDSSData.config['myconfig']['DERParameters']
if DERSetting == 'default':
pvderConfig = self.get_derconfig(DERParameters['default'])
DERArguments = self.get_derarguments(
DERParameters['default'])
elif DERSetting == 'PVPlacement': # for manual feeder config based on PVPlacement
if nodeid in DERParameters['PVPlacement']:
pvderConfig = self.get_derconfig(
DERParameters['PVPlacement'][nodeid])
DERArguments = self.get_derarguments(
DERParameters['PVPlacement'][nodeid])
else:
raise ValueError(
'Distribution node {} not found in config file!'.
format(nodeid))
else:
raise ValueError(
'{} is not a valid DER setting in config file!'.format(
DERSetting))
if 'nodenumber' in OpenDSSData.config['myconfig']:
DERLocation = str(os.getpid()) + '-' + 'bus_' + str(
OpenDSSData.config['myconfig']
['nodenumber']) + '-' + 'node_' + nodeid
else:
DERLocation = str(os.getpid()) + '-' + 'node_' + nodeid
else:
print(
'DERParameters not found in `OpenDSSData` object - using default ratings and parameters!'
)
DERArguments = {}
pvderConfig = {}
SteadyState = True
if SinglePhase:
powerRating = 10.0e3
else:
powerRating = 50.0e3
DERLocation = 'node_' + nodeid
DERArguments.update({'pvderConfig': pvderConfig})
DERArguments.update({'powerRating': powerRating})
DERArguments.update({'SteadyState': SteadyState})
#Va = cmath.rect(DERArguments['VrmsRating']*math.sqrt(2),0.0)
#Vb = utility_functions.Ub_calc(Va)
#Vc = utility_functions.Uc_calc(Va)
a = utility_functions.Urms_calc(
V0['a'], V0['b'], V0['c']) / DERArguments['VrmsRating']
Va = (V0['a'] / a) #Convert node voltage at HV side to LV
Vb = (V0['b'] / a)
Vc = (V0['c'] / a)
DERArguments.update({'identifier': DERLocation})
DERArguments.update({'derConfig': pvderConfig})
DERArguments.update({'standAlone': False})
DERArguments.update({'gridFrequency': 2 * math.pi * 60.0})
DERArguments.update({'gridVoltagePhaseA': Va})
DERArguments.update({'gridVoltagePhaseB': Vb})
DERArguments.update({'gridVoltagePhaseC': Vc})
logging.debug(
'Creating DER instance of {} model for {} node.'.format(
DERModelType, DERLocation))
print('Creating DER instance of {} model for {} node.'.format(
DERModelType, DERLocation))
events = SimulationEvents()
PVDER_model = DERModel(modelType=DERModelType,
events=events,
configFile=DERFilePath,
**DERArguments)
self.PV_model = PVDER_model.DER_model
self.PV_model.LVRT_ENABLE = True #Disconnects PV-DER based on ride through settings in case of voltage anomaly
self.sim = DynamicSimulation(PV_model=self.PV_model,
events=events,
LOOP_MODE=True,
COLLECT_SOLUTION=True)
if DERModelType in self.sim.jac_list:
self.sim.jacFlag = True #Provide analytical Jacobian to ODE solver
else:
self.sim.jacFlag = False
self.sim.DEBUG_SOLVER = False #Check whether solver is failing to converge at any step
self.results = SimulationResults(simulation=self.sim,
PER_UNIT=True)
self.lastSol = copy.deepcopy(self.sim.y0) # mutable,make copy
self.lastT = 0
except Exception as e:
OpenDSSData.log("Failed Setup PVDER at node:{}!".format(nodeid))
def get_derconfig(self, DERParameters):
"""DER config."""
derConfig = {}
for entry in DERParameters:
if entry in templates.VRT_config_template.keys():
derConfig.update({entry: DERParameters[entry]})
return derConfig
def get_derarguments(self, DERParameters):
"""DER config."""
DERArguments = {}
for entry in specifications.DER_argument_spec:
if entry in DERParameters:
DERArguments.update({entry: DERParameters[entry]})
if 'powerRating' in DERArguments:
DERArguments.update(
{'powerRating': DERArguments['powerRating'] * 1e3})
return DERArguments
def prerun(self, gridVoltagePhaseA, gridVoltagePhaseB, gridVoltagePhaseC):
"""Prerun will set the required data in pvder model. This method should be run before
running the integrator."""
try:
# set grid voltage. This value will be kept constant during the run
# as opendss will sync will grid_simulation only once during the run call.
# opendss will solve power flow and will set gridVoltagePhase value.
self.PV_model.gridVoltagePhaseA = gridVoltagePhaseA * math.sqrt(
2) / Grid.Vbase
if templates.DER_design_template[
self.PV_model.DER_model_type]['basic_specs']['unbalanced']:
self.PV_model.gridVoltagePhaseB = gridVoltagePhaseB * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseC = gridVoltagePhaseC * math.sqrt(
2) / Grid.Vbase
self.sim.t = self.lastT + 1 / 120.0
return None
except Exception as e:
OpenDSSData.log("Failed prerun in the pvder model")
def postrun(self, sol, t, KVAbase=50):
"""Postrun will gather results. This method should be run after running the integrator."""
try:
if self.sim.COLLECT_SOLUTION:
self.sim.collect_solution(sol, t)
#np array is mutable, hence make a copy
self.lastSol = copy.deepcopy(sol[-1])
self.lastT = t[-1]
# get S
S = self._getS()
return S * KVAbase
except Exception as e:
OpenDSSData.log("Failed postrun in the pvder model")
def _run(self,
gridVoltagePhaseA,
gridVoltagePhaseB,
gridVoltagePhaseC,
KVAbase=50):
#t = [self.lastT, self.lastT + OpenDSSData.data['DNet']['DER']['DERParameters']['dt']]
t = [self.lastT, self.lastT + 1 / 120.0]
try:
# set grid voltage. This value will be kept constant during the run
# as opendss will sync will grid_simulation only once during the run call.
# opendss will solve power flow and will set gridVoltagePhase value.
self.PV_model.gridVoltagePhaseA = gridVoltagePhaseA * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseB = gridVoltagePhaseB * math.sqrt(
2) / Grid.Vbase
self.PV_model.gridVoltagePhaseC = gridVoltagePhaseC * math.sqrt(
2) / Grid.Vbase
self.sim.t = t
sol, info, ConvergeFlag = self.sim.call_ODE_solver(
self.sim.ODE_model, self.sim.jac_ODE_model, self.lastSol, t)
if self.sim.COLLECT_SOLUTION:
self.sim.collect_solution(sol, t)
#sol,info=odeint(odeFunc,y0,t,Dfun=jac_odeFunc,full_output=1,printmessg=True,hmax = 1/120.,mxstep=50,atol=1e-4,rtol=1e-4)
#sol,info=odeint(odeFunc,y0,t,full_output=1,printmessg=True,hmax = 1/120.,mxstep=50,atol=1e-4,rtol=1e-4)
#np array is mutable, hence make a copy
self.lastSol = copy.deepcopy(sol[-1])
self.lastT = t[-1]
# get S
S = self._getS()
return S * KVAbase
except Exception as e:
#OpenDSSData.log("Failed run in the pvder model")
OpenDSSData.log(
"{}:ODE solver failed between {:.4f} s and {:.4f} s!".format(
self.PV_model.name, t[0], t[-1]))
S = 0
return S
def _getS(self):
try:
return self.sim.PV_model.S_PCC # value at end of PV-DER simulation
except Exception as e:
OpenDSSData.log("Failed getS in the pvder model")