Next Article in Journal
Graphene Plasmonics in Sensor Applications: A Review
Previous Article in Journal
Biosensors—Recent Advances and Future Challenges in Electrode Materials
Previous Article in Special Issue
Numerical Optimization of a Fully Cross-Coupled Rectifier Circuit for Wireless Passive Ultra Low Power Sensor Nodes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low-Cost Automated Design of Compact Branch-Line Couplers

by
Adrian Bekasiewicz
Faculty of Electronics, Telecommunications and Informatics, Gdansk University of Technology, 80-233 Gdansk, Poland
Sensors 2020, 20(12), 3562; https://doi.org/10.3390/s20123562
Submission received: 24 April 2020 / Revised: 29 May 2020 / Accepted: 12 June 2020 / Published: 23 June 2020

Abstract

:
Branch-line couplers (BLCs) are important components of wireless communication systems. Conventional BLCs are often characterized by large footprints which make miniaturization an important pre-requisite for their application in modern devices. State-of-the-art approaches to design compact BLCs are largely based on the use of high-permittivity substrates and multi-layer topologies. Alternative methods involve replacement of transmission-line sections of the circuit, with their composite counterparts, referred to as compact cells (CCs). Due to the efficient use of available space, CC-based couplers are often characterized by small footprints. The design of compact BLCs is normally conducted based on engineering experience. The process is laborious and requires many adjustments of topology followed by manual or, semi-automatic tuning of design parameters. In this work, a framework for low-cost automated design of compact BLCs using pre-defined CCs is proposed. The low cost of the presented design technique is ensured using equivalent-circuit models, space mapping correction methods, and trust-region-based local optimization algorithms. The performance of the framework is demonstrated based on three examples, concerning the design of unequal-power split coupler, comparison of automatically generated compact BLCs, as well as rapid re-design of the coupler for different substrates. Furthermore, the approach has been benchmarked against the state-of-the-art methods for low-cost design of circuits.

1. Introduction

Branch-line couplers (BLCs) belong to the key components of many microwave devices including beamforming networks [1,2], balanced mixers [3,4], amplifiers [5,6], multiplexers [7,8], and others. Conventional BLCs are characterized by large dimensions as they consist of two orthogonal pairs of 90° transmission lines (TLs). Miniaturization is an important pre-requisite to facilitate the application of couplers in modern communication systems, such as portable [9,10], or wearable electronics [11,12], but also Internet of Things (IoT) devices [13,14].
Conventional approaches to couplers miniaturization include utilization of high-permittivity substrates [15], or replacement of TLs with their lumped-element equivalents [16]. An interesting alternative is in the replacement of conventional TLs by so-called compact cells (CCs). The latter are normally implemented using high-impedance lines [17], or a combination of high- and low-impedance sections [16], but also as appropriately folded TLs [18], coupled sections featuring unequal length [19], or fractal-based topologies [20]. Due to the reduced phase velocity (also referred to as slow-wave effect) [21,22,23,24], CCs are shorter than conventional transmission lines which is advantageous for design of miniaturized couplers [24]. The state-of-the-art CC-based BLCs offer size-reduction rates ranging from nearly 40% to over 85% [24,25,26,27,28,29]. On the other hand, a small size is often achieved at the expense of degraded electrical properties (most notably in terms of narrower bandwidth or deviation of power split from the desired value) [26,30]. However, couplers featuring small size and high performance have also been reported in the literature [16,31,32]. Furthermore, owing to low-pass properties of CCs, compact BLCs are often capable of suppressing harmonic frequencies [15,24,28,30,31].
Strategies for design of cell-based miniaturized BLCs fall into two categories: (i) development of CCs tailored for specific problem and (ii) selection of the cells from the pre-defined library. The first approach benefits from high miniaturization rates [24,29,33], but it is also laborious and time-consuming [24,29]. The process follows the concept of cognitive design [34]. It involves a manual experience-driven development of cells, but also maintaining their geometrical consistency with adjacent components. The latter is often achieved through definition of complex geometrical relations between individual cells [35,36]. Reuse of pre-defined CCs simplifies BLC development [16,32,37], but at the expense of more modest miniaturization rates. It is worth noting that, for the majority compact of BLCs, miniaturization is achieved as a by-product of performance-oriented tuning rather than dedicated design approach [31,33,38]. Nonetheless, small couplers, obtained as a result of explicit optimization w.r.t. electrical performance and size, have also been reported [39,40].
The selection of cells that ensure high coupler performance while maintaining its small size is a difficult task. The main challenge involves developing a method for unequivocal comparison of CCs [34]. An attempt to address the problem was proposed in [16], where cells are compared in terms of the slow-wave factor [41]. A more generic method concerning pre-optimization of cells with respect to specifications followed by their miniaturization was considered in [24]. Both approaches disregard the importance of optimization for controlling size and performance of the CC-based coupler. Moreover, the computational cost of the methods is high as they exploit only electromagnetic (EM) simulations for evaluation of designs. As indicated in [32,37], this is undesired because performance of the cell can be approximated via simulations of its equivalent-circuit model.
Performance-related tuning of compact BLCs is often a subject to strict requirements concerning, among others, reflection, isolation, phase shift, or bandwidth. Conventional approaches for tuning of structure performance—based on parametric studies followed by visual inspection of responses—are unsuitable for controlling multiple design criteria at a time [30,33]. From this perspective, the application of numerical optimization is a mandatory step of BLC design [42,43,44,45]. Unfortunately, conventional algorithms require a large number of model evaluations to find the optimum solution. At the same time, due to geometrical complexity (particularly in terms of tightly packed cells and multiple cross-couplings between individual components), accurate evaluation of compact BLCs is possible only by means of high-fidelity EM simulations [16,38]. Consequently, the application of conventional optimization routines for the design of miniaturized couplers is a numerically demanding task [35].
Challenges related to high cost of compact BLC design can be addressed using surrogate-based optimization (SBO) [46,47,48,49,50,51,52]. In SBO, the computational burden of numerical optimization is shifted to a surrogate model. The latter consists of a numerically cheap, less accurate low-fidelity model and an appropriate correction layer, which is determined based on occasionally acquired high-fidelity EM simulation data [48,52]. SBO approaches that are suitable for the design of compact microwave circuits exploit both physics-based and functional approximation surrogates [32,47,48]. They include variants of space mapping (SM) [48,53,54], but also feature-based optimization techniques [55,56]. Surrogate-assisted methods tailored to design of miniaturized structures were also considered in the literature [36,37,38,53,54]. A sequential SM technique, where iterative replacement of conventional TL sections with their CC-based equivalents, is followed by surrogate-based refinement of structure response was proposed in [53]. A more advanced method—where correction of the low-fidelity model was first performed at the level of individual cells, and then at the level of the assembled circuit—was introduced in [54]. Another approach, reported in [32], allows for screening of cells based on their performance w.r.t. reference TLs. Then, the selected CCs undergo further optimization using space-mapping-corrected local approximation models [32]. In [55], a feature-based EM-driven optimization of compact couplers was considered. The core of the method is the translation of original frequency characteristics into a set of carefully selected feature points representing structure behavior. Less-nonlinear dependence of response features on changes of input parameters—as compared to frequency characteristics—improves convergence of the optimization and reduces its cost [55,57]. Surrogate-assisted design can be also realized using variable-fidelity EM simulation models. Despite higher-cost than for circuit-model-based methods [32,53,54], variable-fidelity approaches proved to be useful for generation of compact cells [38].
The above considerations indicate that a framework for low-cost design of compact couplers should embed the mechanisms for: (i) screening the pre-defined library in order to select suitable CCs (geometry- and performance-wise), (ii) ensuring geometrical consistency of assembled structure, and (iii) optimizing its performance/size-related properties for a range of specifications. From the design automation standpoint the main challenges are related to determination of suitable starting point for cell-level optimization, but also development of techniques that enable low-cost unequivocal comparison of CCs [32,38,40]. Standardized approaches for ensuring geometrical consistency of compact couplers, although available, support only structures that combine CCs and meandered TLs [32,37,38]. For other components, manual definition of constraints and mutual relations between individual cells is required [24,35,54]. Finally, state-of-the-art techniques for low-cost design of couplers are either difficult to setup and unsuitable for automation [32,53], or too expensive to justify their use as a tool for, e.g., the optimization-driven comparison of topologies [35,38]. The implication is that existing methods are dedicated for solving specific narrow classes of problems pertinent to design of compact BLCs. In this regard, a more holistic approach oriented for reliability and simplification of the design process is required.
In this work, a framework for unsupervised design of miniaturized branch-line couplers is presented. The method follows the concept of bottom-up design, where the BLC development is realized in the following steps: (i) selection/optimization of compact cells from pre-defined library, (ii) BLC assembly followed by its two-stage optimization, and (iii) surrogate-assisted fine tuning of the coupler. The proposed method exploits circuit and EM simulation models, low-fidelity model correction schemes, and feature-based representation of the structure responses. The response features are used at the cell- and BLC-model levels. In each design stage, the optimization is performed using a gradient algorithm embedded in a trust-region framework. The introduced methodology is validated using three case studies concerning design of miniaturized coupler with unequal power split, comparison of BLC topologies in terms of achievable miniaturization rates, and rapid re-design of coupler for different substrate materials. The tool has been benchmarked against the other numerically efficient design methods in terms of computational cost and quality of obtained solutions. The average cost of BLC development using the presented methodology amounts to only 17.5 min. Moreover, it is capable of generating topologies featuring miniaturization rates of almost 80 percent.
The remainder of the paper is organized as follows. Section 2 discusses the models and algorithms used in the work. Implementation of the developed framework is explained in Section 3, whereas Section 4 contains numerical verification of the proposed methodology and discussion of the results. Section 5 concludes the paper and indicates the directions for future research.

2. Simulation Models and Design Tools

The section describes physics-based simulation models employed in the proposed framework. In order to make the paper self-consistent, a brief explanation of the core methods and algorithms used in the work is also included here. The details concerning implementation of the proposed framework for automated design of couplers are provided in Section 3.

2.1. Transmission Line-Based Sections vs. Compact Cells

Microstrip transmission lines belong to basic building blocks of many conventional microwave circuits, including BLCs. Due to inefficient use of available space, TL-based circuits are often characterized by unacceptably large footprints [24,25]. The problem can be addressed by replacing TLs with compact cells. CCs mimic electrical behavior of transmission lines while offering reduced physical length [16,31]. The cells also offer increased number of degrees of freedom compared to TLs (controlled only through adjustment of electrical length and characteristic impedance [58,59]), which provides additional flexibility in the course of circuit design process [29,38]. As already mentioned, compact BLC can be constructed using off-the-shelf cells (easy to reuse but offer moderate miniaturization rates), or structures developed manually based on the principles of cognitive design (laborious and time consuming but desirable from miniaturization standpoint) [34]. In this work, a balance between simplicity of the design process and achievable size reduction rates is maintained using asymmetrical CCs in appropriate configuration. Specifically, the BLCs considered here combine cells in vertical and horizontal orientations. Such a complementary setup is convenient, as it maintains the geometric consistency of coupler topology using a single parameter that controls separation between the cells. Another benefit is that the cells can be easily interchanged which facilitates rapid generation of various geometries. A conceptual comparison of conventional and miniaturized BLC structures is shown in Figure 1.

