REF: Calculation of impedance in EnergySource.

Author: juanjospina

src/data_model/transformations/ravens2math.jl

@@ -1500,6 +1500,13 @@ function _map_ravens2math_energy_source!(data_math::Dict{String,<:Any}, data_rav
             math_obj["gen_bus"] = bus_obj["bus_i"]
             data_math["bus"]["$(bus_obj["index"])"] = bus_obj
 
+            # Impedance calculation
+            r = get(ravens_obj, "EnergySource.r", 0.0)
+            r0 = get(ravens_obj, "EnergySource.r0", 0.0)
+            x = get(ravens_obj, "EnergySource.x", 0.0)
+            x0 = get(ravens_obj, "EnergySource.x0", 0.0)
+            Z_ABC = _impedance_conversion_ravens_energy_source(data_ravens, ravens_obj, r+x*1im, r0+x0*1im)
+
             branch_obj = Dict(
                 "name" => "_virtual_branch.EnergySource.$name",
                 "source_id" => "EnergySource.$name",
@@ -1511,8 +1518,8 @@ function _map_ravens2math_energy_source!(data_math::Dict{String,<:Any}, data_rav
                 "angmax" => fill(10.0, nconductors),
                 "c_rating_a" => fill(Inf, nconductors),
                 "br_status" => math_obj["gen_status"],
-                "br_r" => _impedance_conversion_ravens_energy_source(data_ravens, ravens_obj, "EnergySource.r", "EnergySource.r0"),
-                "br_x" => _impedance_conversion_ravens_energy_source(data_ravens, ravens_obj, "EnergySource.x", "EnergySource.x0"),
+                "br_r" => real(Z_ABC),
+                "br_x" => imag(Z_ABC),
                 "g_fr" => zeros(nconductors, nconductors),
                 "g_to" => zeros(nconductors, nconductors),
                 "b_fr" => zeros(nconductors, nconductors),

src/data_model/utils_ravens.jl

@@ -78,27 +78,14 @@ function _impedance_conversion_ravens(data_eng::Dict{String,Any}, eng_obj::Dict{
 end
 
 
-"converts impendance in Ohm/m by multiplying by length"
-function _impedance_conversion_ravens_energy_source(data_eng::Dict{String,Any}, eng_obj::Dict{String,Any}, key1::String, key2::String)
-    # Default energy sources considered 3 phases
-    nphases = 3
-    _impedance_matrix = zeros(Float64, nphases, nphases)
-
-    z = get(eng_obj, key1, 0.0)
-    z0 = get(eng_obj, key2, 0.0)
-
-    for i in 1:nphases
-        for j in 1:i
-            if(i==j)
-                _impedance_matrix[i, j] =  z + ((z0 - z)/3)
-            else
-                _impedance_matrix[i, j] = (z0 - z)/3
-                _impedance_matrix[j, i] = (z0 - z)/3
-            end
-        end
-    end
-
-    return _impedance_matrix .* get(eng_obj, "Conductor.length", 1.0)
+"converts impendance in Ohm/m in EnergySource"
+function _impedance_conversion_ravens_energy_source(data_eng::Dict{String,Any}, eng_obj::Dict{String,Any}, z1::Complex, z0::Complex)
+    # TODO : Single-phase
+    a = 1*exp(120*im*π/180)
+    A = [1 1 1; 1 a a^2; 1 a^2 a]
+    Z_012 = [z0 0im 0im; 0im z1 0im; 0im 0im z1]
+    Z_ABC = A^-1 * Z_012 * A
+    return Z_ABC
 end