Design of Wideband 3-way Power Divider/Combiner

with 3D FDTD

Abstract— An optimum topology design of wideband power A wide variety of power dividers/combiners with high
divider/combiner has been developed. The performance is very performance have been modeled in different topologies, with
good, with high transmission and low insertion loss from dc to 20 various methods. An overview of these topologies can be
GHz. 3D-FDTD Method with implementation of Uniaxial
Perfectly Matched Layer (UPML) absorbing boundary conditions
found in [2,3].
is used to simulate the structure. A numerical optimization was In the present paper, a three-way power divider/combiner
successfully used to determine the optimal topology structure for based on microstrip circuit has been analysed, modelled and
design 3-way power divider/combiner and an improvement was simulated with a great performance characteristic and an
achieved compared with some other results of recent publications. optimum topology, the designed power divider/combiner is
The Simulation results obtained in this article show a great operating at ultra wideband from dc to 20 GHz.
performance characteristic improvement of the 3- way power
divider/combiner, the return loss is greater than - 10.0 dB from Several computational techniques were used for the analysis
dc to 20, and better than -14 dB at ku band. The proposed of power dividers/combiners structures. These techniques can
structure has more application in microwave systems which be classified as: finite element method, moment method, finite
require amplifiers operating at higher frequencies and providing difference time domain method. A three-dimensional finite-
high power such as radar or satellite applications. difference time-domain (3D-FDTD) method is used to
simulate the electromagnetic (EM) waves propagation in the
Index Terms— Power Divider/Combiner, Microstrip, FDTD,
proposed power divider, it is simple, easy to implement and a
Topology, UPML.
very efficient solution for wideband problems [4,5] and
Paper Classification particularly well suited to analysis of microwave circuits such
DOI: 10.7251/ELSxxxxxxxx as mixers, multiplexer, power amplifiers, antenna array, power
combiners and dividers, and it has an advantage to give a
broadband solution with a single simulation run.
I. INTRODUCTION Boundary condition for reducing boundary reflections
should be introduced to simulate the proposed structure. The
P OWER dividers/combiners are essential passive devices
for many applications in telecommunication, radars,
satellites, and other systems, such as array antenna power
most used techniques to truncate the limit for open region
modeling problems are the Mur's ABC artificial boundary
amplifiers. They are widely used in applications that require condition, the Liao ABC [5, 6] and various perfectly matched
power higher than the capacity of a single electromagnetic layer PML formulations [7]. The PML concept was initially
source to increase their output power [1,2,3]. They can proposed by Jean Pierre Berenger in 1994 and has been
combine electromagnetic waves from multiple sources to a modified and extended to the uniaxial PML (UPML) [8,9], it
single output port, and can also be used as power dividers have the ability to absorb evanescent waves, more
when used in reverse. The key features of power computationally efficient and is proved to be stable with
divider/combiner are clearly their insertion loss, return loss, greater flexibility [9,11]. In our work, uniaxial perfectly
bandwidth and size. The techniques adopted for the matching layer (UPML) in three dimensions is used as the
optimization of device performance are strictly related to the Mur's ABC to terminate the finite difference time domain
type of power divider/combiner structure and to the selected formulation for electromagnetic wave, it is efficient and can be
realization technology [1,2]. easily implemented with 3D- FDTD method and it will absorb
unwanted reflections and truncating unbounded domains in
numerical simulations.
This paper is divided into the following sections: In the first,
the theoretical analysis and design configurations of the
proposed power divider/combiner will be introduced. In
section II the simulation algorithm and numerical optimization
will be explained, next we give the results obtained and
provide a comparison with other results of recent publications.
Finally, the paper is concluded in section III.
 II. ANALYSIS AND DESIGN OF POWER DIVIDER/COMBINER Where sx ,sy and sz are defined by [6, 7]:
The analysis consists of using a theory to provide complete
explanations with the implementation of a numerical method. x y 
sx  k x  ; sy  ky  ; sz  k z  z (5 a,b,c)
Starting with mathematical proof, we explain how to calculate j j j
the scattering parameters. In first we give a formulation to
calculate the electromagnetic fields that will be used to One approach is to define the constitutive relationship to
calculate the Sij coefficients, then we present the design decouple the frequency dependent terms [6,7], specifically, let:
procedure with a numerical optimization. ~ s ~ ~ s ~ ~ sy ~
Dx   z E x ; D y   x E y ; Dz   Ez (6 a,b,c)
A. FDTD with UPM Formulation sy sy sz
UPML absorbing boundary conditions have been
Then, equation (4) is rewritten as:
implemented within the three-dimensional FDTD method for
solving the scattering problem to simulate a microstrip power ~ ~
 H H y  ~
Divider/combiner. It is efficient and can be easily implemented  z
 s y 0 0  D x 
 ~y z  ~ 
j  0 0   Dy 
in the FDTD method. More information and formulation of
~ = sz
fields in uniaxial medium can be found in [6-9]. In this section  H z H z  (7)
  
we describe how to compute the electric and Magnetic fields
 ~z x 
~
 0 0 s x   Dz 
distribution from FDTD method using UPML boundary  H y Hx 
conditions; next we give a technique for extracting scattering  x y 
parameters.  

The time harmonic Maxwell's curl equations in a three- Substitute sx, sy, and sz from equations (6) into (7), and
dimensional UPML medium can be write as [6]: then we use the inverse Fourier transform by the identity
jf ( )  (  ) f (t ) .This gives an equivalent system of

~ 
~ t
  H   j s H . (1)
~ 
~ time-domain differential equations for equation (7):
  E   js H . (2)
~ ~
 H H y 
Where s is the diagonal tensor defined by:  z
 k y 0 0 D x 
 z   0 0   Dy 
 ~ y
~
 H z H z  =  kz
  t
 s x1  ~z x   0 0 k x   Dz 
0 0  s y 0 0  s z 0 0 ~
    H y Hx 
s  0 sx 0  0 s y1 0   0 sz 0   x y 
0  
 0 s x   0 0 s y   0 0 s z1 
 y 0 0   D x 
s y s z s x1 0 
  0  z 0   Dy 
0 1
  (8)
= 0 s x s z s y1 0  (3) 
 0 0  x   Dz 
 0 0 s x s y s z1 


Used equations (1), (2) et (3), Ampere's law in a matched

UPML is expressed as: Yee lattice" is used to discretize equation (8). This leads to
explicit time step expressions for Dx, Dy and Dz.

~ ~
 H H y  For example, the Dx update is given by [6]:
 z
 s y s z s x1 0 0  ~
 ~y z    Ex 
s   H ~
H z   0 s x s z s y1 0  ~ 

z
 E y   2k y   y t   2t 
 0 s x s y s z1 
n 1
 D  
n
z x  0 ~ 
E Dx
 ~ ~   z i  1 , j ,k  2k   t  x
 
i  1 , j ,k  2k   t 
 
 H y
2 2
Hx  y y y y

 x y  n 1 n 1 n 1 n 1
  Hz 2
 Hz 2 Hy 2
 Hy
i  1 , j ,k i  1 , j  1 ,k i  1 , j ,k  1 i  1 , j ,k  1
(4) ( 2 2 2
 2 2 2 2
)
y z
(9)
 After multiplying both sides by par and substituting A similar procedure can be derived to update the
for et from (5a) and (5c), we obtain the equation: components of Ex, Ey, Hx, Hy and Hz using the same approach.

x ~  ~
(k x  ) Dx   (k z  z ) E x (10) B. Computation of S Parameters
j j
For wideband S-parameter extraction, a Gaussian pulse is
Multiplying both sides by j and transforming into the time used as the source excitation. To calculate the parameter S11
domain we can get: for this structure, two computations are required one to capture
the incident waveform at port 1, and a second to capture the
     total wave (incident + reflected) seen at port 1.
( k x Dx )  x Dx    ( k z E x )  z E x  (11)
t   t   The reflected wave form seen at port 1 is determined by
taking the difference at port 1 between the total waveform and
Similarly, from equations (6b) and (6c), then we get a: incident waveform from the structure simulation. The return
loss parameter S11 is calculated by taking the Fourier
 y    transforms of both incident and reflected then computing the
(k y D y )  Dy    (k x E y )  x E y  (12) ratio defined by [12,13]:
t   t  
 FT ( E port
reflected

y  S11  20Log  (17)
1
   incident 
( k z D z )  z Dz    ( k y E z )   FT ( E port1 
 z 
E (13)
t   t
The insertion loss parameter S21 is computed by recording
the value of the Ez field component at port 2. The Fourier
Discretizing the time in equations (11~13) using standard
transforms are then taken for the incident waveform at port 1
Yee. We have , and . For example, the update is
and the transmitted waveform at port 2.
given by:
S21 is given by the ratio:
n 1  2k z   z t   1 
   E x  
n
E x i  1 , j ,k
2
 2k z   z t 
i  1 , j ,k
2
 ( 2k    t )   FT ( E transmited )
z z
S 21  20Log 
port 2
 (18)
incident
2k   t ) D  FT ( E port1 ) 
 (2k x   x t ) Dx in 1 , j ,k 
n 1
 x x x i  1 , j ,k
2 2 
S31 and S41 parameters are computed in the same manner.
(14)

C. Procedure modelling and Simulation

Where x, y, and z namely the space steps in x y z
direction, and t is the time step. A three-dimensional finite-difference time domain FDTD
simulation code with the uniaxial perfectly matched layer
The maximum time step that may be used is by stability UPML absorbing boundary condition [6, 10-16], has been
restriction of the finite difference equations: used to simulate a 3-way power divider/combiner.
The design technique of the proposed structure sets as
1 (15) follows:
t  .
1 1 1
 max    Step 1: Build geometries of proposed structure.
x 2 y 2 z 2
 Step 2: Define parameters.
 Step 3: Initialize all electric and magnetic fields to be
Where  max is the maximum speed of light. zero.
 Step 4: Excite the system by electric source Ez
The input signal is a Gaussian pulse, it has a frequency  Step 5: Updating equations of electric fields.
spectrum from dc to the desired cut-off frequency by adjusting  Step 6: Enforce boundary conditions on the electric
the width of the pulse [13, 14]. The form of the Gaussian fields.
signal is as follows:  Step 7: Update the magnetic fields from equations.
2
 Step 8: Enforce boundary conditions on the magnetic
( t t0 )
E z  f (t )  e 2 (16) fields.
 Step 9: Set the updated electric and magnetic fields for
Where t0 is the pulse delay and  is period of pulse . the next time step, and perform a new cycle by
 repeating between steps 3 to step 6 until the system
comes to a steady state.
 Step 10: Compute scattering matrix coefficients from
time domain results using procedure described in
section B.

D. Design of 3-way Power Divider/Combiner

We will build a structure based on microstrip to design a 3-

way power divider/combiner with wide operating bandwidth,
this structure has a symmetric layout, it's a passive device used
to divide an input signal into three signals, or to combine more
signal sources when used in reverse. In this section, we discuss
the design of circuit, then we present our numerical results to
describe the performance of the proposed design of 3-way
power divider/combiner. Firstly, we give the multi sectional
layout for proposed microstrip circuit in three dimensions (3D)
as shown in figure 1.

The line width, the section length and the total length are Fig. 1. Structure of the proposed power Divider/combiner
indicated with Wj, Li respectively. Index i is starting from 1
and ending with 6 and j from 1 to 2. Then we use a
computational procedure described in section C. E. Optimization

A microstrip is used to design the proposed power divider/ In order to reduce unwanted reflections and improve the
combiner [16], the structure has a simple form, which has four transmission of the proposed power divider/combiner, a
segments, ports 2-4 are three Output ports which receive numerical optimization was successfully used to determine the
power while port 1 is the signal input port, as shown in Figure optimal topology structure for design 3-way power
1. The power divider/combiner has an attractive property, it is divider/combiner. Below, we will describe the proposed
in planar form, and can be manufactured with simple method for optimization utilized in our work which depends
technology without any complicated mechanical assembly, the on the following steps:
interconnections with other planar circuits are simple. This
structure can operate in a large wideband. • Define of power divider/combiner geometry.
• Calculate scattering parameters using 3D FDTD Simulation.
The power divider/combiner is designed with the following Algorithm described in section B.
specifications: RT Duroid substrate 5880 with thickness 2.41,
• Redefine geometry area of topology structure by gradually
width h=0.794, relative permittivity 2.2 and loss tangent
varying its size.
0.0009. The grid cells sizes in x y z are: nx = 80, ny = 100, nz
• Recalculate scattering parameters.
= 16 for 3-way power divider/combiner. The space steps are:
dx = 0.4064e-3, dy = 0.4233e-3 and dz = 0.2650e-3. • This process goes on until we reach to find the best values
[12,15,16,17, 18]. to obtain an overall acceptable performance.
The optimization flowchart can be summarized in Figure 2.
The input port of the circuit is excited by a source Ez and
equally divided at the three output ports, the interior of the
structure is discretized by an FDTD grid [6, 12].

The reflection and transmission coefficient are

computed from the electrical fields at four ports.
 Fig 3a. Ez Field distribution at t = 53 ps in 3-way Power
Divider/Combiner

Fig .2 Flowchart of optimization

A series of simulations have been used to determine the

optimal layout parameters for a minimum insertion loss and
good values of return loss at wideband. The final optimized
dimensions of the proposed 3-way power divider/ combiner
are given in table 1.

Table1.Dimensions of the power divider/combiner

Structure Parameters Values

(mm)
h 0.794
W 2.413 Fig 3b. Ez Field distribution at t = 105 ps in 3-way Power
L1 7.308 Divider/Combiner
Dimensions
L2 5.055
L3 30.856
L4 10
L5 40.650

Three different images of Ez distribution has been captured

at three different times, when the input signal passes through
the device and equally divided between the outputs ports. It
can be observed that the input signal is divided at port 2, 3 and
port 4 as shown in figures 3 (a, b, c).

Fig 3c. Ez Field distribution at t = 221 ps in 3-way Power

Divider/Combiner
 III. CONCLUSION
The scattering transmission parameters for the proposed In this article, a 3-way power divider/combiner was
circuit are shown in Figure (4), the return loss is better than – designed and simulated using the 3D-FDTD Method. The
11 dB over a frequency band from dc to 20 GHz, the results of simulated results indicate that the proposed power
S21, S31 and S41 show performance levels of -4.50 dB from dc divider/combiner is constructed with excellent wideband
to 20 GHz. At Ku band S11 is better than – 14 dB and S21
performance especially at ku band, has simple structure and
lower than -6 dB as shown in figure 4.
easy to process, with good return loss, excellent amplitude and
low insertion loss. The structure is compact, lightweight and
0 low volume, therefore, is suited for low cost applications, they
can also be used as power combiners when used in reverse.
The design of this power divider/combiner can be easily
-5
Scattering parameters (dB)

extended to various multi-band microwave devices using

FDTD method.
-10
.

-15 REFERENCES
[1] Frederick H. Raab, Peter Asbeck, Steve Cripps, RF and Microwave
-20 Power Amplifier and Transmitter Technologies Part 1 May 2003 High
Frequency Design RF Power Amplifier pp. 22, 54..
[2] T. T. Ballen, M.S. Gupta Power combining efficiency and its
-25 optimization IEE Proceeding-H, Vol. 139, N°. 3, June I992.
S11 [3] D. Abouzahra and K. C. Gupta, Multiport power divider/combiner
circuits using circular microstrip disk configurations, IEEE. Microwave
-30 S21
Theory Tech., vol. 35, pp. 1296-1302, Dec 1987.
S31
[4] Yee, K. S., "Numerical solution of initial boundary value problems
S41
involving Maxwell’s equations in isotropic media," IEEE Trans.
-35
0 2 4 6 8 10 12 14 16 18 20 Antennas Propag., Vol. 14, 302-307, 1966.
[5] A. Taflove and M.E. Brodwin Numerical Solution of Steady-State
Frequency (GHz)
Electromagnetic Scattering Problems Using the Time-Dependent
Maxwell's Equations IEEE Transactions on Microwave Theory and
Fig 4. Scattering parameters for 3-way Power Divider /Combiner Techniques, 23(8), 623-630,1975
[6] Allen Taflove, Susan C. Hagness Computational Electrodynamics: The
Finite-Difference Time-Domain Method, 3rd edition. 2005, Pages 661-
662
F. Table performance comparison with some other [7] Gedney, S.D. (1995). An anisotropic perfectly matched layer absorbing
publications medium for the truncation of FDTD lattices. University of Kentucky,
Lexington, KY. Report EMG-95-006.
A comparative study of these results with some recent works [8] Gedney, S.D. (1996). An anisotropic perfectly matched layer absorbing
medium for the truncation of FDTD lattices. IEEE Transactions on
using various techniques is shown in Table 2 for 3-way power
Antennas and Propagations 44, 1630-1639.
divider/combiner. From Table 2, it is clear that the proposed [9] Gedney, S.D. (1998). The perfectly matched layer absorbing medium. In
power divider/combiner with optimum topology has best A. Taflove (Ed.), Advances in computational electrodynamics: The
specifications such return loss and insertion loss. This structure finite-difference time-domain method. Norwood, MA: Artech House.
is well suited for ultra-wide band telecommunication systems. [10] Gedney, S.D., and Lansing, E (1995). Non uniform orthogonal grids. In
A. Taflove (Ed.), Computational electrodynamics: The finite difference
This comparison with other recent publications show a time-domain method. Norwood, MA: Artech House.
significant improvement achieved a -11 dB in Return Loss at [11] Gedney, S.D., and Taflove, A. (2000). Perfectly matched layer
ultra wideband. absorbing boundary conditions. In A. Taflove and S.C. Hagness (Eds.),
Computational electrodynamics: The finite-difference time domain
method. (2nd Ed.). Norwood, MA: Artech House.
T a b l e 2 . Performance comparison with some other articles [12] Jason S. Ayubi-Moak Global Modeling of Microwave Transistors Using
a Full-band Cellular Monte Carlo/full-wave Maxwell Simulator Arizona
State University, 2008, pp. 70- 79.
[13] Khaled ElMahgoub Fan Yang Atef Elsherbeni Scattering Analysis of
Ref Frequency S11 S21 Periodic Structures Using Finite-Difference Time-Domain Method
(GHz) Constantine A. Balanis, Series Editor May 2012.
This work Fig. 4 2 – 20 -11 -4.5 [14] Berenger, J. P., "Perfectly matched layer for the FDTD solution of
wave structure interaction problems ", IEEE Trans Antennas and
[19] Fig. 3 1-7 -2 -5 Propag. Vol 51, 110-117, 1996.
[20] Fig. 3 0.1 - 10 -4.5 -6 [15] Bouayad, M. & Merad, L.. (2018). 3D FDTD design and comparison
performance of power dividers/combiners in different configurations.
[21] Fig 3 a, d 4-6 -6 -5.3 International Journal of Microwave and Optical Technology. 13. 377-
[22] Fig 8 0.5 – 1.5 -7 -5 385.
[16] G. Mur Absorbing Boundary Conditions for the Finite-Difference
[23] Fig 7a 14 - 23 -14 -5 Approximation of the Time-Domain Electromagnetic-Field Equations
This work Fig. 4 12- 18 -14 -6 IEEE Transactions on Electromagnetic Compatibility , Volume: EMC-
23, Nov. 1981.
[17] David M. Pozar Microwave engineering 4th edition John Wiley & Sons,
Inc. 2011 pp. 147
[18] D.M. Sheen, S.M. Ali, M.D. Abouzahra, and J.A. Kong Application of
the three-dimensional finite-difference time-domain method to the
analysis of planar microstrip circuits IEEE Trans Microwave Theory.
Techn MTT-38 (1990), 849–856. 10.
[19] Chun Zhang, Tong Yang Wang, Yong Bo Xiang and Chang Hong
Liang, Wideband three way out-of-phase microstrip power divider
Chang Hong Liang. Volume: 51, Issue: 5 , 3 5 2015, pp. 404 - 406.
[20] The Novel 3-Way Power Dividers/Combiners Structure and
Design 2006 IEEE Annual Wireless and Microwave Technology
conference.
[21] Haijun Fan, Junping Geng, Member A 3-way Reconfigurable Power
Divider/Combiner, IEEE, VOL. 63, NO. 3 March 2015 pp.986-998.
[22] G Yukio Ikeda,1 Fumio Takeda,2 Kazuya Tuda Planar Three-way
Power Divider Using Short Terminated Stubs in the Center of Isolation
Bridges, Electronics and Communications in Japan, Part 2, Vol. 85, No.
8, January 2002, pp. 1–7
[23] Abdulelah A. Alshehri, Saleh M. Alsaif, Abdulrahman, 3-way cascade
power divider and combiner for satellite Communications IEEE June
2017.