# Multiresolution Time-Domain Analysis of Multiconductor Transmission Lines Terminated in Linear Loads

Authors: Zongliang Tong; Lei Sun; Ying Li; Luis Diaz Angulo; Salvador Gonzalez Garcia; Jianshu Luo
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

This paper derives a multiresolution time-domain (MRTD) scheme for the multiconductor transmission line (MTL) equations based on Daubechies’ scaling functions. The terminations are characterized by a state-variable formulation which allows a general description of the termination networks. For the linear load terminations, a method incorporating the terminal constraints is proposed to work out the scheme at and close to the terminations. The MRTD scheme is implemented with different basis functions for linear components including resistances, inductances, and capacitances. Numerical results show that the MRTD schemes obtain a more stable result than the conventional finite difference time-domain (FDTD) method with a coarse space step.


---

## Body

## 1. Introduction

The multiconductor transmission line (MTL) equations have been thought to be an efficient model for interconnected lines [1]. With the increasing of the data transmission in electric systems, transient analysis of multiconductor transmission lines (MTLs) has significant meaning in circuit design and electromagnetic compatibility (EMC). For the problems in high frequency or with wide bandwidth, the time-domain algorithm can get the solution directly and analyze the integrated response and influence quickly. Several methods have been introduced to analyze the time-domain response, including the SPICE program, the time-domain to frequency-domain (TDFD) transformation method and the finite difference time-domain (FDTD) method [2]. Since the SPICE program can obtain the exact solution for the lossless lines, it has been tried to extend to the lossy case [3–6]. Some closed-form solutions for MTLs have been developed [7–10]. This method often needs convolution integrals to transform the solution from frequency-domain to time-domain. As a direct time-domain numerical method, the FDTD method has been used to solve the MTL equations naturally [11–13]. However, this method requires a fine spatial discretization along the lines for the rapidly varying signals, which leads to large memory requirements to guarantee accuracy.For saving the memory of computing, we introduce the multiresolution time-domain (MRTD) scheme to solve the MTL equations in this paper. The MRTD scheme provides an efficient algorithm for the computation of electromagnetic fields because of its excellent capability to approximate exact solutions with lower sampling rates [14]. There are two kinds of MRTD schemes: one is to expand the variants by scaling functions in space domain, so it is called S-MRTD; the other one is to expand the variants by scaling functions and wavelet functions in space domain, so it is called W-MRTD [14]. Since the Battle-Lemarie wavelet function used in [14] is not compactly supported, the iterative equations in MRTD need to be cut off and this may lead to truncation errors. To overcome this problem, some compactly supported scaling functions have been used to formulate the MRTD schemes [15–19]. To ensure the stability of the solutions and control dispersion errors of the different schemes, the stability condition and space sampling rules for these MRTD schemes have been studied in [20–24]. Furthermore, for MTL equations, the MRTD schemes can obtain an accurate solution with a coarse space interval. The MRTD scheme has been used to solve two-conductor transmission line equations with pure resistive loads [25], and we extend this method to the MTL equations terminated with linear load networks. The linear terminal networks can be described by state-variable characteristic and allow us to model the terminal constraints directly [26, 27].In this paper, we use MRTD scheme to analyze the terminal response of MTLs with linear loads. In Section2, the MRTD scheme is derived based on Daubechies’ scaling functions for the lossy MTL equations. In Section 3, the iterative equations for terminal voltages are derived based on the linear terminal networks. The iterative equations of voltages and currents close to the terminals are modified. In Section 4, some numerical results are presented for the terminal response of both purely resistive terminations and the terminations network including inductances and capacitances. The numerical results are compared to the results obtained with SPICE and FDTD methods.


## 2. MRTD for MTL Equations

