Actually our target is to do an economic dispatch for a given network in presence of Solar PV Power
We want to do the ED for multiple periods,24hr, and we want to use AC OPF to compare the simulation results with the practical results is there any way to do this task? بتاريخ ٢٠١٧/١١/٢١ ٨:١٧ م، كتب "Mohammed Alhajri" <[email protected]>: > Thank you Ray for your clarification > > بتاريخ ٢٠١٧/١١/٢١ ٧:٠٦ م، كتب "Ray Zimmerman" <[email protected]>: > >> From the Introduction section in the MOST User’s Manual >> <http://www.pserc.cornell.edu/matpower/docs/MOST-manual-1.0.pdf> … >> >> *While the problem formulation is general and incorporates a full >> nonlinear AC network model, the current implementation is limited to DC >> power flow modeling of the network. Some work has been done on an AC >> implementation, but it is not yet ready for release.* >> >> Maybe someday … but currently, MOST is limited to a DC model of the >> network. >> >> — Ray >> >> >> >> On Nov 19, 2017, at 12:42 PM, Mohammed Alhajri <[email protected]> wrote: >> >> Hello all >> >> i wanted to compare between AC OPF results and DC OPF results using MOST >> >> My code for AC OPF is: >> >> ++++++++++++++++++++++++++++++++++++++++++++ >> >> define_constants; >> verbose = 1; >> mpopt = mpoption('verbose', verbose); >> mpopt = mpoption(mpopt, 'out.gen', 1); >> mpopt = mpoption(mpopt, 'model', 'AC'); >> mpopt = mpoption(mpopt, 'opf.ac.solver', 'MIPS'); >> mpopt = mpoption(mpopt, 'most.solver', 'MIPS'); >> if ~verbose >> mpopt = mpoption(mpopt, 'out.all', 0); >> end >> >> casefile = 'ex_case3b'; >> mpc = loadcase(casefile); >> xgd = loadxgendata('ex_xgd_ramp', mpc); >> profiles = getprofiles('ex_load_profile'); >> nt = size(profiles(1).values, 1); % number of periods >> >> mdi = loadmd(mpc, nt, xgd, [], [], profiles); >> mdo = most(mdi, mpopt); >> EPg1 = mdo.results.ExpectedDispatch; % active generation >> Elam1 = mdo.results.GenPrices; % nodal energy price >> if verbose >> most_summary(mdo); >> end >> >> >> ++++++++++++++++++++++++++++++++++ >> >> For DC OPF >> >> ++++++++++++++++++++++++++++++++++ >> >> %% set up options >> define_constants; >> verbose = 1; >> mpopt = mpoption('verbose', verbose); >> mpopt = mpoption(mpopt, 'out.gen', 1); >> mpopt = mpoption(mpopt, 'model', 'DC'); >> mpopt = mpoption(mpopt, 'opf.dc.solver', 'MIPS'); >> mpopt = mpoption(mpopt, 'most.solver', 'MIPS'); >> mpopt = mpoption(mpopt, 'most.dc_model', 1); >> if ~verbose >> mpopt = mpoption(mpopt, 'out.all', 0); >> end >> >> casefile = 'ex_case3b'; >> mpc = loadcase(casefile); >> xgd = loadxgendata('ex_xgd_ramp', mpc); >> profiles = getprofiles('ex_load_profile'); >> nt = size(profiles(1).values, 1); % number of periods >> >> %%----- Multiperiod DC OPF (w/ramp) ----- >> mdi = loadmd(mpc, nt, xgd, [], [], profiles); >> mdo = most(mdi, mpopt); >> EPg1 = mdo.results.ExpectedDispatch; % active generation >> Elam1 = mdo.results.GenPrices; % nodal energy price >> if verbose >> most_summary(mdo); >> end >> >> +++++++++++++++++++++++++++++++ >> >> >> i suppose to get in the results, isn't it? >> >> the results are attached, i don't find any difference between them! >> <DC OPF.pdf><AC OPF.pdf> >> >> >>
