A PLC Channel Model for Home Area Networks

Smart meters (SMs) are key components of the smart grid (SG) which gather electricity usage data from residences and businesses. Home area networks (HANs) are used to support two-way communications between SMs and devices within a building such as appliances. This can be implemented using power line communications (PLCs) via home wiring topologies. In this paper, a bottom-up approach is designed and a HAN-PLC channel model is obtained for a split-phase power system which includes branch circuits, an electric panel with circuit breakers and bars, a secondary transformer and the wiring of neighboring residences. A cell division (CD) method is proposed to construct the channel model. Furthermore, arc fault circuit interrupter (AFCI) and ground fault circuit interrupter (GFCI) circuit breaker models are developed. Several HAN-PLC channels are presented and compared with those obtained using existing models.


Introduction
Increasing energy demands worldwide have resulted in significant growth in greenhouse gas emissions.As a consequence, government organizations are making efforts to combat climate change and reduce these emissions [1].One of the approaches is to modernize the power grid infrastructure and use renewable energy sources [2].These distributed energy resources require effective communications for the management of power transmission and distribution.A power grid with bidirectional communication capabilities to meet these purposes is often referred to as a smart grid (SG) [3][4][5].Smart meters (SMs) are part of the smart grid which are deployed at residences in place of traditional electric meters.
In North America and Europe, the legislative approval of smart meters has led to an increase in the number of intelligent devices such as smart appliances with a growing market of billions of dollars [6,7].One of the purposes of SMs is to help both power suppliers and residence owners with power management.Figure 1 shows the average power consumption of residences according to the survey by U.S. Department of Energy in September 2017.The appliances for heating, ventilation and air conditioning (HVAC) and water heating account for approximately half of the power consumed.A smart meter installed in a residence can provide flexibility in managing these appliances to reduce costs.
SMs enable bidirectional communications in a smart grid.They can transmit data to data acquisition centers (DACs) for analysis and receive real time pricing information from the DACs.In a home area network (HAN), SMs collect power usage data from home appliances, send commands to control these appliances, and exchange information with homeowners and tenants.Real time pricing (RTP) allows appliances such as dishwashers, washing machines and dryers to be operated during off-peak times when prices are low [8].This can also be used to reduce HVAC and water heater costs.Energy theft prevention is another function of SMs.With advanced theft detection and encryption algorithms, SMs can detect and report, as well as prevent illegal electricity usage [9].Appliance monitoring can also be provided by regularly checking and reporting the status of appliances.Power suppliers can also benefit from SMs as meter data can be obtained remotely which is more cost efficient than manual reading.Power line communications (PLC) can be employed for SM services via a home area network (HAN) [10].Since the power distribution system was not originally designed for communication purposes, the modeling of power lines as communication channels is critical to the design and analysis of HAN PLC.