In this section, we derive the S-MRTD scheme for MTL equations. Similarly to the FDTD method, the expansions of the currents are shifted by half a discretization interval in space and time with respect to the expansions of the voltages.The homogeneous equations for(m+1)-conductor lines can be written as(1a)∂∂zVz,t+RIz,t+L∂∂tIz,t=0,(1b)∂∂zIz,t+GVz,t+C∂∂tVz,t=0,where the column vectors Vz,t and Iz,t represent the voltages and currents along the m conductors (the other one is the reference conductor) and contain m voltages Viz,t and m currents Iiz,t as shown below:(2a)Vz,t=V1z,t,…,Viz,t,…,Vmz,tT,(2b)Iz,t=I1z,t,…,Iiz,t,…,Imz,tT.m×m matrices R, L, G, and C represent the per-unit-length resistance, inductance, conductance, and capacitance, respectively.The voltage vectorV(z,t) and current vector I(z,t) can be expanded as follows:(3a)Vz,t=∑k,n=-∞+∞Vknφkzhnt,(3b)Iz,t=∑k,n=-∞+∞Ik+1/2n+1/2φk+1/2zhn+1/2t,where Vkn and Ik+1/2n+1/2 are the coefficient vectors for the voltages and currents in terms of scaling functions, respectively. Vkn and Ik+1/2n+1/2 are m dimensions vectors which are represented as follows:(4a)Vkn=V1,kn,…,Vi,kn,…,Vm,knT,(4b)Ik+1/2n+1/2=I1,k+1/2n+1/2,…,Ii,k+1/2n+1/2,…,Im,k+1/2n+1/2T.The index i represents the number of the conductors and the indices n and k are the discrete spatial and temporal index with respect to space and time coordinates via z=kΔz and t=nΔt, where Δz and Δt represent the spatial and temporal discretization intervals in z and t direction. The functions hn(t) and φk(z) are defined as(5)hnt=htΔt-n,(6)φkz=φzΔz-k.The function h(t) is defined as Haar’s scaling function and φ(z) is Daubechies’ scaling function.To derive the MRTD scheme, we need the following formulations:(7)hnt,hn′t=δn,n′Δt,hnt,∂hn′+1/2t∂t=δn,n′-δn,n′+1,where δn,n′ represents the Kronecker symbol.(8)φkz,φk′z=δk,k′Δz,φkz,∂φk′+1/2z∂z=∑i=-LSLS-1aiδk+i,k′,where LS denotes the effective support size of the basis functions. The coefficients a(i) are called connection coefficients and Table 1 lists a(i) for 0≤i≤LS-1 according to [22]. For i>LS-1,   ai=0,  and, for i<0,  ai fulfills the symmetry relation a-1-i=-ai.Table 1
Connection coefficientsa(i) of Daubechies’ scaling functions.
i
Db2
Db3
Db4
0
1.2291666667
1.2918129281
1.3110340773
1
-
0.0937500000
-
0.1371343465
-
0.1560100110
2
0.0104166667
0.0287617728
0.0419957460
3
-
0.0034701413
-
0.0086543236
4
0.0000080265
0.0008308695
5
0.0000108999
6
-
0.0000000041
q
max
⁡
0.7500
0.6844
0.6585Daubechies’ scaling functions satisfy the shifted interpolation property [28](9)φi+M1=δi,0for i integer, where M1=∫-∞+∞zϕ(z)dz is the first-moment of the scaling functions.Making use of (9), the basis function (6) is modified as [16](10)φkz=φzΔz-k+M1.In spite of the support set length of the scaling functions [29], single point sampling of the total voltages and currents can be taken at integer points with negligible error. Taking voltage at spatial point kΔz and at time nΔt, for example, we obtain(11)VkΔz,nΔt=∬-∞+∞∑k′,n′=-∞+∞Vk′n′φk′zhn′tδzΔz-kδtΔt-ndzdt=∑k′,n′=-∞+∞Vk′n′φk-k′+M1hn-n′=Vkn,where δ is the Dirac delta function. Equation (11) means the voltage value at each integer point is equal to the expansion coefficient. The current values have the same character at each half-integer point. Therefore, we use Vkn and Ik+1/2n+1/2 directly to represent the voltages at kΔz,nΔt and the currents at k+1/2Δz,n+1/2Δt in this paper.The modifiedφk(z) in (10) also satisfies integrals (8). Substituting (3a) and (3b) into (1a) and (1b) and applying the Galerkin method to (1a) and (1b), which means using hn(t)φk+1/2(z) as the test function for (1a) and hn+1/2(t)φk(z) as the test function for (1b), we can obtain the following iterative equations for the voltages and currents:(12a)Vkn+1=A1A2Vkn-ΔtΔz∑i=0LS-1aiIk+i+1/2n+1/2-Ik-i-1/2n+1/2,(12b)Ik+1/2n+1/2=B1B2Ik+1/2n-1/2-ΔtΔz∑i=0LS-1aiVk+i+1n-Vk-in,where a(i) is the connection coefficient and A1, A2, B1, and B2 are constant matrices:(13)A1=C+Δt2G-1,A2=C-Δt2G,B1=L+Δt2R-1,B2=L-Δt2R.For the stability of the solution, the temporal and spatial discretization should satisfy the Courant condition [20, 25](14)vmaxΔtΔz≤1∑i=0LS-1ai,where vmax represents the maximum propagation velocity of the modes. For MTLs in inhomogeneous media that have n velocities of propagation of the modes, the Courant condition has to be satisfied by the largest of the mode velocities, and thereby it will be satisfied by the smaller mode velocities.Note the Courant number as(15)q=vΔtΔz.The maximum values ofq required by a stable algorithm can be calculated from the connection coefficients. Table 1 lists qmax⁡ for different Daubechies’ scaling functions. For one-dimensional problems, the stability condition for the conventional FDTD method can be written as [30](16)vmaxΔtΔz≤1.Inequality (16) indicates that, for the conventional FDTD method, the maximum Courant number to guarantee the stability of the differential scheme is qmax⁡=1. In Table 1, all qmax⁡ values for MRTD schemes are smaller than 1, which means the MRTD schemes must obey a stricter stability condition compared to the conventional FDTD method. When using the same space step for MRTD schemes and FDTD method, the MRTD schemes need a smaller time step to guarantee the stability of the algorithm.


## 3. Terminal Formulations for Linear Loads