2.2. Compact Cell and BLC Models

Let Cf.o(xo,xgap,so) be the response (S-parameters versus frequency) of a generic model of the compact cell. Here, xo represents the vector of model input parameters, whereas xgap is a global, fixed parameter used to maintain desired separation between adjacent CCs that constitute the assembled BLC structure. The vector so = [Zo h εr t tanδ]T represents setup parameters, such as impedance and substrate properties (i.e., height, permittivity, metallization thickness, and loss tangent). The subscript parameters f ∈ {l, h, s} and o ∈ {H, V} represent fidelity of the model and configuration of the cell, respectively. In other words, replacement of f with l, h, or s means that the structure is evaluated using either the low-fidelity equivalent circuit model, the high-fidelity EM model, or surrogate model. Similarly, models and parameter vectors with H or V in subscript indicate cells in horizontal and vertical configurations (cf. Figure 1), respectively. For simplicity of notation, let Cf(x) = Cf.o(x) = Cf.o(xo,xgap,so). External dimensions of individual CC are defined as [ao.w ao.h xgap]T = AC(Cf(xo)), where ao.w and ao.h stands for width and height of the cell, whereas AC is the function that calculates the size. Conceptual illustration of the cell in horizontal configuration is shown in Figure 2.
The BLC response is calculated based on the transmission line theory from complex electrical characteristics of individual cells [58,59,60]. A MATLAB-based implementation of the considered coupler model supports seamless and automated switching between the CCs pre-defined in the library, as well as unsupervised evaluation of the resulting topologies. The coupler responses are calculated sequentially between selected pairs of ab ports (see Figure 1 for ports enumeration) as,
B a b ( y ) = V a b ( C f . H ( x H , x g a p , s H ) , C f . V ( x V , x g a p , s V ) )
where y = [xH xV]T. The function Vab realizes the electrical connection that is specific for the selected ab pair [59]. The response of (1) is of the form Bpair.ab(y) = [Baa(y) Bab(y) Bba(y) Bbb(y)]. Note that yielding sufficient information on BLC performance requires evaluation of only ab ∈ {12, 13, 14} pairs, i.e., Bpair.12(y) = [B11(y) B12(y) B21(y) B22(y)], Bpair.13(y) = [B11(y) B13(y) B31(y) B33(y)], and Bpair.14(y) = [B11(y) B14(y) B41(y) B44(y)]. Then the selected responses are aggregated into the following matrix:
B ( y ) = [ B 11 ( y ) B 12 ( y ) B 13 ( y ) B 14 ( y ) ] .
It should be stressed out that (1)–(2) support the high-fidelity, low-fidelity and surrogate-based representations of cells, i.e., B(y) = Bf(y), f ∈ {h, l, s}. The benefits of the considered model are as follows (i) it operates only on electrical responses and thus supports seamless combination of cells developed using various software packages, (ii) it is universal as the type of the circuit under design (here, BLC) can be changed only through modification of Vab function, (iii) it can be rapidly re-used to obtain structure responses for various combinations of cells topologies. An important remark is that (2) neglects the cross-coupling effects between the cells which limits its accuracy [35,37]. These, however, are accounted for in a high-fidelity EM model of the assembled coupler (cf. Figure 1), denoted as BA(y). The latter is used only at the final stage of the design process. The detailed block diagram that summarizes calculation of BLC responses based on characteristics of individual CCs is shown in Figure 3.
The BLC footprint can be calculated as follows,
A B ( y ) = max { a H . w , 2 a V . w + a H . g } ( 2 a H . h + a V . h + 2 w 0 )
where w0 represents the width of coupler input ports (calculated based on transmission line theory for the given substrate parameters) [58,59]. Conceptual illustration of the considered universal coupler model is shown in Figure 4.

2.3. Problem Formulation

The problem concerning design of the individual cell or the coupler can be formulated as the following non-linear minimization task [62],
z = arg min z U ( R ( z ) )
where R(z) is response of the structure under design obtained for the given vector of design parameters z. Here, U represents a scalar objective function. It should be noted that for the cell-level design z = xo and R(z) = Ch.o(x). Similarly, z = y and R(z) = B(y) when coupler-level design is considered. The goal of (4) is to find the optimal design z* through minimization of the objective function.

2.4. Surrogate-Assisted Optimization

Direct solving of (4) is numerically expensive when R(z) is evaluated using high-fidelity representation of the structure at hand. Computational cost of the design process can be substantially reduced using surrogate-based optimization. SBO generates a series of approximations, j = 1, 2, …, to (4) by solving [62],
z ( j + 1 ) = arg min z U ( R s ( j ) ( z ) )
Here, Rs(j)(z) represents the surrogate model response of the structure under design at jth iteration of (5). The surrogate is implemented as the equivalent-circuit model enhanced using the correction layer. The latter is determined based on occasional evaluations of the high-fidelity model which are performed to verify the quality of z(j+1) designs [48,62]. The high-fidelity data is also used to gradually improve the Rs(j)(z) accuracy by updating the correction layer. In this work, the low-fidelity model correction is applied only at the level of individual cells. The correction layer is implemented as a combination of the implicit space mapping (SM) and frequency scaling [62,63]. The surrogate model of the cell Rs(j)(z) = Cs(j)(x) is of the following form,
C s ( j ) ( x ) = C s ( j ) ( x , p ( j ) ) = C c ( x , x c ( j ) , α ( j ) ω )
where p(j) = [xc(j) α(j)]T is a set of parameters used for adjusting the correction layer. The vector xc(j) = [h1 h2εr1 εr2]T represents a set of control variables (also referred to as the preassigned parameters) for implicit space mapping (ISM), i.e., substrate thickness h and permittivity εr. They can be adjusted to reduce discrepancy between the surrogate and the high-fidelity model responses [63,64]. The component α(j)ω = α0(j) + α1(j)ω represents frequency scaling. It is used to “shift” and/or “stretch/squeeze” response of the surrogate w.r.t. the original frequency sweep ω.
In each SBO step, concurrent update of the correction layer parameters is performed. The adjustment is realized in the course of a so-called parameter extraction (PE) [48,64]—a curve fitting process oriented towards matching the low-fidelity model response Cl to the high-fidelity model Ch. Here, PE is defined as the following minimization task:
p ( j ) = arg min p C h ( x ( j ) ) C l ( x ( j ) , p ( j ) ) .
A conceptual illustration of the CC surrogate with highlight on the low-fidelity model and the correction layer is shown in Figure 5. For more detailed discussion on SBO, as well as low-fidelity model correction methods, see [49,51,54,62,63].

2.5. Feature-Based Representation of Structure Responses

One of the main challenges related to development of reliable methods for unsupervised design is that responses of couplers—and the responses of their building blocks—are highly nonlinear functions of frequency [55]. To be considered generic, CCs stored in the library have to support a broad range of substrate and geometry parameters, as well as operate within relatively wide frequency spectrum. Additionally, automated selection of cells should be performed as a result of rigorous optimization process. However, development of objective function that accurately represents cell performance for a wide range of design scenarios and input parameters is difficult. Yet another issue associated with optimization-based cell-selection is the risk that algorithm will get stuck in poor local optimum. The difficulties concerning performance adjustment become even more pronounced when cells are used for construction of structures with complex response characteristics, such as BLCs.
Mentioned challenges can be addressed by converting frequency responses of the structure into an appropriate domain. Here, this task is realized by reformulating the frequency characteristics to the form of response features [57]. The latter represent each key feature of the structure frequency response using a pair of carefully selected coordinates. Feature points can be defined as follows. Let F(z) = P(R(z)) be the response of the model R expressed in the form of feature points, where P is the function that realizes the transformation. The feature-based response of the structure at hand is given by,
F ( z ) = [ ω ( R ( z ) ) l ( R ( z ) ) ] T = [ ω l ] T = [ ω 1 ω N l 1 l N ] T
where ωn and ln, n = 1, …, N represent the frequency point and its corresponding value of the response (also referred to as level), respectively. As shown in Figure 6, the nth coordinate can be defined either for a fixed frequency point, or for a desired level of circuit response. This flexibility makes features a convenient tool for tracking the key changes of structure responses. Another advantage is less non-linear dependency of the features as a function of input parameters compared to original frequency characteristics (cf. Figure 7 and Figure 8). The mentioned properties make this form of data representation an efficient tool for addressing challenges related to unsupervised design of compact couplers. In particular, description of circuit responses in terms of response features improves convergence of the optimization process and reduces its computational cost, as compared to design tasks that involve direct evaluation of frequency responses [56,57].

2.6. Optimization Engine