Power Line Communications for Home Area Networks
The growing interest in smart home services has prompted research on broadband power line communications (BB-PLC) for high-speed data transmission.IEEE Std 1901-2010 provides networking protocols for broadband communications using power lines [11] while IEEE Std 1901.2-2013 for narrowband power line communications (NB-PLC) assures coexistence with broadband communications [12].NB-PLC is suitable for the collection of power usage data and so can be used to provide SM services within businesses and residences via HANs.Modulation schemes for PLC communications include binary phase shift keying (BPSK), frequency shift keying (FSK), spread frequency shift keying (S-FSK), spread spectrum (SS), orthogonal frequency division multiplexing (OFDM) or a combination of these techniques [13].Media access control (MAC) protocols include carrier sense multiple access (CSMA) with collision avoidance (CSMA/CA) and time division multiple access (TDMA).Single-input single-output (SISO), single-input multiple-output (SIMO), multiple-input single-output (MISO) and multiple-input multiple-output (MIMO) techniques have been developed for wireless systems, but are not included in the IEEE BB-PLC and NB-PLC standards.
Previous research on PLC-based smart metering has focused on the links between SMs and DACs [14][15][16][17].However, there has been little research on PLC for HANs.Thus, further study is required for real home wiring topologies to obtain accurate channel models.Several PLC channel models have been proposed in the literature based on experimental results for specific topologies and so cannot be generalized.In [18], a channel model was developed based on an NEC compliant topology including a panel.However, this model was constructed as a cascaded two-port network with the same path-loss in the phase and neutral conductors, which ignores the fact that they connect to different parts of the panel.In [19][20][21], a HAN-PLC model with a panel was given, but only the bonding strap to ground was considered within the panel.The model in [22] has only five branch circuits which is insufficient to model a HAN network.The panel was modeled as a black box so the results cannot be generalized either.In [23], the influence of normal and ground fault circuit interrupter (GFCI) circuit breakers on PLC was investigated which only from one part of the HAN topology.
An impedance carry back (ICB) method for PLC channel modeling was given in [24].Figure 2 shows a simple example of this method.The backbone refers to the direct path between the transmitter and receiver, while outlets o 1 , o 2 and o 3 , and nodes n 1 and n 2 are on an indirect path connected to the backbone at n b .The impedance of the path is calculated following a bottom-up approach from o 1 and o 2 to n 2 , and finally to n b as illustrated in Figure 2 from step (a) to (c).Unfortunately, the use of indirect paths in developing a channel model is difficult for large topologies with a long backbone and many indirect paths.To tackle this problem, a simple and flexible bottom-up approach that can be used for any topology is developed here.In this approach, a generalized cell division (CD) method is applied to the split-phase power system in North American residences.The contribution of this paper is in the area of power line communications (PLC) in the context of home area networks (HANs).A wiring topology is modeled for North American residences with a split-phase power supply.A smart meter is considered as the PLC transmitter or receiver which is seldom discussed in existing channel models.The parameters of the topology components are obtained and their impedances are derived.A detailed description of the panel is given for accurate modeling.Furthermore, models for arc fault circuit interrupter (AFCI) and GFCI circuit breakers are presented.These components can be replaced or the parameters modified to obtain channel models for new topologies.Compared to other channel modeling methods, cell division is simpler and can be used to model any wiring topology including those that are large and complex.A channel model obtained with the proposed method is compared with those in the literature to verify its accuracy.
The remainder of this paper is organized as follows.Section 2 gives an overview of HAN-PLC channel modeling.The cell division method is introduced in Section 3. In Sections 4 and 5, this method is used to develop HAN-PLC channel models and compared with existing modeling techniques.Finally, Section 6 provides some concluding remarks.

HAN-PLC Wiring Topology Modeling
PLC channels can be modeled using either a top-down approach or a bottom-up approach [25].Top-down models are based on parameters from extensive measurements.As wiring topologies differ, the results could not be reproduced.Furthermore, the accuracy of the measurements can have a significant effect on the performance.Conversely, a bottom-up approach uses actual topology parameters to construct the model.Since the components can easily be changed, these channel models are generic and can be adapted to any topology.This section gives an overview including the component parameters used in the bottom-up approach.

Wiring Topology of HAN-PLCs
In this paper, the wiring topology is modeled based on the National Electrical Code (NEC) and the American Wire Gauge (AWG) [26] standards for North American residences.A typical topology can be divided into three parts consisting of the topology above the SM, the electric panel up to the SM and the branch circuits.Figure 3 shows the first part of this topology.A secondary transformer delivers power to residences using a three-conductor service entrance cable (SER) [18] with AWG 4/0 conductors.From this cable, a SER with AWG 2/0 conductors is connected to the panel through the meter.The SERs above and below the SM are labeled L A and L B , respectively.SER L A is between the SM and AWG 4/0 conductors.L a 1 , L a 2 and L a 3 are the conductors in L A corresponding to phase one, neutral and phase two, respectively.SER L B is between the SM and panel and L b 1 , L b 2 and L b 3 correspond to L a 1 , L a 2 and L a 3 .The second part of the topology is the panel up to the smart meter shown in Figure 4. Two AWG 2/0 phase conductors connect the main breaker to the corresponding hot bars, where 120 V single and 240 V double-pole circuit breakers are connected to phase conductors in the branch circuits.The AWG 2/0 neutral conductor connects to a bonding strap with neutral bars at the ends.An AWG 6 bare conductor extends from the bonding strap to a ground rod, so the neutral bars have zero potential.Both the neutral and ground conductors are connected to the neutral bars.
The third part of the topology is the branch circuits which are classified as either individual, lighting or small appliance (SA) circuits.AWG 6, 8, 10, 12 and 14 conductors are used in the model corresponding to 50, 40, 30, 20 and 15 A branch circuits, respectively.An individual circuit (single or split-phase), supports one appliance which typically has high power consumption.Split-phase circuits are used for high power appliances such as a range, range top, washing machine, dryer, or water heater that are above 2000 VA.A single-phase circuit supports appliances with a relatively low average power, but large initial power may be needed to start internal motors.Lighting and SA circuits have multiple outlets.Modems are connected to outlets to enable communications and provide information on the associated devices.

Branch Circuits and the Topology Above the Panel
The per unit length (p.u.l.) resistance R, inductance L, conductance G and capacitance C of conductors are considered first.These parameters are determined by the physical properties of conductors [26] including the material (copper or aluminum), the number and diameter of conductors in cables and the number of strands, and the material and thickness of the conductor insulation.The p.u.l.resistance of a conductor is [27 where r is the radius of the conductor, σ is the conductivity, δ = 1 πµσ f is the skin depth, µ is the magnetic permeability and f is the frequency.For copper or aluminum, µ is equal to the vacuum magnetic permeability µ 0 = 4π × 10 −7 H/m.The conductivity σ is the reciprocal to the resistivity Ωm for copper, and 3.45 × 10 −8 Ωm for aluminum [26].
For simplicity, frequency dependent parameters are abbreviated such that R represents R( f ).For a pair of conductors, the inductance (in the differential mode where L in is the inner self inductance, L out is the outer self inductance, and M is the mutual inductance.If δ ≥ r, then L in = µ 8π [28].When δ < r, for a circular conductor The p.u.l.outer self inductance of a conductor is where l is the length of the conductor.The p.u.l.mutual inductance is where d is the distance between the conductors.In a four-conductor cable, two conductors can be in adjacent or diagonal positions [19] with a difference in d of √ 2. If the two conductors have equal length, and l is much greater than r and d, then ) 4 , when δ < r. ( The p.u.l.capacitance and conductance satisfy LC = µξ and LG = µσ [27].The dielectric constant is ξ = ξ 0 ξ r , where ξ 0 = 8.859 × 10 −12 F/m is the vacuum dielectric constant, and ξ r is the relative dielectric constant which is 2.3 for AWG 2/0 to 3 conductors (polyethelene) and 2.55 for AWG 4 to 14 conductors (nylon polyamide).For conductors with multiple strands (e.g., 7 strands for AWG 3 to 6, and 19 strands for AWG 2/0 and 4/0), R is multiplied by a correction factor X R [27,29].The characteristic impedance Z c and propagation constant γ are [24] where ω = 2π f .For each conductor, the end closer to the transmitter is the input, and the farther end is the output.The output impedance of a conductor Z out is determined by the impedance of all components at the output.For instance, if the conductor is between an outlet and appliance which is on, the appliance impedance is the output impedance.If the appliance is off or the outlet is open (i.e., no modem or appliance), then Z out = ∞.With N parallel impedances Z 1 , Z 2 , . . ., Z N at the output NEC recommends the length of a branch circuit should accommodate a maximum 3% voltage drop [26].For individual circuits, the voltage drop V d of a conductor is [30] where R s is the p.u.l.DC resistance of the conductor [26].The minimum length L min is 6 ft [31], and the maximum length L max should be less than 100 ft.The maximum current I max is 0.8 times the amperage rating of the corresponding circuit breaker.L max can be obtained from (7) considering the maximum voltage drop.For lighting and small appliance circuits [22] where N is the number of outlets and N max = I max /1.5 if the rated current of each outlet is 1.5 A [26].
The farthest outlet from the circuit breaker corresponds to n = 1, and l n is the length of the conductor between outlets n and n + 1, or between outlet N and the circuit breaker [18], so I n = 1.5 × n.NEC recommends that for a branch circuit, the distance from the circuit breaker to the closest outlet should not exceed 70 ft for an AWG 12 conductor (SA circuit), or 50 ft for an AWG 14 conductor (lighting circuit).The distance between outlets is 0 to 12 ft.The input impedance of a conductor is [32] and this is used as the output impedance of other conductors.The transfer function (TF) of a conductor is the ratio of the voltage at the output V o to the voltage at the input V i given by [18]

Appliance Modeling
Appliance impedances are the output impedances of the corresponding outlet conductors.The home appliances considered here are given in Table 1 [31].They can be classified as resistive, reactive or linear periodically time varying (LPTV) (types 1 to 3, respectively) [25].LPTV appliances have either commuted (3-1) or harmonic (3-2) impedance variations.There are seven types of circuits (a to g), which are split phase individual circuits, single phase individual circuits, lighting circuits, kitchen SA circuits, bedroom, study and living room (BSL) SA circuits, laundry area SA circuits, and bathroom SA circuits.P r is the power of an appliance and the impedance of the corresponding resistive load is where U = 120 V for a single phase circuit and U = 240 V for a split phase circuit.The impedance of a reactive load is obtained from the parallel RLC circuit model [25] as where R s = U 2 P f P r is the resistance at resonance.The power factor P f is between 0.7 and 1 for reactive loads.The quality factor Q f is an indication of frequency selectivity and is typically between 5 and 25.
The resonant frequency f 0 = ω 0 2π is between 25 kHz and 200 kHz [33,34].LPTV loads have impedance variations caused by non-linear elements such as thyristors which can be obtained using the approach in [25].

Secondary Transformer Modeling
In [35], impedance measurements for secondary transformers with 10 kVA to 50 kVA capacities were given.These measurements (with no cables connected), show that the real and imaginary parts of the impedance, R t and X t respectively, are proportional to the frequency.The ratio of X t to R t is approximately constant between 5 kHz and 20 kHz.Thus, the impedance can be expressed as where at 0 Hz, R t is between 0 and 1 Ω and X t is 0 Ω, and increase with frequency by (0.054, 0.081) Ω/kHz and (0.4, 1) Ω/kHz, respectively.Accurate transformer models can be obtained if the material and structure of the windings are known.

Topology Inside the Panel
The conductors in the panel are modeled differently from the rest of the topology.In branch circuits and the topology above the panel, the conductors are closely packed in cables and sealed by insulation.The length l is far greater than the cross section dimensions r and d.In the panel, the conductors are further apart and the cross section can be either circular or rectangular.The latter type comprises bars and the bonding strap, which are not sealed, and the cross section dimensions are comparable to the lengths.The p.u.l.impedance of a conductor in the panel is where R cp is the p.u.l.resistance.For circular conductors, R cp is the same as in branch circuits.For the rectangular case [36] where W and T are the width and thickness of the conductor, respectively.The imaginary part of Z cp is X cp = 2π f L cp where L cp = L in + L out is the p.u.l.inductance.The inner self inductance is where the outer self inductance L out is the same as in the circular case.The corresponding TF is The parameters of the rectangular conductors are summarized in Table 2. Circuit Breaker Modeling Thermal magnetic breakers are widely used in North America and are referred to as normal breakers.In the past 20 years, advanced AFCI or GFCI breakers have been developed to provide fault current detection and protection [37][38][39].In [23], circuit breakers with various ampere ratings were modeled, but only single-pole normal breakers were considered.In the following, both single-pole and double-pole normal and advanced circuit breakers are modeled.
Figure 5 shows the general model for a breaker.The normal breaker structure is shown on the left of the dashed line.It contains a bare copper wire, a bimetallic strip with copper and steel, and a single copper strip.For the bare copper conductor, r is determined by the breaker amperage and l is 2 in.The width and length of the strips are 0.5 and 1.75 in, respectively.In the main breaker, the thickness is 0.08 in while in branch circuit breakers it is 0.04 in.An AFCI or GFCI breaker includes the right part with coils and two conductors.One phase conductor of length 1.75 in connects the single copper strip to the corresponding branch circuit.The other is a neutral conductor of length 17.5 in which is between the branch circuit and the neutral bar.In an AFCI breaker, sensing coils T1 and T2 detect series and parallel arcing.A GFCI breaker only has a T2 coil to detect parallel arcing.A double-pole circuit breaker can be considered as two parallel single-pole breakers.For AFCI or GFCI double-pole breakers, three conductors are used and are monitored by the T1 and T2 coils.The impedance and transfer functions of the conductors within these breakers can be obtained using ( 9) and ( 10) or ( 14) and ( 17) with the parameters given above.

Bottom-Up Approach with the Cell Division Method
In this section, the cell division (CD) method is proposed as an efficient means of modeling topologies for HAN-PLC.Different from the ICB approach [24] which models indirect paths, this method divides a wiring topology into cells.Each cell is comprised of several basic components including appliances, modems and the smart meter, conductors and the secondary transformer.Each cell comprises only one junction which is a node connecting three conductors.The conductors are divided into two groups, junction-to-junction (J-J) and junction-to-unknown (J-U).A J-J conductor connects cells of two adjacent junctions.A J-U conductor is between a junction and component such as an open outlet, an appliance with a modem, a single appliance, the smart meter, or another J-U conductor.The bottom-up approach with the CD method consists of the following steps.
1. Cell division.The topology is divided into cells.2. Determine the conductor inputs and outputs.The Tx and Rx are chosen from the modems and smart meter to determine the inputs and outputs of all conductors.3. Impedance computation.For branch circuits without the Tx cell, the impedances are computed from the ends to the panel.The topology above the smart meter is modeled using the same approach.For the branch circuit with the Tx cell, the impedances are computed from the panel and the end of the circuit to the Tx cell.4. Transfer function computation.The backbone for TF computation is as follows.The cells are sorted.The Rx cell is considered first.The backbone is then included.The channel TF is then obtained.
The TF of backbone H b is then obtained as where i = 1 corresponds to the Rx cell and N c is the number of TF cells, and H J−U and H J−J are the TFs of the J-U and J-J conductors, respectively.A simple topology is shown in Figure 6.For the five cells, the impedance is computed from left to right using the bottom-up approach and only the TFs of J-U and J-J conductors named in the figure are computed.Assume that the smart meter is the Tx.Consider the topology above the panel as shown in Figure 3.If L b 1 and L b 2 are used which correspond to phase one and neutral, then the input impedance Z b of the current residence is obtained from random topologies which can be used to represent the impedances of other residences.The input impedance Z a of L A is then computed for the topology above the SM including the other residences, SER cables and secondary transformer.The TF of the topology above and including the SM is where Z sm is the impedance of the smart meter, which is the modem impedance assumed to be 50 Ω as in [25].The TF for the complete topology is then

Channel Models
To obtain a channel model for a wiring topology, a database is used which includes the characteristic impedances and propagation constants of the conductors, the cross section areas of conductors in the branch circuits, and the dimensions of the conductors in the panel.Frequencies from 3 kHz to 500 kHz are considered in the analysis [12].The frequency dependent parameters are represented by 256 discrete values uniformly distributed is this interval.In a topology, the SM is assumed to be the Tx and the Rx is chosen from modems in branch circuits.The transfer functions are determined for the channels to each of the modems.

Wiring Topology Parameters
The branch circuits are considered first and the appliances are either on or off with independent probabilities of 0.5.There are 9 individual circuits for types a and b in Table 1 with 5 double-pole breakers and 4 single-pole breakers accordingly.The power P r is uniformly distributed in the given range and used to infer the impedances of the appliances, conductors and circuit breakers.The lighting circuits are determined by the floor area.A medium size house with a floor area of 2000 ft 2 is considered.According to the NEC, the minimum lighting is 3 VA/ft 2 so 6000 VA is required.Lighting circuits are 120 V and 15 A with a maximum load of 15 × 120 × 80% = 1440 VA considering a 80% load rating for overcurrent protection.There can be 1440/180 = 8 outlets on a circuit as the outlet rating is 180 VA.Thus, 4 lighting circuits are used for general lighting and AFCI breakers recommended.Another circuit with 4 lights for 2 bathrooms is assumed and GFCI breakers recommended.Each switch controls a type c appliance such as a fluorescent or incandescent lamp with independent probabilities of 0.5.
Small appliance (SA) circuits correspond to circuit types d to g in Table 1 and each supports a maximum of 10 outlets.NEC suggests at least two circuits for kitchens, plus one for the laundry area and one for bathrooms.Thus, two kitchen SA circuits each with 10 outlets are considered, and each bathroom and laundry circuit has 4 outlets.GFCI circuit breakers are recommended for these circuits.There are 8 different appliances which can be located on the two kitchen SA circuits, 3 on the laundry SA circuit and 4 on the bathroom SA circuit.Four bedrooms, a living room and a study are assumed with 5, 10 and 5 outlets in each room, respectively, which requires 4 BSL SA circuits.There are 10 different appliances which can be located on these circuits.A TV, laptop and smart phone charger are randomly positioned in each room with independent probabilities of 0.5.One live clock is placed on a random SA outlet in the house.Furthermore, one humidifier or dehumidifier with equal probability of 0.5 is randomly located on an SA circuit.An iron is located on a laundry or BSL SA circuit and a hair dryer is located on a bathroom or BSL SA circuit.The parameters for small, medium and large house sizes including the numbers of appliances are summarized in Table 3.
Within the panel, the circuit breakers are placed evenly on both sides of the hot bars starting from the top.The length of the cable (AWG 2/0) between the panel and SM is between 6 ft and 10 ft.The distance from the SM to the cable (AWG 4/0) between residences is between 40 ft and 50 ft.From 5 and 20 residences are connected along with the secondary transformer using AWG 4/0 cable with a distance of 50 ft to 70 ft between them.

Modeling Results
The insertion loss (IL) refers to the signal loss introduced by the topology and can be expressed as To better illustrate the impact of the topology components, first only normal breakers are considered.AFCI and GFCI breaker models will be included later.MATLAB was used to obtain the transfer functions for 1000 random topologies for each of the small, medium and large house sizes.The execution times were 233 s, 322 s and 382 s, respectively.The average ILs of the TFs for each house size are given in Figure 7 and they are summarized in Table 4.Only the medium house size is considered in the remainder of this section.Figure 8 shows the ILs of 8 randomly generated topologies.To better understand the three parts of the topology (branch circuits, the panel up to the SM and above the SM), the average IL for each part is given in Figure 9 and summarized in Table 5. SA and individual circuits are now compared.The average IL of the SA and individual circuit channels is shown in Figure 10.The average IL difference is between −17.0 dB and −8.8 dB with a mean of −14.4 dB.The effect of normal and advanced circuit breakers is now examined.The average IL using only normal breakers and using AFCI and GFCI breakers as recommended by NEC [26] is shown in Figure 11.The difference is between −7.7 dB and −2.7 dB with a mean difference of −5.9 dB.

Channel Model Evaluation
The models obtained in the previous section for three house sizes with different floor areas illustrate the flexibility of the proposed approach to channel modeling.Once the channel model of a topology has been obtained, the parameters can easily be changed to obtain a new model.Figure 7 shows the average IL for the three house sizes.the IL is highest at 3 kHz but decreases rapidly with frequency and is lowest at 500 kHz.The primary reason for the IL at low frequencies is the reactive and LPTV appliances with resonant frequencies between 25 kHz and 200 kHz.At 3 kHz their impedances are small which leads to high power dissipation for PLC signals.The appliance impedances are higher between 25 kHz and 200 kHz, so the IL is lower in this region.As the frequency increases further, the skin effect of the conductors becomes a factor and results in higher losses.The IL of neighboring residences and the secondary transformer also increases.The influence of the secondary transformer is significant between 400 kHz and 500 kHz as shown in Figure 9.The average IL of the small house is approximately 5.9 dB better than the medium house and 7.6 dB better than the large house.Thus, both the frequency and topology have a significant influence on the IL, and a larger house has a higher IL.
Figure 8 shows the IL for eight randomly generated medium house size topologies.While there are variations in the ILs, some peaks and notches are similar due to the same components in the topologies.The IL of the three parts of the topology is given in Figure 9 and indicates that the topology above the SM has a significant effect on the channel.Figure 11 illustrates the effect of replacing the normal breakers with advanced breakers and shows that the advanced breakers increase the IL.This is primarily because the advanced breakers have additional neutral conductors to detect fault currents.

Channel Model Comparison
The Cañete wiring topology is a simple topology with 3 appliances as shown in Figure 12a [25].This was modeled using the ABCD matrix method in [40].The appliance parameters are given in Table 6 and the corresponding impedances |Z i | are shown in Figure 13.At junction J i , the input impedance of conductor L i is so the ABCD matrix of this topology is The transfer function H T without considering Z S is [32] With the CD method, there are three cells in the topology as shown in Figure 12b.The impedance computation starts from cell 1 and proceeds to cell 3.The TFs of conductors H J−U (1), H J−J (1), H J−J (2), H J−U (3) were obtained and their product H b was computed using (18).Figure 14 shows that this provides the same IL as H T obtained from (24).However, the proposed method is flexible in modeling complex topologies, whereas the ABCD matrix method is not efficient as it uses a single matrix to represent the entire topology which is impractical.

Conclusions
In this paper, a power line communications (PLC) channel model was developed for a home area network (HAN) in the smart grid (SG).The use of PLC for smart meter (SM) communications and the associated standards were discussed.The wiring topology of a split-phase power system was modeled in detail and the parameters of the components within the topology were provided.Cell division (CD) was proposed as a bottom-up approach to obtain a channel model for a topology.This method uses the actual parameters of the topology components and thus provides an accurate model.Performance results were presented for three topology sizes, the three parts of the topology, normal breakers and advanced breakers, and individual and small appliance (SA) circuits.These show that the CD method is an efficient means of modeling HAN topologies.The CD method was also used to develop a model for the well known Cañete wiring topology, and this was compared with another method implemented on the same topology.The results indicate that the proposed method provides accurate channel models.This method can also be used to obtain channel models for other electrical systems such as electric vehicles or solar panels, and for any building topology.
As with other bottom-up approaches to developing channel models, the parameters of the components within the wiring topology must be obtained.However, these parameters can be estimated or randomly generated to obtain a more general model.Furthermore, the proposed model is efficient and flexible so that models with different components and parameters can easily be obtained.In the future, the models obtained can be used to evaluate HAN performance with PLC considering various modulation and medium access control techniques.The proposed technique can also be used to model the noise presents in a HAN.The influence of transient events (including plug in and removal events), can also be considered in the future.

Figure 1 .
Figure 1.The average power consumption of residences.

Figure 2 .
Figure 2.An example of the impedance carry back method using three steps in (a), (b) and (c), respectively [24].

Figure 3 .
Figure 3.The topology above the panel.

Figure 4 .
Figure 4.The topology of a 200 A panel.

Figure 6 .
Figure 6.A simple topology divided into cells.

Figure 7 .
Figure 7.The average IL of the channels in the 1000 topologies for each house size.

Figure 8 .
Figure 8.The IL of 8 random channels.

Figure 9 .Table 5 .
Figure 9.The average IL of the three parts of the home topology.Table 5.The IL of the Three Parts of the Home Topology.Part Minimum (dB) Maximum (dB) Mean (dB) Branch circuits −18.8 −9.1 −15.68 Panel up to the SM −10.5 −4.5 −9.21 SM and above −66.2−27.9 −39.67

Figure 10 .
Figure 10.The average IL of the 9000 individual circuit channels and 8000 SA circuit channels in the 1000 topologies.

Figure 11 .
Figure 11.The average IL with normal breakers and AFCI and GFCI breakers.

Figure 12 .
The Cañete wiring topology modeled using (a) the ABCD matrix method, and (b) the cell division method.

Table 6 .Figure 13 .
Figure 13.The impedances of the three appliances.

Figure 14 .
Figure 14.The TFs of the Cañete wiring topology using two methods.

Table 2 .
Parameters of Rectangular Conductors in the Panel.

Table 3 .
Parameters for the Three House Sizes.

Table 4 .
The IL for the Three House Sizes.