In this section, we consider the iterative equations at the source and the load. Since the iterative equations in (12a) and (12b) are not suitable for the terminations, we need to derive the iterative equations with the terminal conditions. As shown in Figure 1, we assume the length of the lines to be L, the lines are divided uniformly into NDZ segments with the space interval Δz, and the total solution time is divided into NDT steps with the uniform time interval Δt. Since we have expanded the voltages V(z,t) at the integer points and currents I(z,t) at the half-integer points in (3a) and (3b), the voltages and currents will appear alternatively along the lines as shown in Figure 1.Figure 1
Discretization of the voltages and currents.As the iterative equations (12a) and (12b) indicate, there are two parts in the iterative equations which we need to update. One part is the voltage iterative equations at the terminals. Since all the currents are located at the inner points along the lines, we do not need to analyze the currents at the terminals. The other part is the voltages and currents near the terminals in which the voltages are Vin and VNDZ-in, for i=1,2,…,LS-1, and the currents are Ii+1/2n+1/2 and INDZ-i+1/2n+1/2 for i=0,1,…,LS-2. The iterative equations of these voltages and currents contain some terms exceeding the index range.Taking voltages, for example, to derive the modified iterative equations, the modified iterative equations for currents can be obtained by the similar steps. Using the relation of the connection coefficients [17](17)∑i=0LS-12i+1ai=1substituting (17) into (12a) gives(18)∑i=0LS-12i+1aiVkn+1=A1∑i=0LS-12i+1aiA2Vkn+1+∑i=0LS-1Δt2i+1Δzai2i+1Ik+i+1/2n+1/2-Ik-i-1/2n+1/2.Considering the corresponding terms withi, we can decompose (18) as [18](19a)a0Vkn+1=A1a0A2Vkn-a0ΔtΔzIk+1/2n+1/2-Ik-1/2n+1/2,(19b)3a1Vkn+1=A13a1A2Vkn-3a1Δt3ΔzIk+3/2n+1/2-Ik-3/2n+1/2⋮(19c)2i+1aiVkn+1=A12i+1aiA2Vkn+1-2i+1aiΔt2i+1ΔzIk+i+1/2n+1/2-Ik-i-1/2n+1/2⋮(19d)2LS-1aLS-1Vkn+1=A12LS-1aLS-1A2Vkn+1-2LS-1aLS-1Δt2LS-1ΔzIk+LS-1/2n+1/2-Ik-LS+1/2n+1/2for i=0,1,…,LS-1.Equation (19c) indicates that each subequation in (19a), (19b), (19c), and (19d) is a central difference scheme of the FDTD method with the spatial interval (2i+1)Δz. Summing up all the subequations in (19a), (19b), (19c), and (19d) we can obtain (18), so we can view the MRTD schemes as the weighted mean of the conventional FDTD method, and the weighting number is (2i+1)a(i).By viewing the MRTD schemes as the weighted mean of the conventional FDTD method, (18) also indicate that the MRTD schemes contain LS differential terms in each iterative equation. For the conventional FDTD method, the iterative equations only contain 1 differential term for one-dimensional problems. Ignoring the addition terms in the iterative equations, that means the computational complexity of MRTD is LS times of conventional FDTD under the same grid numbers.We will use this relationship between FDTD method and MRTD scheme to derive the voltage iterative equations at the terminals and the voltages and currents close to the terminals. Firstly, we takeV0n+1 as an example to derive the voltage iterative equations at the terminals.DecomposingV0n+1 as (19a), (19b), (19c), and (19d), we can obtain each subequation as follows:(20)2i+1aiV0n+1=A12i+1aiA2V0n-2i+1aiΔt2i+1ΔzIi+1/2n+1/2-I-i-1/2n+1/2for i=0,1,…,LS-1.As shown in Figure1, the term I-i-1/2n+1/2 (for i=0,1,…,LS-1) in (20) does not exist, its subscript for spatial discretization exceeds the index range, and therefore the central difference scheme is not suitable for (20). Using the forward difference scheme to replace the central difference scheme, we obtain(21)2i+1aiV0n+1=A12i+1aiA2V0n-2i+1aiΔt2i+1Δz/2Ii+1/2n+1/2-I0n+1/2for i=0,1,…,LS-1.Summing up all the subequations in (21), we obtain the iterative equations at the terminals(22)V0n+1=A1A2V0n-2ΔtΔz∑i=0LS-1aiIi+1/2n+1/2-I0n+1/2.In (22), the term I0n+1/2 represents the current expansion coefficient at the point 0Δz, and since the terminal network is located at the point 0Δz, the formulation of I0n+1/2 must satisfy the constraint condition from the terminal network.In the case of a linear network, the terminal can be characterized in a state-variable form [13](23a)ddtXt=MXt+NUtwith an associated output relation(23b)Yt=OXt+PUt+QddtUt.The vectorX contains the state variables of the lumped network. These are typically the inductor currents and capacitor voltages in that network or some subset of those variables. The vectorU contains the independent sources in the network (the inputs), and the vectorY contains the designated outputs (currents and/or voltages) of the network.M,N,O,P, andQ are constant coefficient matrices, which are all determined by the specific linear network.Define the output in (23a) and (23b) as IS and IL, and then(24a)ddtXSt=MSXSt+N0V0t+NSSSt,(24b)ISt=OSXSt+P0V0t+PStSSt+Q0ddtV0t+QSddtSSt(25a)ddtXLt=MLXLt+NNDZVNDZt+NLSLt,(25b)ILt=OLXLt+PNDZVNDZt+PLtSLt+QNDZddtVNDZt+QLddtSLt,where S and L represent the networks at the source and the load.Since the variablesIS(t), XS(t), and SS(t) are located at the integer point, we can expand them as (3a) and (3b):(26)ISt=∑n=-∞+∞ISn+1/2hn+1/2t,XSt=∑n=-∞+∞XSn+1/2hn+1/2t,SSt=∑n=-∞+∞SSnhnt.Here, we expand IS(t) and XS(t) at the half-integer point to correspond with the current I0n+1/2 and this will lead to an explicit expression for the iterative equations below. V0(t) is also needed in (24a) and (24b), space sampling the voltage at the point 0Δz, and we can obtain(27)V0t=Vz,t,φ0z=∑n=-∞+∞V0nhnt.Substituting (26)–(27) into (24a) and (24b), taking hn(t) as the test function for (24a) and hn+1/2(t) as the test function for (24b), and applying the Galerkin method to (24a) and (24b), we obtain(28a)XSn+1/2-XSn-1/2Δt=MSXSn+1/2+XSn-1/22+N0V0n+NSSSn,(28b)ISn+1/2=OSXSn+1/2+P0V0n+1+V0n2+PSVSn+1+VSn2+Q0V0n+1-V0nΔt+QSSSn+1-SSnΔt.Substituting (28b) into (22)(29)V0n+1=A3A4V0n-∑i=0LS-12aiIi+1/2n+1/2+∑i=0LS-12aiOSXSn+1/2+∑i=0LS-1aiPSVSn+1+VSn+2ΔtQSSSn+1-SSn,where(30)A3=ΔzΔtC+Δz2G+∑i=0LS-1aiP0+2ΔtQ0-1,(31)A4=ΔzΔtC-Δz2G-∑i=0LS-1aiP0-2ΔtQ0,(32)XSn+1/2=1m-Δt2MS-11m+Δt2MS+ΔtN0V0n+NSSSn,where 1m represents m dimension unity matrix.Note that the expression ofXS in (32) is a little different from the iterative equations in [2, 13]. The expression of state variable XSn+1/2 in this paper falls behind the terminal voltage V0n+1 half a temporal step, while in [2, 13] the state variable and the terminal voltage are at the same time point. So the state variable and the terminal voltage are coupling in [2, 13], and there is a matrix inversion process to obtain the state variables and the terminal voltages. On the other hand, these two variables are decoupled in (28a) and (28b). That means we separate the iterative equations of the state variables and the terminal voltages, and this separation reduces the matrix inversion process to obtain the state variables and the terminal voltages.Following a similar procedure, we obtain(33)VNDZn+1=A5A6VNDZn+∑i=0LS-12aiINDZ-i-1/2n+1/2-∑i=0LS-12aiOLXLn+1/2-∑i=0LS-1aiPLVLn+1+VLn+2ΔtQLSLn+1-SLn,where(34)A5=ΔzΔtC+Δz2G+∑i=0LS-1aiPNDZ+2ΔtQNDZ-1,A6=ΔzΔtC-Δz2G-∑i=0LS-1aiPNDZ-2ΔtQNDZ,XLn+1/2=1m-Δt2ML-11m+Δt2ML+ΔtNNDZVNDZn+NLSLn.In the case of resistance terminations, the generalized Thevenin equivalent representation in (24a) and (24b) is written as [2, 13](35)IS=-GSV0+GSVS,IL=GLVNDZ+GLVL,where GS=RS-1 and GL=RL-1.Hence,XS=XL=0, MS=N0=NS=OS=Q0=QS=0, ML=NNDZ=NL=OL=QNDZ=QL=0, P0=GS, PS=GS, PNDZ=GL, and PL=-GL.The iterative equations of the terminal voltages can be reduced as(36)V0n+1=A7A8V0n-2RS∑i=0LS-1aiI1/2n+1/2+VSn+1+VSn,VNDZn+1=A9A10VNDZn+2RL∑i=0LS-1aiINDZ-1/2n+1/2+VLn+1+VLn,where(37)A7=ΔzΔtRSC+Δz2RSG+∑i=0LS-1ai1m-1,A8=ΔzΔtRSC-Δz2RSG-∑i=0LS-1ai1m,A9=ΔzΔtRLC+Δz2RLG+∑i=0LS-1ai1m-1,A10=ΔzΔtRLC-Δz2RLG-∑i=0LS-1ai1m.For the iterative equations of voltages and currents near the terminals, we will use the truncation method to remove the terms whose index exceeds the range. This method is based on the point to view the MRTD scheme as the weighted mean of the conventional FDTD method. TakingVkn+1 as an example, for k=1,2,…,LS-1, decomposing the iterative equation of Vkn+1 as (19a), (19b), (19c), and (19d), each subequation can be written as(38)2i+1aiVkn+1=A12i+1aiA2Vkn+1-2i+1aiΔt2i+1ΔzIk+i+1/2n+1/2-Ik-i-1/2n+1/2for k=1,2,…,LS-1.Notice that ifk>i, there does not exist a term exceeding the index range in subequation (38). That means the first k terms of the subequations in (38) do not contain the exceeding index range terms, while the terms exceeding the index range all appear in the last LS-k subequations. By viewing MRTD schemes as the weighted mean of the conventional FDTD method, the last LS-k subequations are suitable for forming the iterative equations in MRTD schemes. Using the part summation to approximate the whole summation, we update the iterative equation of Vkn+1 by using the first k subequations. The iterative equation becomes(39)Vkn+1=A1A2Vkn-∑i=0k-12i+1ai-1ΔtΔz∑i=0k-1aiIk+i+1/2n+1/2-Ik-i-1/2n+1/2for k=1,2,…,LS-1.Using the same method, we obtain the iterative equations of voltages and currents close to the terminals. The voltage iterative equations close to the load are(40)Vkn+1=A1A2Vkn-∑i=0NDZ-k-12i+1ai-1ΔtΔz·∑i=0NDZ-k-1aiIk+i+1/2n+1/2-Ik-i-1/2n+1/2for k=NDZ-LS+1,NDZ-LS+2,…,NDZ-1.The current iterative equations close to the source are(41)Ik+1/2n+1/2=B1B2Ik+1/2n-1/2-∑i=0k2i+1ai-1ΔtΔz∑i=0kaiVk+i+1n-Vk-infor k=0,1,…,LS-2.The current iterative equations close to the load are(42)Ik+1/2n+1/2=B1B2Ik+1/2n-1/2-∑i=0NDZ-k-12i+1ai-1ΔtΔz·∑i=0NDZ-k-1aiVk+i+1n-Vk-infor k=NDZ-LS+1,NDZ-LS+2,…,NDZ-1.