The optimization engine used in this work for solving (5) and (7) is a gradient algorithm embedded within a trust-region (TR) framework. It takes the form of [44],
z ( j + 1 ) = arg min z : z z ( j ) r ( j ) U ( G ( j ) ( z ) )
where G(j)(z) is the first-order Taylor expansion model constructed at the level of response features,
G ( j ) ( z ) = F f ( z ( j ) ) + J ( z ( j ) ) ( z z ( j ) )
Here, Ff = P(Rf) represents the feature-based responses of the structure at hand obtained at either the low- (f = l), high-fidelity (f = h), or surrogate model level (f = s), respectively. The Jacobian J(z(j)) is obtained using a large-step finite differentiation as [44]
J ( z ( j ) ) = [ d 1 1 ( F J . f ( z ( j ) + e 1 ) F J . f ( z ( j ) ) ) d D 1 ( F J . f ( z ( j ) + e D ) F J . f ( z ( j ) ) ) ] T
The responses FJ.f = P(RJ.f) used for construction of the Jacobian are obtained from the low-fidelity RJ.f = Rl or surrogate RJ.f = Rs model simulations. The vector ek = [0 … dk … 0]T, k = 1, …, D, is the perturbation size w.r.t. kth design parameter (D corresponds to the length of the z vector, i.e., the number of geometry parameters that represent the structure at hand) and d = [d1dkdD]T is the vector of perturbations obtained for the structure under design. Note that, depending on the design stage, the vector of response features Ff can be obtained from the responses of either Cf, or BA models (cf. Section 2.2). Perturbations for construction of the low-fidelity-model-based Jacobian are set to 10−3.
In each iteration of (9), the gain ratio ρ—which represents actual change of the objective function obtained for the high-fidelity model versus the one predicted by the Taylor-expansion model—is calculated as [44]:
ρ = U ( F f ( z ( j + 1 ) ) ) U ( F f ( z ( j ) ) ) U ( G ( j ) ( z ( j + 1 ) ) ) U ( G ( j ) ( z ( j ) ) ) .
The gain ratio is used to assess quality of the solution obtained from (9). When ρ > 0, the design z(j+1) is accepted and used as a starting point for the next iteration. Otherwise, it is rejected and z(j) is re-used for optimization. The coefficient ρ is also used to update the TR radius r(j) as follows,
r ( j + 1 ) = { max ( η 1 z ( j + 1 ) z ( j ) , r ( j ) ) ,   ρ > 0.9 η 2 z ( j + 1 ) z ( j ) ,   ρ < 0.3 r ( j ) ,   otherwise
where the scaling factors for update of the TR radius are set to η1 = 2.5 and η2 = 0.25 [44]. The default initial radius is r(0) = 1.
The algorithm is terminated when the obtained objective function value U(Ff(z(j+1)) is below 0 or when either of the following conditions is fulfilled: (i) the Euclidean distance between consecutive optimal designs, (ii) TR radius, or (iii) average change of objective function value for three consecutive iterations:
z ( j + 1 ) z ( j ) ε
r ( j + 1 ) ε
i = j 1 j + 1 | U ( F f ( z ( i 1 ) ) ) U ( F f ( z ( i ) ) ) | ε .
Here, the user-defined threshold value is set to ε = 10−2. It should be noted that the numerical cost of (9) is only D + 1 evaluations of the structure per successful iteration. Additional simulations are required for ρ < 0.
The optimization algorithm can be summarized as follows:
(1)
Set j = 0, z(j) = z0, r(j) = 1;
(2)
Evaluate Ff(z(j)) and select model Fl = P(Rl) for construction of Jacobian;
(3)
(Optional) Perform cell-level PE at z(j) as described in Section 2.5; select model Fs = P(Rs) for Jacobian construction;
(4)
(Optional) If j = 0, estimate r(j) as explained in Section 3.3;
(5)
Generate perturbations around z(j), construct the Jacobian and the G(j) model;
(6)
Solve (9) to obtain a temporary solution ztmp;
(7)
Evaluate Ff(ztmp), calculate ρ as in (12) and adjust radius r(j+1) as in (13);
(8)
If ρ > 0 set z(j+1) = ztmp; otherwise set z(j+1) = z(j), j = j + 1 and go to Step 6;
(9)
If U(Ff(z(j+1)) < 0 or either of conditions (14)–(16) is satisfied then END; otherwise set j = j + 1 and go to Step 2.

3. Automated Surrogate-Assisted Design of Compact Branch-Line Couplers

This section describes the proposed framework for automated surrogate-assisted design of couplers. In particular it provides a detailed discussion of each design step, i.e., cell-level optimization, two-stage design of low-fidelity BLCs, and surrogate-assisted tuning of high-fidelity coupler model. Each of the mentioned steps build up on the methods and algorithms of Section 2. The last part of the section contains summary of the framework. An important remark is that the main prerequisite for application of the presented framework to unsupervised design of couplers is availability of a feasible user-defined specifications (e.g., center frequency, bandwidth, power-split ratio, etc.) that are required to determine electrical parameters of individual cells used for coupler construction [32], as well as for controlling the objective functions. It should be reiterated that, to ensure low computational cost, each step of the design process that involves optimization is performed using a gradient algorithm embedded in the TR framework.

3.1. Optimization of the Pre-Defined Cells

As already mentioned, the conventional BLC comprises two pairs of horizontal and vertical TL sections with electrical length θ0 and characteristic impedances of ZH, and ZV, respectively. TLs parameters, that are appropriate for obtaining the desired coupler performance, can be determined based on the transmission line theory [58,65,66].
In the first design step, optimization of pre-defined CCs is performed. The process is oriented towards searching for the electrical properties of TLs that would normally be used for construction of conventional coupler. The feature-based objective function is given as:
U C ( F f ( x ) ) = max ( l 3 l 3 max , 0 ) + max ( | l 4 θ 0 | l 4 max , 0 ) + ω 1 ω 0 ω 2 ω 0 .
As shown in Figure 6, ω1 and ω2 represent frequency points at which reflection of the cell has its minimal value and where the phase shift is equal to the desired value θ0, respectively. The coordinates l3 and l4 represent the reflection and phase shift at the center frequency ω0. The user-defined parameters l3max = 0.003 and l4max = 0.1° are used to activate/deactivate the first two components of (17). Note that cells optimization is performed at the low-fidelity model level (i.e., f = l; cf. Section 2.2). The feature-based response is calculated individually for each cell, i.e., Ff(x) = Ff.V(x) = P(Cf.V(xo)), or Ff(x) = Ff.H(x) = P(Cf.H(xo)).
The starting point for cell-level optimization is selected from a set of random samples generated as follows:
X = { x | x = ( u b l b ) y + l b , y R D , 0 y 1 } .
Here, lb and ub represent the lower and upper bounds on geometry parameters of the cell at hand. The symbol “○” represents component-wise multiplication. The starting point for cell optimization is obtained from:
x init = arg min x X ( U C ( F l ( x ) ) ) .
In other words, xinit is the design from the set X for which the objective function (17) has the lowest value. The cell is optimized by solving (9). Due to utilization of equivalent-circuit model representation, the cost of finding the initial and optimal solutions is negligible. Finally, the parameters of the optimized cells xH* and xV* (cf. Section 2.2) are concatenated and used as a starting point for coupler design.
It should be emphasized that the role of cell-level optimization is twofold. On one hand, the procedure rejects the cells that are incapable of fulfilling the user-defined specifications. On the other hand, it allows to find optimal design parameters of selected CCs.

3.2. Two-Stage BLC Optimization

Upon completion of selection and optimization procedure, the CCs are used to construct the BLC model (cf. Section 2.2) which then undergoes a two-stage design process. In the first step, the coupler is optimized w.r.t. performance requirements. The design criteria include: (i) maximization of bandwidth, (ii) maintaining minimum of S11 and S41 at ω0, and (iii) minimization of power-split error at ω0. The objective function is defined as,
U B .1 ( F ( y ) ) = β 1 ( B W + B 0 ) + ω 5 ω 0 ω 6 ω 0 ω 7 ω 0 + ( max ( | Δ C Δ C 0 | Δ C max , 0 ) Δ C max ) 2 + + β 2 ( max ( { l 5 , l 6 } S max , 0 ) S max ) 2
where F(y) = P(Bf(y)). The feature coordinates ωn and ln, n = 1, 2, …, 9 are defined as shown in Figure 6. The parameters ωl = max(ω1, ω2) and ωh = min(ω3, ω4) represent lower and upper corner frequency for which |S11| and |S41| are both below the specified threshold (here, −20 dB). The bandwidth is defined as BW = 2min{ωhω0, ω0ωl,}. Parameters ΔC0, ΔCmax, and Smax denote desired power-split imbalance at the center frequency, acceptable power-split error—i.e., deviation of ΔC = |l8l9| from the desired value of ΔC0 (here, 0.2 dB)—and maximum allowed reflection level of S11 and S41 at ω0 (here, −26 dB). The scaling coefficients are set to β1 = 10 and β2 = 1000. The particular values are determined based on the numerical experiments.
The starting point for performance-wise coupler optimization is y0 = [xH* xV*]T (cf. Section 3.1). The final design is found through minimization of (20) by solving (9). The obtained solution y1* is set as the starting point for a second-stage of the low-fidelity-model-based BLC design, i.e., minimization of coupler size. The objective function is given by:
U B .2 ( F ( y ) ) = A B ( y ) + β ( max ( U B .1 , 0 ) ) 2 .
Here, β = 10. Due to representation of the BLC using equivalent-circuit model, the cost of minimizing (20)–(21) is low. The final design (here, denoted as y2*) is used as a starting point for surrogate-assisted optimization of the assembled coupler, as described below.

3.3. Surrogate-Assisted Optimization

The final design stage involves surrogate-assisted optimization of the assembled coupler. In each SBO iteration, the feature-based linear model (10) is constructed from a single EM response of the BA model that represents the assembled BLC. The Jacobian J is calculated based on the responses of the BLC surrogate model (f = s). In other words, in (10), Ff(z(i)) is calculated based on the BA response obtained at design z(i), whereas J(z(i)) is constructed from Bs simulations—determined based on corrected responses of compact cells—obtained around the z(i) design. The surrogate models of horizontal and vertical cells that constitute the Bs model are individually corrected to the high-fidelity model level using implicit SM and frequency scaling (cf. Section 2.4). The rationale behind correction of individual cells is that PE can be carried out using relatively small number of control parameters, which is beneficial for fast convergence of TR-based optimization. Furthermore, the BLC surrogate, composed of individually corrected CCs, is characterized by good generalization capability [54].
At the beginning of SBO, the initial TR radius is estimated. The proposed estimation process is realized in a separate TR loop that exploits linear approximation model G constructed from Bs responses. The algorithm can be summarized as follows:
(1)
Set j = 0, rtmp(0) = 0.5;
(2)
Set y(j) = y2*;
(3)
Find y(j+1) through minimization of the objective function (20) by solving (9);
(4)
Evaluate Fs(y(j+1)) = P(Bs(y(j+1))) and calculate ρ;
(5)
If ρ > 0.5, set rtmp(j+1) = 0.5(rtmp(j) + 1); otherwise set rtmp(j+1) = 0.5(rtmp(j) + 0.1);
(6)
If j = 5 set r(0) = rtmp(j+1) and exit; otherwise set j = j + 1 and go to Step 2.
It should be reiterated that the surrogate model Bs does not account for cross-couplings effects between the cells and therefore exhibits a certain level of inaccuracy (most notably in the form of frequency shift, or power-split discrepancy w.r.t. the BA model). This is accounted for through refinement of the estimated radius as r(0) = r(0)sF, where sF is the scaling factor that averages the discrepancy between ω5, ω6 features obtained for the surrogate and fine model. The optimization process is oriented towards minimization of (20) using the version of the TR algorithm that employs both optional steps (cf. Section 2.6). Moreover, additional constraint is introduced to (9) in order to ensure that the footprint of optimized BLC does not exceed AB(y2*) (cf. Section 2.2 and Section 3.2). The vector y* obtained after completing SBO is the final solution of unsupervised design process.

3.4. Summary of the Design Framework

The proposed design framework can be summarized as follows:
(1)
Define coupler-level specifications and determine the electrical parameters of individual BLC sections that correspond to the design requirements;
(2)
Find the starting point for cell-level design as described in Section 3.1;
(3)
Minimize (17) to select and optimize CCs that are capable of fulfilling the imposed specifications;
(4)
Use the selected CCs for construction of the Bc model. Sequentially minimize (20) and (21) to obtain y1* and y2*;
(5)
Perform SBO using BA and Bs models to obtain the final BLC design y* (cf. Section 3.3).
The user-defined figures for automated design include electrical length of cells, center frequency, power-split imbalance, bandwidth, and substrate properties. It should be noted that characteristic impedances of CCs required to obtain selected ΔC0 can be calculated from [65],
Z H = Z 0 C 0 / ( 1 + C 0 ) ,   Z V = Z 0 C 0
where C0 = 10ΔC0/10 and Z0 = 50 Ω. Note that for power-split imbalance ΔC0 = 0, the BLC features equal (3-dB) power-split. It should be reiterated that, owing to embedding each design step (including PE of cells) into the TR loop, the overall cost of coupler design using the proposed framework is low. Successful iterations of surrogate-assisted BLC refinement require only one EM simulation per CC (for PE) and one EM evaluation of the assembled coupler. Unsuccessful steps require additional simulations of BA (PE is performed only for successful iterations).

4. Numerical Results

The main focus of this section is description of the library of pre-defined cells followed, as well as numerical validation of the proposed design framework. The methodology is demonstrated based on three examples concerning (i) design of a compact coupler with unequal power-split, (ii) comparison of BLCs characterized by different topologies w.r.t. size/performance trade-off, and (iii) rapid re-design of structures for various substrate parameters. The section is summarized by the discussion and comparison of the proposed framework against the methods from the literature. It should be noted that all numerical experiments have been performed on a dual Intel Xeon E5540 machine with 32 GB of RAM. Moreover, unless stated otherwise, the dimensions of all structures considered in this section are expressed in mm.

4.1. Database of Compact Cells

The topologies of the CCs used in this work are shown in Figure 9. The library contains four horizontal CH(i), i = 1, 2, 3, 4, and four vertical CV(j), j = 1, 2, 3, 4, cells. Each is implemented at the level of both equivalent-circuit and EM simulation models. As already mentioned in Section 2, the cells are described using a set of geometry variables xo = [xo.g1 xo.g2 …]T, setup parameters so, and dimension xgap that represents distance between cells embedded into the BLC circuit. When necessary, the cells also implement the correction parameters p = [xc(j) α(j)]T. The meaning of all mentioned variables as well as their role within the cell model is illustrated in Figure 5. Apart from containing information on cell orientation, the database encompasses ranges of design parameters for which CC geometries remain feasible, functions for calculating cells dimensions, etc. For simplicity of notation, the coupler composed of ith horizontal and jth vertical cell is denoted as B(i,j). The high-fidelity EM models and low-fidelity equivalent-circuit models of cells are implemented in Keysight ADS [67].

4.2. Design of Unequal-Split BLC

The first example is the design of an unequal-power-split coupler implemented on a Rogers RO4003C substrate (εr = 3.38, h = 0.813 mm, tanδ = 0.0021). The following design specifications are considered: center frequency ω0 = 1.5 GHz, power-split imbalance ΔC0 = 3 dB, and bandwidth BW ≥ 150 MHz. The coupler is composed of cells CH(4) and CV(3) (cf. Figure 9). Their characteristic impedances—calculated from (22)—are ZH = 40.8 Ω and ZV = 70.6 Ω, respectively. The gap is set to xgap = 0.3 mm.
The structure was designed according to the methodology outlined in Section 3.4. In the first step, the starting points for CCs optimization were obtained using procedure of Section 3.1 and the designs xH* = [2.49 0.41 2.29 0.49 0.66 0.48]T and xV* = [3.38 0.2 0.2 0.82 1 1]T were found through minimization of (17). Comparison of cells responses at the initial and optimized designs is shown in Figure 10. Next, a two-stage optimization of the coupler was performed as described in Section 3.2. The initial design was set to y0 = [xH* xV*]T. The design y1* = [2.6 0.37 1.44 2.16 0.9 0.28 3.71 0.2 0.28 0.39 1 0.98]T, obtained by minimization of (20) using (9), was set as a starting point for size-oriented coupler optimization. The final design y2* = [2.69 0.33 1.45 2.16 0.9 0.28 3.64 0.2 0.25 0.36 0.99 0.98]T was found through minimization of (21). The coupler footprints at y0, y1* and y2* are 295 mm2, 285 mm2, and 276 mm2, respectively. Low-fidelity model responses of the structure at y0 and y2* designs are compared in Figure 11a. It should be noted that the geometry obtained directly from optimization of individual CCs already offers acceptable characteristics. The two stage optimization further improved the performance and reduced the BLC footprint. In the last step, surrogate-based optimization of the assembled coupler BA that accounts for cross-couplings between individual cells—was performed as described in Section 3.3. The high-fidelity optimized design y* = [1.97 0.29 1.68 2.31 0.84 0.39 3.98 0.36 0.34 0.2 0.78 0.9]T was found after only four SBO iterations. Responses of the high-fidelity EM model at y2* and y* designs are compared in Figure 11b. The dimensions of the optimized structure are 10 mm × 28 mm = 280 mm2, which corresponds to 75.5% size reduction compared to conventional BLC coupler (dimensions: 31.8 mm × 36 mm = 1144.8 mm2).
Geometry of the optimized structure, as well as comparison of the BLC responses obtained from EM simulations performed using Keysight Momentum and CST Microwave Studio (time-domain solver) [68] are shown in Figure 12. The results are in good agreement. It should be reiterated that the coupler design was performed without user supervision. The overall computational cost of BLC design corresponds to just 12.1 evaluations of the BA model (around 11 min of CPU-time) including 4 BA simulations, a total of 6 EM simulations of CH.h(4) and CV.h(3) required for SBO refinement, as well as some computational overhead for PE and TR-based optimization.
The proposed design approach was benchmarked against the state-of-the-art SBO algorithms [44,48,62]. The considered methods include: (i) direct TR-based optimization of BA model, (ii) two versions of implicit SM executed at the level of Bl model, (iii) two versions of ISM combined with frequency scaling, and (iv) the version of algorithm of Section 3 that does not exploit the mechanism for automated estimation of the initial radius at the beginning of SBO. The figures considered for comparison include optimization cost, as well as bandwidth, power-split, phase-shift, and size of the final designs. For fair comparison, the design y2* was set as a starting point for each benchmark algorithm. The results shown in Table 1 indicate that ISM-based routines failed to obtain satisfactory solutions due to divergence. Algorithm (i) yields the best results in terms of structure size and bandwidth. The designs obtained using (iii) are characterized by the largest footprints and highest deviations of ΔC from the target value. The proposed design approach yields relatively small design while ensuring decent electrical performance. It should be noted that the bandwidth of designs obtained using benchmark algorithms is broader compared to the one determined using the proposed algorithm. However, they also feature worsened power-split imbalance. The reason is that the objective function used for BLC optimization aggregates performance-related design objectives into a weighted sum [45,69]. From the objective function standpoint, bandwidth enhancement compensates the degradation of power-split imbalance.
The proposed design optimization approach provides the lowest computational cost among algorithms that yield acceptable designs. For the considered example, the routine of Section 3 is about 85% faster compared to (i), but also its cost is up to 82% lower compared to methods based on algorithm (iii). Moreover, the proposed algorithm is over 32% faster compared to the method that does not involve automatic determination of the initial TR radius. It should be emphasized that significant computational savings w.r.t. algorithms (i)–(iii) result from embedding the parameter extraction process into the TR framework, which dramatically reduces the number of required simulations. As can be seen from Table 1, this change reduces the number of low-fidelity model evaluations required for SBO by one to two orders of magnitude. Another advantage of the proposed algorithm is that parameter extraction is performed at the level of individual cells of the coupler, which improves generalization of the BLC surrogate compared to (i)–(iii).

4.3. Size/Performance Comparison of Compact BLCs

The second example concerns comparison of compact BLCs with respect to their performance and size. A total of 16 topologies, representing all combinations of cells from the library of Section 4.1, are used for analysis. The design specifications are as follows. The couplers are to be implemented on a Rogers RO4003C substrate (εr = 3.38, h = 0.813 mm, tanδ = 0.0021).
The performance-related design objectives include center frequency of 1 GHz, equal-power split, and at least 90 MHz bandwidth. Based on (22), the equal-split (ΔC0 = 0 dB) is obtained for ZV = 50 Ω and ZH = 35.4 Ω, respectively. The gap between the adjacent cells is fixed to xgap = 0.3 mm. First, all CCs from the database were optimized as described in Section 3.1. Due to insufficient performance, the cell CV(3) was rejected by the algorithm from the design process. Next, a set of 12 coupler topologies composed of CH(j), i ∈ {1, 2, 3, 4}, horizontal and CV(j), j ∈ {1, 2, 4}, vertical cells that feature satisfactory performance was selected for two-stage optimization. Finally, the resulting BLC designs were refined to the high-fidelity model level using surrogate-assisted optimization.
Table 2 shows comparison of the couplers in terms of bandwidth, power split imbalance, phase shift, size, and miniaturization. To provide more unequivocal measure of structures properties, a performance-to-size coefficient (PtS) was defined as,
P t S ( i , j ) = ( B W ( i , j ) 200 max ( Δ C ( i , j ) Δ C max Δ C max , 0 ) 2 ) 1 0.1 A B ( i , j )
where BW(i,j), ΔC(i,j) and AB(i,j) denote bandwidth, power-split imbalance, and size of the optimized BA(i,j) coupler. The parameter ΔCmax = 0.2 dB (cf. Section 3.2). Owing to being proportional to quality of the BLC with respect to design requirements, PtS can be used to sort BLCs w.r.t. their performance. Structures with negative PtS are considered unacceptable from the specification standpoint. The results of Table 2 indicate that only five of twelve couplers, i.e., BA(1,2), BA(1,4), BA(1,1), BA(2,1), and BA(3,1), feature acceptable performance. The remaining ones are characterized by either too narrow bandwidth or unacceptably high power-split imbalance. It should be noted that even the designs with PtS > 0 may slightly violate requirement concerning ΔCmax. As already mentioned, this is due to aggregation of objectives in (20) using weighted sum. Therefore, slight violation of the requirements can be considered acceptable. Frequency characteristics and geometries of couplers with satisfactory performance are shown in Figure 13, whereas their detailed dimensions are gathered in Table 3.
Another important remark is that the performance of the optimized BLCs vary significantly with topology. In other words, compact dimensions of a coupler often come at the expense of performance deterioration [26,70]. A satisfactory trade-off between the size and the performance of the BLC can be achieved, provided that appropriate combination of cells is selected. The results are on par with findings from the literature references, indicating that size- and performance-related objectives are conflicting [16,37,38].
The average design cost is only 18.2 min of CPU-time per coupler. The cost accounts for two-stage optimization of the assembled structure low-fidelity model (~3 min per design), and surrogate-assisted refinement to the EM model level. The overall design cost of 12 BLCs amounts to 3.4 hours. It should be reiterated that all the couplers were designed without external supervision.

4.4. BLC Re-Design for Substrates with Different Parameters

The last example concerns re-design of the coupler for substrates with different parameters. The structure of choice consists of CH(2) and CV(4) cells. The gap between CCs is set to xgap = 0.2 mm. Geometry of the coupler is shown in Figure 14. The design specifications are as follows: ω0 = 0.8 GHz, ΔC0 = 0 dB and BW ≥ 80 MHz. The structure is re-designed for operation on Arlon AD250C (εr = 2.5, h = 0.762 mm, tanδ = 0.0014), Arlon AD300C (εr = 2.97, h = 0.762 mm, tanδ = 0.002), and FR-4 (εr = 4.3, h = 1 mm, tanδ = 0.02) substrates, respectively.
For each substrate the coupler was designed using the framework of Section 3. Dimensions of the optimized structures and their performance figures are gathered in Table 4 and Table 5, respectively. Figure 14 shows frequency responses of the obtained designs. Differences in dimensions of the individual designs indicate that the couplers cannot be accurately scaled using transmission line theory methods [71]. Despite similar performance characteristics, two of the optimized designs slightly violate the requirements concerning bandwidth. The obtained results are considered acceptable (especially having in mind complexity of the objective function used for their optimization). The average computational cost of BLC re-design is around 17.4 min of CPU-time.

4.5. Discussion and Measurements

The proposed design framework employs several mechanisms that make it unique compared to other database-inspired solutions considered in the literature [16,32,37,72]. First of all, it reduces the involvement of the user in the design process, merely to choose the specifications and substrate parameters. Moreover, due to integration of design steps into the TR loop, the method generates optimal designs at a cost that is substantially lower than conventional SBO approaches [32,53,54]. The difference is especially pronounced when high-fidelity EM models, used for evaluation of structure performance are efficient [53]. Finally, the algorithm supports adjustment of substrate parameters, and automated generation of the initial designs for the cell-level optimization.
As indicated in Table 6, the proposed framework provides improvement in terms of computational cost, flexibility of operation and design automation as compared to other algorithms based on the concept of re-using cells from the library [16,32,72]. For instance, in [16], the cells are represented only at the level of EM models, which negatively affects the cost of their tuning. Moreover, the method does not permit changing cell properties (substrate parameters, input impedance), which limits its flexibility in terms of cells application scenarios. Finally, the method in [16] is, at best, semi-automatic as it does not provide standardized mechanisms for cells/couplers tuning. In [72], a relatively low design cost is maintained using surrogate-assisted optimization. The method, however, does not provide mechanisms for rejection of structures with inferior performance. Instead, it yields a set of designs representing trade-off between selected performance figures to compare miniaturized circuits. On the other hand, in real-life one is normally interested in finding a single design that meets the specifications. In this regard, the redundancy of the generated solutions is undesirable from the standpoint of computational cost unless designer’s priorities are not clearly defined [45,72]. In [32], a methodology dedicated to design bandwidth-enhanced couplers is presented. Although the method offers mechanisms for adjusting the cells for broadband operation, their selection is still driven by the user. Furthermore, the method is expensive in terms of CPU-time due to application of numerically inefficient EM representation of coupler.
Despite using pre-defined cells, the results of Section 4.3 show that the proposed framework is capable of producing designs that are competitive (miniaturization-wise) even with manually developed BLCs. Table 7 provides benchmark of the coupler BA(3,1) against compact BLCs from the literature. For fair comparison, dimensions of structures have been expressed in terms of a guided wavelength λg calculated for the given center frequency and substrate parameters. Moreover, the bandwidth BW—expressed in percent to account for different operating frequencies of considered circuits—is defined for |S11| and |S41| being both below −20 dB. The structure of Section 4.3 offers not only substantial miniaturization rate of almost 80% but also maintains nearly 10% bandwidth, which is more than for majority of considered BLCs. It should be noted that the couplers of [18,20] have been miniaturized using fractal curves. For considered operation frequencies, practical application of fractals is limited only to thin substrate materials (h < 0.2 mm). The reason is that thin substrates are required to obtain relatively narrow low-impedance TLs, which can be folded by more than one iteration of the selected fractal shape.
As shown in Section 4.4, the proposed methodology can also be used for rapid re-design of BLCs w.r.t. given substrate parameters. Although, the techniques for scaling microwave and antenna structures at a cost corresponding to only a few EM simulations have been reported in the literature [71], the prerequisite for their application is availability of an inverse model. The latter is identified based on a set of individually optimized training points. Since the number of reference designs grows quickly with the number of figures being of interest from the structure re-design standpoint, construction of inverse surrogate might be impractical when it is not intended for multiple re-use. In this work, re-design of the structure for different substrates was performed at an average cost corresponding to about 19 high-fidelity EM model simulations (~17 min of CPU-time) which is low having in mind that, in [71], up to 13 independently optimized reference designs were used for constructing inverse surrogate dedicated for scaling of couplers w.r.t. only one figure of interest.
The coupler BA(3,1) of Section 4.3 has been fabricated and measured. Photograph of the manufactured structure is shown in Figure 15a, whereas comparison of simulation and measurement results is given in Figure 15b. The characteristics are in good agreement and the bandwidths are similar. Slight discrepancies between the obtained responses result from fabrication tolerances and the circuit assembly components which were not accounted for in the simulation models.
It is worth noting that low-pass properties of compact cells make them useful for construction of BLCs with harmonic suppression capability [15,24,28,31,32]. The mechanism is demonstrated in Figure 15a where broadband responses of the BA(3,1) structure are compared against characteristics of the conventional BLC. Although the CCs used for construction of BA(3,1) circuit have been optimized only for the center frequency of 1 GHz (cf. Section 3 and Section 4.3), the coupler features |S21| and |S31| below −20 dB within 2.2 GHz to 4.3 GHz range. At the same time, the conventional structure is not capable of suppressing unwanted frequencies.
In summary, the obtained numerical results indicate that the proposed framework can be used not only for unsupervised design of BLCs, but also for the comparison of compact topologies, and re-purposing the circuits for various substrates. Re-implementation of the SBO mechanisms around the TR-based algorithm ensures substantial reduction of optimization cost, as compared to conventional surrogate-assisted routines. The proposed mechanisms might be useful for design of other—more complex—structures with modular topology. Finally, the measurement results confirm the correctness of the used EM simulation models, whereas broadband simulations indicate the capability of the considered structures in suppressing harmonic frequencies.

5. Conclusions

A framework for low-cost automated design of compact branch-line couplers has been presented. The proposed methodology is based on a bottom-up design concept that involves automated selection of cells, bi-stage optimization of assembled BLC with respect to performance and size, and surrogate-assisted refinement of the coupler to the high-fidelity model level. To ensure the low operation cost, each step of the design process is embedded in a trust-region framework. Furthermore, the framework exploits feature-based representation of structure responses—at cell- and coupler-levels—and space-mapping-correction of low-fidelity models responses. It should be emphasized that the high-fidelity EM representation of BLC is used only during the final stage of the design process.
The performance of the proposed framework was demonstrated using three test cases concerning rapid optimization of compact coupler, size-oriented comparison of BLC geometries, and rapid re-design of coupler for various substrates. The average design cost using the presented methodology amounts to just 17.5 min of CPU-time. The approach proved to be useful for generating topologies featuring miniaturization rates of nearly 80%. The numerical performance and quality of the obtained designs have been positively benchmarked against other routines from the literature. The selected simulation results have been confirmed by measurements of fabricated BLC prototype.
Although the proposed methodology has been demonstrated only on BLC couplers, it is suitable for variety of circuits characterized by modular topologies such as antenna feeding networks, multi-band or broadband couplers, and others. Further work will focus on the adaptation of the framework for design of mentioned components.

Funding

This work was supported in part by the National Science Centre of Poland Grant 2018/31/B/ST7/02369.

Acknowledgments

The author would like to thank Slawomir Koziel, (Reykjavik University, Iceland) for making Keysight ADS and CST Microwave Studio available.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Tajik, A.; Alavijeh, A.S.; Fakharzadeh, M. Asymmetrical 4 × 4 butler matrix and its application for single layer 8 × 8 butler matrix. IEEE Trans. Ant. Prop. 2019, 67, 5372–5379. [Google Scholar] [CrossRef]
  2. Titz, D.; Ferrero, F.; Pilard, R.; Laporte, C.; Jan, S.; Ezzeddine, H.; Gianes, F.; Gloria, D.; Luxey, C. New wideband miniature branchline coupler on IPD technology for beamforming applications. IEEE Trans. Compon. Packag. Manuf. Technol. 2014, 4, 911–921. [Google Scholar] [CrossRef]
  3. Shahroury, F.R.; Wu, C. The design of low LO-Power 60-GHz CMOS quadrature-balanced self-switching current-mode mixer. IEEE Microw. Wirel. Comp. Lett. 2008, 18, 692–694. [Google Scholar] [CrossRef]
  4. de Padova, A.; Longhi, P.E.; Colangeli, S.; Ciccognani, W.; Limiti, E. Design of a GaN-on-Si single-balanced resistive mixer for Ka-band satcom. IEEE Microw. Wirel. Comp. Lett. 2019, 29, 56–58. [Google Scholar] [CrossRef]
  5. Gillick, M.; Robertson, I.D.; Joshi, J.S. Coplanar waveguide two-stage balanced MMIC amplifier using impedance-transforming lumped-distributed branchline couplers. IEE Proc. Microw. Ant. Prop. 1994, 141, 241–245. [Google Scholar] [CrossRef]
  6. Tseng, C.; Chang, C. Improvement of return loss bandwidth of balanced amplifier using metamaterial-based quadrature power splitters. IEEE Microw. Wirel. Comp. Lett. 2008, 18, 269–271. [Google Scholar] [CrossRef]
  7. Han, S.K.; Ha, Y.-S.; Kim, J.; Kang, K.-Y. Characteristics of high power capability for high temperature superconducting microwave multiplexer. Proc. Asia-Pacific Microw. Conf. 1997, 1, 109–112. [Google Scholar]
  8. Deferm, N.; Reynaert, P. CMOS Front Ends for Millimeter Wave Wireless Communication Systems; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  9. Qiu, S.; Wang, H.; Li, J.; Zhao, H.; Wang, Z.; Wang, J.; Wang, Q.; Plettemeier, D.; Bärhold, M.; Bauer, T.; et al. Towards wearable-inertial-sensor-based gait posture evaluation for subjects with unbalanced gaits. Sensors 2020, 20, 1193. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, W.; Yu, C.; Tu, C.; Lyu, Z.; Tang, J.; Ou, S.; Fu, Y.; Xue, Z. A survey on hand pose estimation with wearable sensors and computer-vision-based methods. Sensors 2020, 20, 1074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Altaf, A.; Seo, M. Dual-band circularly polarized dielectric resonator antenna for WLAN and WiMAX applications. Sensors 2020, 20, 1137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zainudin, N.; Latef, T.A.; Aridas, N.K.; Yamada, Y.; Kamardin, K.; Rahman, N.H.A. Increase of input resistance of a normal-mode helical antenna (NMHA) in human body application. Sensors 2020, 20, 958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tian, Q.; Lin, Y.; Guo, X.; Wang, J.; AlFarraj, O.; Tolba, A. An identity authentication method of a MIoT device based on radio frequency (RF) fingerprint technology. Sensors 2020, 20, 1213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zilani, T.A.; Al-Turjman, F.; Khan, M.B.; Zhao, N.; Yang, X. Monitoring movements of ataxia patient by using UWB technology. Sensors 2020, 20, 931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Hou, J.-A.; Weng, Y.-H. Design of compact 90 and 180 couplers with harmonic suppression using lumped-element bandstop resonators. IEEE Trans. Microw. Theory Tech. 2010, 58, 2932–2939. [Google Scholar] [CrossRef]
  16. Kurgan, P.; Filipcewicz, J.; Kitlinski, M. Development of a compact microstrip resonant cell aimed at efficient microwave component size reduction. IET Microw. Ant. Prop. 2012, 6, 1291–1298. [Google Scholar] [CrossRef]
  17. Tang, C.-W.; Chen, M.-G. Synthesizing microstrip branch-line couplers with predetermined compact size and bandwidth. IEEE Trans. Microw. Theory Tech. 2007, 55, 1926–1934. [Google Scholar] [CrossRef]
  18. Ghali, H.; Moselhy, T.A. Miniaturized fractal rat-race, branch-line, and coupled-line hybrids. IEEE Trans. Microw. Theory Tech. 2004, 52, 2513–2520. [Google Scholar] [CrossRef]
  19. Gai, C.; Jiao, Y.; Zhao, Y. Compact dual-band branch-line coupler with dual transmission lines. IEEE Microw. Wirel. Comp. Lett. 2016, 26, 325–327. [Google Scholar] [CrossRef]
  20. Chen, W.-L.; Wang, G.-M. Exact design of novel miniaturised fractal-shaped branch-line couplers using phase-equalising method. IET Microw. Ant. Prop. 2008, 2, 773–780. [Google Scholar] [CrossRef]
  21. Zhou, C.; Yang, H.Y.D. Design considerations of miniaturized least dispersive periodic slow-wave structures. IEEE Trans. Microw. Theory Tech. 2008, 56, 467–474. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, F. Miniaturized and harmonics-rejected slow-wave branchline coupler based on microstrip electromagnetic bandgap element. Microw. Opt. Technol. Lett. 2009, 51, 1080–1084. [Google Scholar] [CrossRef]
  23. Wang, J.; Wang, B.-Z.; Guo, Y.-X.; Ong, L.C.; Xiao, S. A compact slow-wave microstrip branch-line coupler with high performance. IEEE Microw. Wirel. Comp. Lett. 2007, 17, 501–503. [Google Scholar] [CrossRef]
  24. Kurgan, P.; Filipcewicz, J.; Kitlinski, M. Design considerations for compact microstrip resonant cells dedicated to efficient branch-line miniaturization. Microw. Opt. Technol. Lett. 2012, 54, 1949–1954. [Google Scholar] [CrossRef]
  25. Eccleston, K.W.; Ong, S.H.M. Compact planar microstripline branch-line and rat-race couplers. IEEE Trans. Microw. Theory Tech. 2003, 51, 2119–2125. [Google Scholar] [CrossRef] [Green Version]
  26. Liao, S.-S.; Sun, P.-T.; Chin, N.-C.; Peng, J.-T. A novel compact-size branch-line coupler. IEEE Microw. Wirel. Comp. Lett. 2005, 15, 588–590. [Google Scholar] [CrossRef]
  27. Liao, S.-S.; Peng, J.-T. Compact planar microstrip branch-line couplers using the quasi-lumped elements approach with nonsymmetrical and symmetrical T-shaped structure. IEEE Trans. Microw. Theory Tech. 2006, 54, 3508–3514. [Google Scholar] [CrossRef]
  28. Mondal, P.; Chakrabarty, A. Design of miniaturised branch-line and rat-race hybrid couplers with harmonics suppression. IET Microw. Ant. Prop. 2009, 3, 109–116. [Google Scholar] [CrossRef]
  29. Tseng, C.-H.; Chang, C.-L. A rigorous design methodology for compact planar branch-line and rat-race couplers with asymmetrical T-structures. IEEE Trans. Microw. Theory Tech. 2012, 60, 2085–2092. [Google Scholar] [CrossRef]
  30. Choi, K.-S.; Yoon, K.-C.; Lee, J.-Y.; Lee, C.-K.; Kim, S.-C.; Kim, K.-B.; Lee, J.-C. Compact branch-line coupler with harmonics suppression using meander T-shaped line. Microw. Opt. Technol. Lett. 2014, 56, 1382–1384. [Google Scholar] [CrossRef]
  31. Selga, J.; Vélez, P.; Coromina, J.; Fernández-Prieto, A.; Bonache, J.; Martín, F. Harmonic suppression in branch-line couplers based on slow-wave transmission lines with simultaneous inductive and capacitive loading. Microw. Opt. Technol. Lett. 2018, 60, 2374–2384. [Google Scholar] [CrossRef]
  32. Kurgan, P.; Koziel, S. Design of high-performance hybrid branch-line couplers for wideband and space-limited applications. IET Microw. Ant. Prop. 2016, 10, 1339–1344. [Google Scholar] [CrossRef]
  33. Lian, G.; Wang, Z.; He, Z.; Zhong, Z.; Sun, L.; Yu, M. A new miniaturized microstrip branch-line coupler with good harmonic suppression. Prog. EM Research Lett. 2017, 67, 61–66. [Google Scholar] [CrossRef] [Green Version]
  34. Mroczka, J. The cognitive process in metrology. Measurement 2013, 46, 2896–2907. [Google Scholar] [CrossRef]
  35. Kurgan, P.; Koziel, S. Fast surrogate-assisted simulation-driven optimization of compact microwave hybrid couplers. Eng. Opt. 2016, 48, 1109–1120. [Google Scholar] [CrossRef]
  36. Koziel, S.; Kurgan, P. Rapid design of miniaturised branch-line couplers through concurrent cell optimisation and surrogate-assisted fine-tuning. IET Microw. Ant. Prop. 2015, 9, 957–963. [Google Scholar] [CrossRef]
  37. Kurgan, P.; Koziel, S. Selection of circuit geometry for miniaturized microwave components based on concurrent optimization of performance and layout area. AEU-Int. J. Electron. Comm. 2019, 108, 287–294. [Google Scholar] [CrossRef]
  38. Kurgan, P.; Kitlinski, M. Novel doubly perforated broadband microstrip branch-line couplers. Microw. Opt. Technol. Lett. 2009, 51, 2149–2152. [Google Scholar] [CrossRef]
  39. Koziel, S.; Kurgan, P. Compact cell topology selection for size-reduction-oriented design of microstrip rat-race couplers. Int. J. Rf Microw. CAE 2018, 28, 1–9. [Google Scholar] [CrossRef]
  40. Kurgan, P.; Koziel, S. Surrogate-assisted EM-driven miniaturization of wideband microwave couplers by means of co-simulation low-fidelity models. Int. J. Rf Microw. CAE 2017, 281–287. [Google Scholar] [CrossRef]
  41. Caloz, C.; Itoh, T. Electromagnetic Metamaterials: Transmission Line Theory and Microwave Applications; Wiley-IEEE Press: Hoboken, NJ, USA, 2005. [Google Scholar]
  42. Nocedal, J.; Wright, S. Numerical Optimization, 2nd ed.; Springer: New York, NY, USA, 2006. [Google Scholar]
  43. Conn, A.; Scheinberg, K.; Vincente, L.N. Introduction to Derivative-Free Optimization; MPS-SIAM Series on Optimization; University City Science Center: Philadelphia, PA, USA, 2009. [Google Scholar]
  44. Conn, A.; Gould, N.I.M.; Toint, P.L. Trust-Region Methods; MPS-SIAM Series on Optimization; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2000. [Google Scholar]
  45. Deb, K. Multi-Objective Optimization Using Evolutionary Algorithms; John Wiley & Sons: New York, NY, USA, 2001. [Google Scholar]
  46. Gorissen, D.; Zhang, L.; Zhang, Q.J.; Dhaene, T. Evolutionary neuro-space mapping technique for modeling of nonlinear microwave devices. IEEE Trans. Microw. Theory Tech. 2011, 59, 213–229. [Google Scholar] [CrossRef] [Green Version]
  47. Cheng, Q.S.; Bandler, J.W.; Koziel, S. Tuning space mapping: The state of the art. Int. J. RF Microw. CAE 2012, 22, 639–651. [Google Scholar] [CrossRef]
  48. Bandler, J.W.; Cheng, Q.S.; Dakroury, S.A.; Mohamed, A.S.; Bakr, M.H.; Madsen, K.; Søndergaard, J. Space mapping: The state of the art. IEEE Trans. Microw. Theory Tech. 2004, 52, 337–361. [Google Scholar] [CrossRef]
  49. Forrester, A.I.J.; Keane, A.J. Recent advances in surrogate-based optimization. Prog. Aerospace Sci. 2009, 45, 50–79. [Google Scholar] [CrossRef]
  50. Queipo, N.V.; Haftka, R.T.; Shyy, W.; Goel, T.; Vaidynathan, R.; Tucker, P.K. Surrogate-based analysis and optimization. Prog. Aerospace Sci. 2005, 41, 1–28. [Google Scholar] [CrossRef] [Green Version]
  51. Simpson, T.W.; Pelplinski, J.D.; Koch, P.N.; Allen, J.K. Metamodels for computer-based engineering design: Survey and recommendations. Eng. Comput. 2001, 17, 129–150. [Google Scholar] [CrossRef] [Green Version]
  52. Koziel, S.; Leifsson, L. (Eds.) Surrogate-Based Modeling and Optimization: Applications in Engineering; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  53. Bekasiewicz, A.; Kurgan, P.; Kitlinski, M. A new approach to a fast and accurate design of microwave circuits with complex topologies. IET Microw. Ant. Prop. 2012, 6, 1616–1622. [Google Scholar] [CrossRef]
  54. Bekasiewicz, A.; Koziel, S.; Pankiewicz, B. Accelerated simulation-driven design optimisation of compact couplers by means of two-level space mapping. IET Microw. Ant. Prop. 2015, 9, 618–626. [Google Scholar] [CrossRef]
  55. Koziel, S.; Bekasiewicz, A. Fast simulation-driven feature-based design optimization of compact dual-band microstrip branch-line coupler. Int. J. Rf Microw. CAE 2016, 26, 13–20. [Google Scholar] [CrossRef]
  56. Zhang, C.; Feng, F.; Zhang, Q.; Bandler, J.W. Enhanced cognition-driven formulation of space mapping for equal-ripple optimisation of microwave filters. IET Microw. Ant. Prop. 2018, 12, 82–91. [Google Scholar] [CrossRef] [Green Version]
  57. Koziel, S.; Bandler, J.W. Rapid yield estimation and optimization of microwave structures exploiting feature-based statistical analysis. IEEE Trans. Microw. Theory Tech. 2015, 63, 107–114. [Google Scholar] [CrossRef]
  58. Pozar, D.M. Microwave Engineering, 4th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  59. Hong, J.-S.; Lancaster, M.J. Microstrip Filters for RF/Microwave Applications; John Wiley & Sons: Hoboken, NJ, USA, 2001. [Google Scholar]
  60. Rosłoniec, S. Liniowe Obwody Mikrofalowe: Metody Analizy i Syntezy; Wydawnictwa Komunikacji i Łączności: Warszawa, Poland, 1999. [Google Scholar]
  61. Selga, J.; Coromina, J.; Vélez, P.; Fernández-Prieto, A.; Bonache, J.; Martín, F. Miniaturised and harmonic-suppressed rat-race couplers based on slow-wave transmission lines. IET Microw. Ant. Prop. 2019, 13, 1293–1299. [Google Scholar] [CrossRef]
  62. Koziel, S.; Echeverría-Ciaurri, D.; Leifsson, L. Surrogate-based methods. In Computational Optimization, Methods and Algorithms; Koziel, S., Yang, X.S., Eds.; Series: Studies in Computational Intelligence; Springer: New York, NY, USA, 2011; pp. 33–60. [Google Scholar]
  63. Cheng, Q.S.; Bandler, J.W.; Koziel, S. An accurate microstrip hairpin filter design using implicit space mapping. IEEE Microw. Mag. 2008, 9, 79–88. [Google Scholar] [CrossRef]
  64. Bandler, J.W.; Cheng, Q.S.; Nikolova, N.K.; Ismail, M.A. Implicit space mapping optimization exploiting preassigned parameters. IEEE Trans. Microw. Theory Tech. 2004, 52, 378–385. [Google Scholar] [CrossRef]
  65. Gustrau, F. RF and Microwave Engineering: Fundamentals of Wireless Communications; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  66. Matthaei, G.L.; Young, L.; Jones, E.M.T. Microwave Filters, Impedance-Matching Networks, and Coupling Structures; Artech House: Norwood, MA, USA, 1980. [Google Scholar]
  67. Agilent (Keysight) ADS, Ver. 2011.10; Agilent Technologies: Santa Rosa, CA, USA, 2011.
  68. CST Microwave Studio, Ver. 2013; Dassault Systems: Vélizy-Villacoublay, France, 2013.
  69. Coello, C.A.C.; Lamont, G.B.; van Veldhuizen, D.A. Evolutionary Algorithms for Solving Multi-Objective Problems, 2nd ed.; Springer: New York, NY, USA, 2007. [Google Scholar]
  70. Tseng, C.-H.; Chen, H.-J. Compact Rat-Race Coupler Using Shunt-Stub-Based Artificial Transmission Lines. IEEE Microw. Wirel. Comp. Lett. 2008, 18, 734–736. [Google Scholar] [CrossRef]
  71. Koziel, S.; Bekasiewicz, A. Expedited geometry scaling of compact microwave passives by means of inverse surrogate modeling. IEEE Trans. Microw. Theory Tech. 2015, 63, 4019–4026. [Google Scholar] [CrossRef]
  72. Koziel, S.; Bekasiewicz, A. Rapid simulation-driven multiobjective design optimization of decomposable compact microwave passives. IEEE Trans. Microw. Theory Tech. 2016, 64, 2454–2461. [Google Scholar] [CrossRef]
  73. Cui, H.; Wang, J.; Li, J.-L. Compact Microstrip Branch-line Coupler with Wideband Harmonic Suppression. ACES J. 2012, 27, 766–771. [Google Scholar]
Figure 1. Comparison of microstrip equal-split branch-line couplers in terms of size: (a) A conventional structure [58], as well as circuits (b) miniaturized using standard symmetrical cells (size reduction: 30%) [53], (c) manually developed cells (reduction: 84%) [24], and (d) a combination of horizontal (black) and vertical (dark gray) CCs (reduction: 75%). Due to the implementation of couplers on different substrates, miniaturization rates are calculated based on dimensions expressed in terms of guided wavelength [61]. Numbers in circles next to the conventional design denote the enumeration scheme of BLC ports that is used in the work.
Figure 1. Comparison of microstrip equal-split branch-line couplers in terms of size: (a) A conventional structure [58], as well as circuits (b) miniaturized using standard symmetrical cells (size reduction: 30%) [53], (c) manually developed cells (reduction: 84%) [24], and (d) a combination of horizontal (black) and vertical (dark gray) CCs (reduction: 75%). Due to the implementation of couplers on different substrates, miniaturization rates are calculated based on dimensions expressed in terms of guided wavelength [61]. Numbers in circles next to the conventional design denote the enumeration scheme of BLC ports that is used in the work.
Sensors 20 03562 g001
Figure 2. Conceptual illustration of a compact cell (here, in horizontal configuration): (a) low-fidelity model Cl(x) and (b) high-fidelity EM-based model Ch(x). Response of the cell is controlled using input variables vector x, whereas parameters xgap and s = [ZH h εr t tanδ]T maintain separation between adjacent cells and control simulation conditions.
Figure 2. Conceptual illustration of a compact cell (here, in horizontal configuration): (a) low-fidelity model Cl(x) and (b) high-fidelity EM-based model Ch(x). Response of the cell is controlled using input variables vector x, whereas parameters xgap and s = [ZH h εr t tanδ]T maintain separation between adjacent cells and control simulation conditions.
Sensors 20 03562 g002
Figure 3. Calculation of structure responses (here, BLC) based on frequency characteristics of its building blocks. Note that all the steps enclosed in a light gray rectangle are controlled by the design framework of Section 3. Hence, the role of the user boils down to the choice of specifications and substrate-related parameters.
Figure 3. Calculation of structure responses (here, BLC) based on frequency characteristics of its building blocks. Note that all the steps enclosed in a light gray rectangle are controlled by the design framework of Section 3. Hence, the role of the user boils down to the choice of specifications and substrate-related parameters.
Sensors 20 03562 g003
Figure 4. Conceptual illustration of the BLC model (here, B(y) = Bh(y)). The structure response is obtained using transmission-line theory based on characteristics resulting from simulations of individual cells (here, represented in the form of high-fidelity EM models). The model supports cells implemented in the form of EM, equivalent-circuit, or surrogate models. Consistency of coupler geometry is ensured by using cells in horizontal and vertical configurations.
Figure 4. Conceptual illustration of the BLC model (here, B(y) = Bh(y)). The structure response is obtained using transmission-line theory based on characteristics resulting from simulations of individual cells (here, represented in the form of high-fidelity EM models). The model supports cells implemented in the form of EM, equivalent-circuit, or surrogate models. Consistency of coupler geometry is ensured by using cells in horizontal and vertical configurations.
Sensors 20 03562 g004
Figure 5. Surrogate model of individual CC—conceptual illustration. The correction layer parameters p = [xc α]T are adjusted in the course of parameter extraction process. Cell geometry and substrate parameters at the input port are maintained using x, xgap, and s parameters.
Figure 5. Surrogate model of individual CC—conceptual illustration. The correction layer parameters p = [xc α]T are adjusted in the course of parameter extraction process. Cell geometry and substrate parameters at the input port are maintained using x, xgap, and s parameters.
Sensors 20 03562 g005
Figure 6. Representation of structure characteristics in the form of feature points: (a) individual cell and (b) coupler. The squares represent feature points defined w.r.t. level (i.e., frequency point for the given level l). The circles are features defined w.r.t. frequency (i.e., response level for the given frequency point ω). Parameters ω0 and θ0 represent the center frequency and the desired phase shift.
Figure 6. Representation of structure characteristics in the form of feature points: (a) individual cell and (b) coupler. The squares represent feature points defined w.r.t. level (i.e., frequency point for the given level l). The circles are features defined w.r.t. frequency (i.e., response level for the given frequency point ω). Parameters ω0 and θ0 represent the center frequency and the desired phase shift.
Sensors 20 03562 g006
Figure 7. Family of CC characteristics (left-hand side) obtained for an arbitrarily selected line segment z in the search space and their corresponding feature-based response (right-hand side) defined w.r.t. level (here, minimum of S11 in frequency).
Figure 7. Family of CC characteristics (left-hand side) obtained for an arbitrarily selected line segment z in the search space and their corresponding feature-based response (right-hand side) defined w.r.t. level (here, minimum of S11 in frequency).
Sensors 20 03562 g007
Figure 8. Family of BLC characteristics (left-hand side) obtained for an arbitrarily selected line segment z in the search space and their corresponding feature-based responses (right-hand side) defined w.r.t. level (squares) and frequency (circles). Marker colors on frequency plot correspond to marker colors on the selected feature plots.
Figure 8. Family of BLC characteristics (left-hand side) obtained for an arbitrarily selected line segment z in the search space and their corresponding feature-based responses (right-hand side) defined w.r.t. level (squares) and frequency (circles). Marker colors on frequency plot correspond to marker colors on the selected feature plots.
Sensors 20 03562 g008
Figure 9. The database of compact cells used for automated BLC design. Dotted lines represent sub-sections of the substrate characterized by different permittivity/thickness pairs—each index k = 1, 2, 3, corresponds to εrk/hk of the given sub-section. Four horizontal (a) CH(1), (b) CH (2), (c) CH (3), (d) CH (4) and four vertical (e) CV(1), (f) CV (2), (g) CV (3), (h) CV (4) cells are implemented. The structures are designed in the form of an equivalent-circuit and EM models.
Figure 9. The database of compact cells used for automated BLC design. Dotted lines represent sub-sections of the substrate characterized by different permittivity/thickness pairs—each index k = 1, 2, 3, corresponds to εrk/hk of the given sub-section. Four horizontal (a) CH(1), (b) CH (2), (c) CH (3), (d) CH (4) and four vertical (e) CV(1), (f) CV (2), (g) CV (3), (h) CV (4) cells are implemented. The structures are designed in the form of an equivalent-circuit and EM models.
Sensors 20 03562 g009
Figure 10. Low-fidelity model reflection (–) and phase (···) responses of individual cells at the randomly generated initial (gray) and optimized (black) designs: (a) cell CH(4) and (b) cell CV(3).
Figure 10. Low-fidelity model reflection (–) and phase (···) responses of individual cells at the randomly generated initial (gray) and optimized (black) designs: (a) cell CH(4) and (b) cell CV(3).
Sensors 20 03562 g010
Figure 11. BLC responses vs. frequency at various stages of the design process: (a) low-fidelity model-based characteristics at y0 (gray) and y2* (black) and (b) EM-based responses of assembled coupler at y2* (gray) and y* (black). The characteristics are marked as follows: S11 (–), S21 (– –), S31 (–·), S41 (···).
Figure 11. BLC responses vs. frequency at various stages of the design process: (a) low-fidelity model-based characteristics at y0 (gray) and y2* (black) and (b) EM-based responses of assembled coupler at y2* (gray) and y* (black). The characteristics are marked as follows: S11 (–), S21 (– –), S31 (–·), S41 (···).
Sensors 20 03562 g011
Figure 12. Optimized BLC: (a) geometry and (b) comparison of the high-fidelity EM model BA responses: S11 (–), S21 (– –), S31 (–·), S41 (···) obtained using Keysight ADS (gray) and CST Studio (black).
Figure 12. Optimized BLC: (a) geometry and (b) comparison of the high-fidelity EM model BA responses: S11 (–), S21 (– –), S31 (–·), S41 (···) obtained using Keysight ADS (gray) and CST Studio (black).
Sensors 20 03562 g012
Figure 13. Frequency responses and geometries (in-scale) of optimized couplers (high-fidelity EM models) featuring acceptable performance: (a) BA(1,2), (b) BA(1,4), (c) BA(1,1), (d) BA(2,1), and (e) BA(3,1).
Figure 13. Frequency responses and geometries (in-scale) of optimized couplers (high-fidelity EM models) featuring acceptable performance: (a) BA(1,2), (b) BA(1,4), (c) BA(1,1), (d) BA(2,1), and (e) BA(3,1).
Sensors 20 03562 g013
Figure 14. Frequency characteristics and geometries (in-scale) of the BA(2,4) coupler re-designed to work on: (a) Arlon AD250C, (b) Arlon AD300C, and (c) FR-4 substrates.
Figure 14. Frequency characteristics and geometries (in-scale) of the BA(2,4) coupler re-designed to work on: (a) Arlon AD250C, (b) Arlon AD300C, and (c) FR-4 substrates.
Sensors 20 03562 g014
Figure 15. Branch-line coupler BA(3,1) of Section 4.3: (a) photograph of the manufactured circuit prototype, (b) comparison of simulated (gray) and measured (black) responses, and (c) comparison of the broadband frequency characteristics (black) against the responses of conventional BLC (gray).
Figure 15. Branch-line coupler BA(3,1) of Section 4.3: (a) photograph of the manufactured circuit prototype, (b) comparison of simulated (gray) and measured (black) responses, and (c) comparison of the broadband frequency characteristics (black) against the responses of conventional BLC (gray).
Sensors 20 03562 g015
Table 1. Surrogate-Assisted BLC Design: Benchmark Results.
Table 1. Surrogate-Assisted BLC Design: Benchmark Results.
Design Approach *Number of EvaluationsCostBW [MHz]ΔC @ω0 [dB]∠(S21/S31) @ω0 [°]Size [mm2]
BA (min)Total [min]
(i)84 BA84 (77)771903.1489.9266.8
(ii)1,#1240 Bl
2 BA
36.1 (33.1)
2 (1.8)
34.9N/AN/AN/AN/A
(ii) 2,#2749 Bl
6 BA
80 (73.3)
6 (5.5)
78.8N/AN/AN/AN/A
(iii) 12144 Bl
5 BA
62.4 (57.2)
5 (4.6)
61.81802.7590.1282.6
(iii) 2,$1371 Bl
3 BA
39.8 (36.6)
3 (2.8)
39.41903.2890.1299.7
(iv)296 Cl
8 Ch
55 Bl
7 BA
4.3 (3.9)
5.1 (4.7)
1.6 (1.5)
7 (6.4)
16.51803.2390.2276.5
This work221 Cl
6 Ch
39 Bl
4 BA
3.2 (3.0)
3.8 (3.5)
1.1 (1.0)
4 (3.7)
11.21703.0290.3276.5
* Initial design is set to y2*; average cost of obtaining y2* using method of Section 3 is ~3 min. 1 8 preassigned parameters that represent substrate thickness. 2 16 preassigned parameters that represent substrate thickness and permittivity. # Terminated due to divergence. $ Terminated due to divergence after yielding quasi-optimal solution.
Table 2. Performance comparison of optimized compact BLCs.
Table 2. Performance comparison of optimized compact BLCs.
Coupler TopologyBW [MHz]ΔC @ω0 [dB]∠(S21/S31) @ω0 [°]Dimensions [mm × mm]Size [mm2]Size Reduction [%] *PtS #
BA(3,2)760.6690.510.0 × 36.836982.1−26.0
BA(4,1)450.5090.121.6 × 25.755473.1−7.47
BA(3,4)890.4487.79.31 × 43.340380.5−4.99
BA(2,2)100.4090.014.1 × 30.643279.0−4.33
BA(2,4)560.4092.312.1 × 40.649576.0−3.02
BA(4,2)200.3390.217.5 × 23.741479.9−1.66
BA(4,4)890.3589.320.8 × 31.064468.8−0.28
BA(1,2)910.2889.914.9 × 35.252474.61.09
BA(1,4)1000.2187.117.2 × 40.269266.41.44
BA(1,1)940.2390.113.3 × 42.957372.21.55
BA(2,1)960.1589.912.8 × 41.553374.11.80
BA(3,1)990.2490.17.72 × 54.642179.62.20
* Calculated w.r.t. conventional 3dB (ΔC = 0) BLC operating at ω0 = 1 GHz and implemented on RO4003C with dimensions of 45.9 mm × 44.9 mm (footprint of 2061 mm2). # Only couplers with PtS > 0 feature acceptable performance w.r.t. specifications.
Table 3. Dimensions of couplers featuring acceptable performance.
Table 3. Dimensions of couplers featuring acceptable performance.
Coupler TopologyxH.1xH.2xH.3xH.4xH.5xH.6xV.1xV.2xV.3xV.4xV.5xV.6
BA(1,2)2.421.151.162.750.611.622.530.362.560.700.560.86
BA(1,4)1.222.201.142.910.571.6211.32.300.732.071.790.84
BA(1,1)2.080.721.332.820.361.275.034.870.820.621.32N/A *
BA(2,1)1.421.432.501.000.44N/A *5.124.640.611.260.96N/A *
BA(3,1)0.657.312.220.92N/A *N/A *8.022.490.600.580.97N/A *
* Parameter is not present in the structure model.
Table 4. Dimensions of the BA(2,4) coupler re-designed for different substrates.
Table 4. Dimensions of the BA(2,4) coupler re-designed for different substrates.
SubstratexH.1xH.2xH.3xH.4xH.5xH.6xV.1xV.2xV.3xV.4xV.5xV.6
Arlon AD250C1.871.872.501.031.00N/A *19.571.891.321.911.101.44
Arlon AD300C1.180.392.491.950.66N/A *14.331.201.011.341.241.76
FR-41.721.482.501.000.74N/A *13.701.661.391.721.081.03
* Parameter is not present in the structure model.
Table 5. Performance characteristics of the re-scaled BA(2,4) designs.
Table 5. Performance characteristics of the re-scaled BA(2,4) designs.
SubstrateBW [MHz]ΔC @ω0 [dB]∠(S21/S31) @ω0 [°]Dimensions [mm × mm]Size [mm2]
Arlon AD250C800.1288.615.5 × 58.9915
Arlon AD300C720.0488.413.5 × 42.2569
FR-4700.0988.314.0 × 48.9686
Table 6. Comparison of database-inspired design frameworks.
Table 6. Comparison of database-inspired design frameworks.
FrameworkSize ReductionCostFlexibilityAutomation
[16]+
[32]+++
[72]++++
This work+++++++
Table 7. Comparison of coupler BA(3,1) with benchmark structures.
Table 7. Comparison of coupler BA(3,1) with benchmark structures.
Structureω0 [GHz]εrh [mm]BW [%]Dimensions [mm × mm]Dimensions [λg × λg]Size [λg2]Miniaturization [%] *
[28]0.54.301.589.662.5 × 69.10.19 × 0.210.039235.9
[26]2.44.700.806.711.6 × 13.80.17 × 0.210.036141.0
[30]2.12.540.503.822.2 × 14.90.23 × 0.150.034144.2
[20]2.42.200.134.713.9 × 14.20.15 × 0.160.023661.3
[31]1.010.21.2711.314.6 × 20.20.13 × 0.180.022463.3
[73]1.02.940.768.126.6 × 30.90.14 × 0.160.021764.5
[18]2.43.480.1711.08.90 × 8.900.12 × 0.120.014077.1
[24]1.03.500.519.217.7 × 18.30.10 × 0.100.009983.8
[29]0.93.380.518.414.7 × 22.70.07 × 0.110.008086.9
This work1.03.380.819.97.70 × 56.40.04 × 0.300.013279.6
* Calculated w.r.t. conventional 3-dB BLC implemented on RO4003C with footprint of 0.0611 λg2.

Share and Cite

MDPI and ACS Style

Bekasiewicz, A. Low-Cost Automated Design of Compact Branch-Line Couplers. Sensors 2020, 20, 3562. https://doi.org/10.3390/s20123562

AMA Style

Bekasiewicz A. Low-Cost Automated Design of Compact Branch-Line Couplers. Sensors. 2020; 20(12):3562. https://doi.org/10.3390/s20123562

Chicago/Turabian Style

Bekasiewicz, Adrian. 2020. "Low-Cost Automated Design of Compact Branch-Line Couplers" Sensors 20, no. 12: 3562. https://doi.org/10.3390/s20123562

APA Style

Bekasiewicz, A. (2020). Low-Cost Automated Design of Compact Branch-Line Couplers. Sensors, 20(12), 3562. https://doi.org/10.3390/s20123562

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop