In design consideration of microwave and millimeter-wave planar circuits, efficient simulation tools have been required for visualizing the operation characteristics of passive circuit components. Also, to clarify the signal propagation in the circuit system is very important in terms of engineering and educational effects. Microwave simulators are widely used as the supporting tool of microwave circuit design and development. For analyzing the circuit components on planar circuits, several numerical techniques are described about the 3-D electromagnetic field analysis modeling . Nowadays, the method of moments (MoM)  and the finite-difference time-domain method (FDTD) , are mostly used among microwave engineers, and the SNAP simulator  is SPICE like one. On the other hand, there are various microwave circuit components using parallel coupled lines and loop resonators . A square loop line with parallel coupled lines is used as ring resonators and traveling-wave loop directional filters. Here, in order to solve the transient problem of the microwave circuit components composed of transmission lines, we provide a numerical analysis method introduced systematic mixed even and odd modes modeling for coupled lines. This method based on a modified central difference method  can be applied to the time domain analysis of the microwave circuit components constructed with the parallel tightly coupled lines. For ring resonators and traveling-wave loop directional filters composed of such components ,, to clarify the circuit operation, the transient behavior of the voltage and current waves of transmission line networks will be represented with dynamic expressions. In order to elucidate the mechanism of the circuit operation, the behavior of propagating signals on the transmission line network has been represented in . In this article, we describe the utilization of the modified central difference method incorporating internal boundary treatments. Using this simulation technique, the time and frequency domain properties of the ring resonator and traveling-wave multi-loop directional coupler filters are analyzed and demonstrated.
For obtaining accurate operation characteristics of circuit components, simulation techniques with processing multiple reflections are required. The transient behavior of the voltage and current waves on the ring resonator and traveling-wave loop directional filters is demonstrated to obtain the power division and isolation properties of the directional couplers by using a numerical analysis model. The represented transient phenomena include the effects of multiple reflection waves caused by line discontinuities and parallel coupling. The operating mechanism of the circuit components can be easily confirmed by the visualization of the computed voltage and current solutions. The transient responses along the transmission lines are represented with the variations of the instantaneous voltage and current distributions including all the multiple reflections. Additionally, using the input and output responses extracted from the voltage and current solutions for a Gaussian pulse excitation, the frequency responses are obtained by using the fast Fourier transform.
2.1. Modeling of a single transmission line
The configuration of a ring resonator and traveling-wave loop directional filter is shown in Figure 1. By applying numerical even- and odd-mode analysis model to parts of parallel coupled lines, the solution of the whole system can be solved by single transmission line analysis, with the assumption of TEM mode propagation. The 1-D modeling of the single transmission line has reported in . Telegrapher's equations can be generally written by
The derived update equations in the case of lossless line are as follows:
where the suffices
2.2. Modeling of coupled transmission lines
Here, as shown in Figure 2, we carry out the even- and odd-mode numerical analysis for the parts of the coupled lines in Figure 1. In the considered parallel coupled lines model,
2.3. Boundary treatments
In this modeling, the boundary treatments of both the sides of the coupled line need to calculate at each of time step. For the voltages and currents at the boundaries of the coupled line part as shown in Figure 3, the following scattering matrices expressed by the reflection coefficient are used to compute the reflected and transmitted quantities.
Similarly, be done to the right side one. Further, the boundary treatments for each port of the coupled lines with mismatched load are carried out as follows:
Here, the initial conditions of each transmission line are given as follows:
As a final processing, the line voltage and current solutions of the coupled transmission lines are numerically computed by the following equations.
Preprocessing: Decision of line parameters
Postprocessing: Visualization of results
The solutions of the basic coupled lines as shown in Figure 2 are computed by the numerical even and odd modes analysis. The voltage and current solutions of the equivalent even- and odd-mode impedance lines of uniform parallel coupled lines are shown in Figure 4 (a)(b)(e)(f). Consequently, by synthesizing these solutions, namely using equations (11)-(14), the solutions for the drive line and sense line are obtained as shown in Figure 4 (c)(d)(g)(h). Subsequently, Figure 5 and Figure 6 show the input and output voltage waveforms extracted from the voltage solutions of the basic coupled lines with and without matched loads, respectively, from 0 to 800 time steps. The used design parameters as follows:
Transient solutions of the even-mode and odd-mode equivalent lines and the solution of the basic coupled line synthesized with their solutions: Voltage (upper) and Current (lower).
3. Wave propagation in the transmission line networks
In this section, several simulation results are shown for representing the wave propagation along the transmission line network. The transient voltage and current responses of the ring resonator and traveling-wave multi-loop directional coupler filters are analyzed and demonstrated with the time domain analysis.
3.1. Transient behavior of ring resonator filter
First, it is shown that the transient voltage and current responses of the ring resonator filter in Figure 1 (a) are demonstrated in the time domain. The designed parameters for this band stop filter (BSF) were used as follows:
3.2. Time responses of traveling-wave multi-loop directional filters
Next, we show the time responses for the traveling-wave double-loop 3dB directional coupler as shown in Figure 1 (b). Figure 10 shows the variation of instantaneous voltage and current distributions on the directional coupler in the steady state after 1500 time steps. The isolation at port 4 and the equal power division into ports 2 and 3 are observed for the 3-dB coupler at the center frequency. Note that the envelope denotes the standing waves on the transmission lines, in which they appear at only the coupled line area. In Figure 11, the input and output voltage waveforms at each port are represented for 2000 time steps. The CPU time took in a few seconds. The following designed parameters were used:
Figure 12 shows the time responses for the traveling-wave triple-loop directional filter. The variation of instantaneous voltage and current distributions on the directional filter in the steady state are illustrated after 3000 time steps. The isolation at ports 4 and 2, and full power into port 3 is observed for the directional filter at the center frequency
4. Frequency responses
As an application of this method, it is shown that frequency responses are obtained from the transient responses by the fast Fourier transform technique. Using the input and output responses extracted from the transient pulse solutions, we obtain each voltage
Figure 14 (a) illustrates the voltage pulse responses at each port for one loop resonator filter. Figure 14 (b) shows the return loss and insertion loss of the frequency characteristics, which has the characteristic of very sharp notch filter at the center frequency. The same line parameters in the previous section were used.
Figure 15 (a) illustrates the voltage pulse responses at each port for the traveling-wave double-loop 3-dB directional coupler. And, Figure 15 (b) shows the return loss and insertion loss of the frequency characteristics, which has the characteristic of band pass filter with 3dB coupling at the center frequency. Similarly, the result of the traveling-wave triple-loop directional filter is shown in Figure 16. It can be seen that from the frequency response, the sharp skirt characteristics compared with double-loop resonators is obtained. The simulated results of the present method are in good agreement with those of S-NAP design software package .
Thus, by using the presented simulation tool, the circuit designers can efficiently obtain the design parameters to improve the properties of the ring resonator and traveling-wave loop directional filters in both the time and frequency domains. Also, by the visual expression of the solutions, the signal propagation can be make it easy to understand the operation characteristics on the microwave circuit components composed of the coupled transmission lines. Finally, the implemented program is of compact algorithm by the Visual BASIC language, and considerably saves CPU time.
Visualizing the signal propagation in transmission line circuits is very important for understanding the operation mechanism of the circuit systems. It has been shown that the presented time domain simulation method can analyze and demonstrate the transient behaviors of the voltage and current waves of the ring resonator and traveling-wave loop directional filters by dynamic expression. Moreover, we have shown the frequency responses of the multi-loop coupled line filters by using the FFT technique. As a brief test tool, this method introduced the systematic mixed even and odd modes model for the part of parallel coupled lines is useful to confirm the operation characteristics of the circuit components consisted of loop resonator with coupled lines. This run can process on a small size PC system sufficiently.