## 4. Numerical Results

### 4.1. Ribbon Cables

In this section, we consider a three-wire ribbon cable as shown in Figure2. The total length of the line is 2 m. One wire is chosen as the reference conductor and noted as Wire 0, and other wires are noted as Wire 1 and Wire 2, respectively. The source is located at the left side of Wire 2, and the source resistance is 50 Ω. The left side of Wire 1 is a 50 Ω resistance. The right side of Wire 1 and Wire 2 are both resistances with 50 Ω, and this model comes from [2], as shown in Figure 2(a).Figure 2
A three-wire ribbon cable.
(a)
Line dimensions and terminals
(b)
Cross-sectional dimensionsThe cross section of the three-wire ribbon cable is shown in Figure2(b). The center-to-center separation of the wires is d=50 mils (1 mil = 2.54×10-5 m). The wires are identical and are composed of wires with radius of rw=7.5 mils and polyvinyl chloride (PVC) insulations of thickness t=10 mils and relative dielectric constant εr=3.5. Following the method in [2], the per-unit-length inductance matrix L and capacitance matrix C are computed as(43)L=0.74850.50770.50771.0154μH/mC=37.432-18.716-18.71624.982pF/m.UsingL and C matrices above we can compute the mode velocities of propagation. They are vm1=2.324×108 m/s and vm2=2.5106×108 m/s giving one-way time delays of the modes of TDm1=8.606 ns and TDm2=7.996 ns. The characteristic impedances of the modes are ZCm1=91.96 Ω and ZCm2=254.93 Ω.We use a ramp function as the source shown in Figure3, the rising time of the ramp function is τr=2 ns, and the total time is T=200 ns. The amplitude is 1 V and the initial value of the ramp function is 0 V. The bandwidth of the source is approximate to BW=1/τr=500 MHz, so the highest frequency of the incident wave is approximately equal to 500 MHz. Taking 1/10 of the wavelength of the highest frequency of the incident wave as the segment length, the segment length Δz is 4.468×10-2 m. This space step is often suitable for FDTD method, so we use this segment length for numerical computing first. Since we divide the line into NDZ segments uniformly, the space discretization number is NDZ = 45. For the stability of the algorithm, taking q=0.5, we can compute(44)NDT=NDZ×v×Tq×L.Figure 3
Waveform of the source voltage.In Figures4–12, we use Dbi-MRTD to represent the computed results of MRTD scheme by using Dbi scaling functions as basis functions. Since the SPICE program can compute the response exactly for the lossless lines and the FDTD method has been proved to be an effective algorithm for the time-domain response, we compare the computed results of MRTD scheme with SPICE and FDTD methods to validate the effectiveness of the MRTD schemes.Figure 4
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and Db2-MRTD. The space discretization number NDZ = 45 and Courant numberq=0.5.Figure 5
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and MRTD: (a) for Db3-MRTD and (b) for Db4-MRTD. The space discretization number NDZ = 45 and Courant numberq=0.5.
(a)
(b)Figure 6
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and MRTD: (a) for Db3-MRTD and (b) for Db4-MRTD. The space discretization number NDZ = 12 and Courant numberq=0.5.
(a)
(b)Figure 7
A printed circuit board.
(a)
Cross-sectional dimensions
(b)
Line dimensions and terminalsFigure 8
Waveform of a trapezoid pulse.Figure 9
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 398 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 10
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 795 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 11
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 199 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 12
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 398, and the Courant numbers areqFDTD=1, qDb3-MRTD=0.6844, and qDb4-MRTD=0.6585 as shown in Table 1.
(a)
Voltages atS1
(b)
Voltages atL1Figure4 shows the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and Db2-MRTD. It indicates that Db2-MRTD does not agree with the computed results of the SPICE program and the FDTD method. The reason for this phenomenon is that the scaling function of the Db2-wavelet does not have a high enough vanishing moment. When we use Daubechies’ scaling functions to expand the voltages and currents in the MTL equations, the vanishing moment decides the accuracy of the approximation. When using scaling functions with low vanishing moments as basis functions to approximate voltages and currents, such as Db2-wavelet’s scaling functions, the approximation may introduce a larger error, while when using the scaling functions with high vanishing moments as basis functions, the approximation can be more accurate.Figure5 shows computed results of Db3-MRTD and Db4-MRTD. Both the two figures show the corresponding computing results by SPICE program, FDTD method, and MRTD schemes. The results of FDTD method show some oscillation on the edge of transition, while the computing results of Db3-MRTD and Db4-MRTD schemes show a better solution on the edge of transition.To validate the advantage of the MRTD schemes in space sampling, we use a coarse space interval to compute the terminal response. Taking the segment lengthΔz to be Δz=λ/4, the discretization will become NDZ = 12.Figure6 shows that both the computed results of FDTD method and MRTD schemes exhibit some oscillation on the edge of transition with the coarse space interval. The oscillation for the MRTD schemes appears almost at the first “stage” with a smaller amplitude; however, the oscillation for the FDTD method appears at several “stages” of the computed results with a larger amplitude. This indicates that we can use the MRTD schemes to obtain more accurate results than the conventional FDTD method when using the same grid.Table2 shows the runtime for different schemes with different space discretization numbers. The numerical results in Table 2 show that the runtime of Db3-MRTD is almost 5 times that of FDTD method, and the runtime of Db4-MRTD is almost 7 times that of FDTD method. Since the effective support sizes of Db3’s and Db5’s scaling functions are 5 and 7 as shown in Table 1, the times of the runtime between MRTD and FDTD are corresponding to the size of the effective support of the basis functions. These results agree with the analysis in Section 3.Table 2
Runtime for different schemes with different space discretization numbers.
Different schemes
Runtime for NDZ = 12
Runtime for NDZ = 45
FDTD
0.0239 s
0.1085 s
Db3-MRTD
0.1146 s
0.4693 s
Db4-MRTD
0.1855 s
0.7688 sThe numerical results in Figures5 and 6 and Table 2 show that the MRTD schemes need more time to obtain more stable results when using the same space step with the conventional FDTD method. This indicates that the MRTD schemes allow obtaining more accurate results by sacrificing the computing efficiency.


### 4.2. Printed Circuit Boards

In this section, we compute the terminal response for a lossy printed circuit board. Three conductors of rectangular cross section of widthw=5 mils and thickness t=0.5 mils are separated by s=5 mils and placed on one side of a silicon substrate having εr=12 and thickness h=5 mils [13], as shown in Figure 7(a), noting the reference as Wire 0, and the other two lines are noted as Wire 1 and Wire 2.Following the method in [2], the per-unit-length inductance and capacitance matrix are computed as (45)L=0.8057560.5387710.5387711.07754μH/mC=117.791-58.8956-58.895671.8544pF/m.These give mode velocities in the lossless case of vm1=1.25809×108 m/s and vm2=1.47934×108 m/s.Ignoring the dielectric loss, we can obtainG=0. The per-unit-length DC resistance of each conductor is computed as rdc=1/(σwt)=10.6897 Ω/m; just considering the DC resistance, we obtain(46)R=10.78970010.7897Ω/m.The total length of the line is50 cm and the line is terminated with networks as shown in Figure 7(b). These networks at the source and load can be described by state variables as mentioned in Section 3. The state variables in (24a) and (24b) and (25a) and (25b) are listed below as [13](47)XSt=0,MS=0,N1=0,NS=0,OS=0,P1=1RNE001RS,PS=1RNE001RS,Q1=0,QS=0,SSt=0VSt,XLt=iLt0,ML=-LR0,NNDZ=1L0,NL=0,OL=10,PNDZ=0001RL,PL=0,QNDZ=C000,QL=0,where C=100 pF, L=1 μH, R=10 Ω, RS=50 Ω, RL=50 Ω, and RNE=5 Ω and iLt is the inductor current.The source is a trapezoidal pulse shown in Figure8, and the amplitude is 1 V. The rise time/fall time of the pulse are τr=τf=100 ps and pulse width is τ=500 ps. The total computing time is 40 ns. The bandwidth of the input is approximately equal to 1/τr=10 GHz, taking the segment length as Δz=λ/10. Using the small mode velocity vm1 at 10 GHz, we can get NDZ=398. Choosing the Courant number q satisfying the stability condition, the time discretization number NDT can be obtained from (44).Figures9(a) and 9(b) show terminal voltages of Wire 1 at the source and load by the FDTD method and the MRTD schemes. Since Db2-MRTD scheme is not an effective algorithm for MTL equations, we will just present the computed results of Db3-MRTD and Db4-MRTD. Both Figures 9(a) and 9(b) show the agreement of these two kinds of method. It can be seen from Figure 9(a) that there exist more oscillations on the edge of each transition from the computed result of FDTD method, which means the FDTD method does not perform well with the rapid signals. On the other hand, the computed results also show that the MRTD schemes could lead to a more stable result for the rapid signals. And, for the gradually changing signals, both the two methods show good performance as shown in Figure 9(b).To obtain a more accurate result, we could use a smaller space step and time step. Figures10(a) and 10(b) show the computed results when we increase the space discretization number NDZ to 795. Both the figures show fewer oscillations on the edge of transition.Figures11(a) and 11(b) show the computed results when we decrease the space discretization step to Δz=λ/5. Figure 11(a) also indicates that the MRTD schemes could obtain a more stable solution with a coarse space step compared to the conventional FDTD method.Table3 shows the runtime for Figures 9–11. The numerical results show that MRTD schemes need more runtime to obtain more stable results. The runtime of Db3-MRTD is almost 5 times that of FDTD method, and the result is corresponding to the analysis in Section 3. The runtime of Db4-MRTD is almost 6 times that of FDTD method, while the effective support size of the Db4’s scaling functions is 7. The runtime of Db4-MRTD is less than the effective support size of the Db4’s scaling functions. The numerical results in Figures 9–11 and Table 3 also indicate that the MRTD schemes allow obtaining more accurate results by sacrificing the computing efficiency.Table 3
Runtime for different schemes with different space discretization numbers.
Different schemes
Runtime for NDZ = 199
Runtime for NDZ = 398
Runtime for NDZ = 795
FDTD
0.2060 s
0.5999 s
2.0181 s
Db3-MRTD
0.9581 s
2.7241 s
8.7192 s
Db4-MRTD
1.2724 s
3.6920 s
11.9143 sThe oscillations may be caused by the oversampling of the time step, so we compute the terminal voltages by different Courant numbers. For FDTD, we takeq=1 and, for MRTD schemes, we take the maximum Courant number for each MRTD scheme. Figure 12 shows the results. Since the MTLs have different mode velocity in each wire, the “magic step” for one dimension does not exist, and when we take Δt=Δz/vmax for the FDTD method, the oscillations do not disappear. Figure 12(a) shows that there exist more oscillations for the computing results of the FDTD methods with different Courant numbers. It means that the MRTD schemes are more stable than the FDTD methods for rapid signals.


## 4.1. Ribbon Cables

In this section, we consider a three-wire ribbon cable as shown in Figure2. The total length of the line is 2 m. One wire is chosen as the reference conductor and noted as Wire 0, and other wires are noted as Wire 1 and Wire 2, respectively. The source is located at the left side of Wire 2, and the source resistance is 50 Ω. The left side of Wire 1 is a 50 Ω resistance. The right side of Wire 1 and Wire 2 are both resistances with 50 Ω, and this model comes from [2], as shown in Figure 2(a).Figure 2
A three-wire ribbon cable.
(a)
Line dimensions and terminals
(b)
Cross-sectional dimensionsThe cross section of the three-wire ribbon cable is shown in Figure2(b). The center-to-center separation of the wires is d=50 mils (1 mil = 2.54×10-5 m). The wires are identical and are composed of wires with radius of rw=7.5 mils and polyvinyl chloride (PVC) insulations of thickness t=10 mils and relative dielectric constant εr=3.5. Following the method in [2], the per-unit-length inductance matrix L and capacitance matrix C are computed as(43)L=0.74850.50770.50771.0154μH/mC=37.432-18.716-18.71624.982pF/m.UsingL and C matrices above we can compute the mode velocities of propagation. They are vm1=2.324×108 m/s and vm2=2.5106×108 m/s giving one-way time delays of the modes of TDm1=8.606 ns and TDm2=7.996 ns. The characteristic impedances of the modes are ZCm1=91.96 Ω and ZCm2=254.93 Ω.We use a ramp function as the source shown in Figure3, the rising time of the ramp function is τr=2 ns, and the total time is T=200 ns. The amplitude is 1 V and the initial value of the ramp function is 0 V. The bandwidth of the source is approximate to BW=1/τr=500 MHz, so the highest frequency of the incident wave is approximately equal to 500 MHz. Taking 1/10 of the wavelength of the highest frequency of the incident wave as the segment length, the segment length Δz is 4.468×10-2 m. This space step is often suitable for FDTD method, so we use this segment length for numerical computing first. Since we divide the line into NDZ segments uniformly, the space discretization number is NDZ = 45. For the stability of the algorithm, taking q=0.5, we can compute(44)NDT=NDZ×v×Tq×L.Figure 3
Waveform of the source voltage.In Figures4–12, we use Dbi-MRTD to represent the computed results of MRTD scheme by using Dbi scaling functions as basis functions. Since the SPICE program can compute the response exactly for the lossless lines and the FDTD method has been proved to be an effective algorithm for the time-domain response, we compare the computed results of MRTD scheme with SPICE and FDTD methods to validate the effectiveness of the MRTD schemes.Figure 4
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and Db2-MRTD. The space discretization number NDZ = 45 and Courant numberq=0.5.Figure 5
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and MRTD: (a) for Db3-MRTD and (b) for Db4-MRTD. The space discretization number NDZ = 45 and Courant numberq=0.5.
(a)
(b)Figure 6
Comparison of the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and MRTD: (a) for Db3-MRTD and (b) for Db4-MRTD. The space discretization number NDZ = 12 and Courant numberq=0.5.
(a)
(b)Figure 7
A printed circuit board.
(a)
Cross-sectional dimensions
(b)
Line dimensions and terminalsFigure 8
Waveform of a trapezoid pulse.Figure 9
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 398 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 10
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 795 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 11
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 199 and the Courant number isq=0.5.
(a)
Voltages atS1
(b)
Voltages atL1Figure 12
Computing results of terminal voltages for the lossy PCB by FDTD method and MRTD schemes. The space discretization number is NDZ = 398, and the Courant numbers areqFDTD=1, qDb3-MRTD=0.6844, and qDb4-MRTD=0.6585 as shown in Table 1.
(a)
Voltages atS1
(b)
Voltages atL1Figure4 shows the time-domain response of the left-end crosstalk of the ribbon cable by SPICE program, FDTD, and Db2-MRTD. It indicates that Db2-MRTD does not agree with the computed results of the SPICE program and the FDTD method. The reason for this phenomenon is that the scaling function of the Db2-wavelet does not have a high enough vanishing moment. When we use Daubechies’ scaling functions to expand the voltages and currents in the MTL equations, the vanishing moment decides the accuracy of the approximation. When using scaling functions with low vanishing moments as basis functions to approximate voltages and currents, such as Db2-wavelet’s scaling functions, the approximation may introduce a larger error, while when using the scaling functions with high vanishing moments as basis functions, the approximation can be more accurate.Figure5 shows computed results of Db3-MRTD and Db4-MRTD. Both the two figures show the corresponding computing results by SPICE program, FDTD method, and MRTD schemes. The results of FDTD method show some oscillation on the edge of transition, while the computing results of Db3-MRTD and Db4-MRTD schemes show a better solution on the edge of transition.To validate the advantage of the MRTD schemes in space sampling, we use a coarse space interval to compute the terminal response. Taking the segment lengthΔz to be Δz=λ/4, the discretization will become NDZ = 12.Figure6 shows that both the computed results of FDTD method and MRTD schemes exhibit some oscillation on the edge of transition with the coarse space interval. The oscillation for the MRTD schemes appears almost at the first “stage” with a smaller amplitude; however, the oscillation for the FDTD method appears at several “stages” of the computed results with a larger amplitude. This indicates that we can use the MRTD schemes to obtain more accurate results than the conventional FDTD method when using the same grid.Table2 shows the runtime for different schemes with different space discretization numbers. The numerical results in Table 2 show that the runtime of Db3-MRTD is almost 5 times that of FDTD method, and the runtime of Db4-MRTD is almost 7 times that of FDTD method. Since the effective support sizes of Db3’s and Db5’s scaling functions are 5 and 7 as shown in Table 1, the times of the runtime between MRTD and FDTD are corresponding to the size of the effective support of the basis functions. These results agree with the analysis in Section 3.Table 2
Runtime for different schemes with different space discretization numbers.
Different schemes
Runtime for NDZ = 12
Runtime for NDZ = 45
FDTD
0.0239 s
0.1085 s
Db3-MRTD
0.1146 s
0.4693 s
Db4-MRTD
0.1855 s
0.7688 sThe numerical results in Figures5 and 6 and Table 2 show that the MRTD schemes need more time to obtain more stable results when using the same space step with the conventional FDTD method. This indicates that the MRTD schemes allow obtaining more accurate results by sacrificing the computing efficiency.


## 4.2. Printed Circuit Boards

In this section, we compute the terminal response for a lossy printed circuit board. Three conductors of rectangular cross section of widthw=5 mils and thickness t=0.5 mils are separated by s=5 mils and placed on one side of a silicon substrate having εr=12 and thickness h=5 mils [13], as shown in Figure 7(a), noting the reference as Wire 0, and the other two lines are noted as Wire 1 and Wire 2.Following the method in [2], the per-unit-length inductance and capacitance matrix are computed as (45)L=0.8057560.5387710.5387711.07754μH/mC=117.791-58.8956-58.895671.8544pF/m.These give mode velocities in the lossless case of vm1=1.25809×108 m/s and vm2=1.47934×108 m/s.Ignoring the dielectric loss, we can obtainG=0. The per-unit-length DC resistance of each conductor is computed as rdc=1/(σwt)=10.6897 Ω/m; just considering the DC resistance, we obtain(46)R=10.78970010.7897Ω/m.The total length of the line is50 cm and the line is terminated with networks as shown in Figure 7(b). These networks at the source and load can be described by state variables as mentioned in Section 3. The state variables in (24a) and (24b) and (25a) and (25b) are listed below as [13](47)XSt=0,MS=0,N1=0,NS=0,OS=0,P1=1RNE001RS,PS=1RNE001RS,Q1=0,QS=0,SSt=0VSt,XLt=iLt0,ML=-LR0,NNDZ=1L0,NL=0,OL=10,PNDZ=0001RL,PL=0,QNDZ=C000,QL=0,where C=100 pF, L=1 μH, R=10 Ω, RS=50 Ω, RL=50 Ω, and RNE=5 Ω and iLt is the inductor current.The source is a trapezoidal pulse shown in Figure8, and the amplitude is 1 V. The rise time/fall time of the pulse are τr=τf=100 ps and pulse width is τ=500 ps. The total computing time is 40 ns. The bandwidth of the input is approximately equal to 1/τr=10 GHz, taking the segment length as Δz=λ/10. Using the small mode velocity vm1 at 10 GHz, we can get NDZ=398. Choosing the Courant number q satisfying the stability condition, the time discretization number NDT can be obtained from (44).Figures9(a) and 9(b) show terminal voltages of Wire 1 at the source and load by the FDTD method and the MRTD schemes. Since Db2-MRTD scheme is not an effective algorithm for MTL equations, we will just present the computed results of Db3-MRTD and Db4-MRTD. Both Figures 9(a) and 9(b) show the agreement of these two kinds of method. It can be seen from Figure 9(a) that there exist more oscillations on the edge of each transition from the computed result of FDTD method, which means the FDTD method does not perform well with the rapid signals. On the other hand, the computed results also show that the MRTD schemes could lead to a more stable result for the rapid signals. And, for the gradually changing signals, both the two methods show good performance as shown in Figure 9(b).To obtain a more accurate result, we could use a smaller space step and time step. Figures10(a) and 10(b) show the computed results when we increase the space discretization number NDZ to 795. Both the figures show fewer oscillations on the edge of transition.Figures11(a) and 11(b) show the computed results when we decrease the space discretization step to Δz=λ/5. Figure 11(a) also indicates that the MRTD schemes could obtain a more stable solution with a coarse space step compared to the conventional FDTD method.Table3 shows the runtime for Figures 9–11. The numerical results show that MRTD schemes need more runtime to obtain more stable results. The runtime of Db3-MRTD is almost 5 times that of FDTD method, and the result is corresponding to the analysis in Section 3. The runtime of Db4-MRTD is almost 6 times that of FDTD method, while the effective support size of the Db4’s scaling functions is 7. The runtime of Db4-MRTD is less than the effective support size of the Db4’s scaling functions. The numerical results in Figures 9–11 and Table 3 also indicate that the MRTD schemes allow obtaining more accurate results by sacrificing the computing efficiency.Table 3
Runtime for different schemes with different space discretization numbers.
Different schemes
Runtime for NDZ = 199
Runtime for NDZ = 398
Runtime for NDZ = 795
FDTD
0.2060 s
0.5999 s
2.0181 s
Db3-MRTD
0.9581 s
2.7241 s
8.7192 s
Db4-MRTD
1.2724 s
3.6920 s
11.9143 sThe oscillations may be caused by the oversampling of the time step, so we compute the terminal voltages by different Courant numbers. For FDTD, we takeq=1 and, for MRTD schemes, we take the maximum Courant number for each MRTD scheme. Figure 12 shows the results. Since the MTLs have different mode velocity in each wire, the “magic step” for one dimension does not exist, and when we take Δt=Δz/vmax for the FDTD method, the oscillations do not disappear. Figure 12(a) shows that there exist more oscillations for the computing results of the FDTD methods with different Courant numbers. It means that the MRTD schemes are more stable than the FDTD methods for rapid signals.


## 5. Conclusion

In this paper, we derived the MRTD scheme for multiconductor transmission lines with linear terminations. By using state variant to describe the linear terminations, a method of incorporating the terminal constraints into the MRTD iterative equations is proposed to update the iterative equations which contain some terms whose indexes exceed the index range in MRTD schemes. Using different scaling functions as basis functions, the MRTD schemes are implemented for pure resistance terminations and linear network terminations including inductances and capacitances. The numerical results of MRTD schemes are compared to SPICE model and conventional FDTD, and the MRTD schemes are proved to be an effective algorithm for MTL equations and show a more stable performance with a coarse space step. However, the numerical results also show that the MRTD schemes need more time to obtain more stable results when compared to the conventional FDTD method. It indicates that the MRTD schemes allow obtaining more accurate results by sacrificing the computing efficiency.



---

*Source: 9845702-2017-03-20.xml*