A New Approach Based on Iterative Method for the Characterization of a Micro-Strip Line with Thick Copper Conductor

In this work, we applied two electromagnetic models for the characterization of a planar structure including a flat, thick copper conductor. Indeed the first model is consisted by modeling two metal ribbons without bulkiness, placed one above the other at a distance of h2 equal to the thickness of the thick conductor. This approach has been implemented and tested by the iterative method. The results of simulations have been compared with those calculated by the Ansoft HFSS software, and they are in good concordance, validating the method of analysis used. The second model is based on the calculation of the effective permittivity of the medium containing the thick conductor. This medium consists of a metallic region of complex relative permittivity ′′  ε , the rest of this medium is filled with air εr2 = 1. The effective permittivity εeff calculated from these two relative permittivity εr2 and ′′  ε . Comparing the simulation results of this new formulation of the iterative method with those calculated by the software Ansoft HFSS shows that they are in good matching which validates the second model.


Introduction
In Radio frequency, most devices are made in micro strip technology [1] [2].This technology became the best known and most used, this is due to its flat nature, ease of manufacturing, low cost, easy integration with circuits in the solid state, good heat dissipation structure that is used as good mechanical support etc.In general and in various research works concerning modeling and study of these structures, most researchers assumed that they possessed metallic conductors without thickness [3] [4].This simplifying assumption decreases the accuracy of the results of analytical methods used for the characterization of these structures.Several methods were used to characterize the influence of the thickness of the conductor used, such as mode-matching method [5], method of lines [6], spectral domain method [7] and conformal mapping method [8].For the study of structures with flat and thick conductors such as micro strip line in this document, we have taken the first model proposed in [6].A new formulation of the iterative method FWCIP (Fast Wave Concept Iterative Process) was made to extend the study of planar structures with thick flat conductor.This method is based on the concept of wave [9].It is developed for the simple planar layer modeling of structures [10] or multi-layer [11], and even arbitrarily complex shape.This is an easy method to implement due to the absence of test functions.It is always convergent and has considerable execution speed due to the FMT (Fast Modal Transform).

Iterative Method F.W.C.I.P
This method is well suited to the calculation of planar structures.In fact, TE and TM modes are used in the iterative method as digital basis of spectral domain in which the FFT.Subsequently, the concept of fast wave is introduced to reflect the boundary conditions and continuity of relationships in different parts of the interface Ω in terms of waves.The method involves determining an effective relationship to link the incident and reflected waves in different dielectric layers expressing thoughts in modal domain and the boundary conditions and continuity, expressed in terms of waves in spatial domain.The iterative process is then used to move from one field to another using the FMT thus accelerating the iterative process and then the convergence of the method.The use of the FMT requires the pixel description of the different regions of the dielectric interfaces [12].Thus the electromagnetic behavior of a planar structure will be described by writing the boundary conditions and continuity of the tangential fields on each pixel containing the interface circuitry to study.This integral formulation retains the advantages well known iterative methods including ease of implementation and speed of execution.

Skin Effect Phenomenon
An approach based on skin effect phenomenon for the modeling of a micro-strip line with thick copper conductor [3]- [5].In fact at high frequency the skin effect phenomenon occurs and the current only flows on the periphery of the conductor.According to the law of Biot and Savart a moving current in a conductor generates a magnetic field H around it.If an alternating current flows through a conductive wire, fillers vibrate and the field H varies which creates an induced current loop that opposes to the change of current in the conductor.Therefore, the sum of the alternating current with that of the loop always lowers the center of the driver while the two currents are added to the surface (see Figure 1).The analytical method used is based on the Fast Wave Concept Iterative Process (FWCIP).

Presentation of the Studied Structure
The studied structure is composed by a dielectric substrate of thickness H, on one of its two faces is deposited a metal strip of copper of thickness T, the other face is the ground plane of the structure.In fact the current in hyper frequency flows only on the surface of the conductor (see Figure 2).We modeled the thick metal tape by two metal strips without bulkiness placed one above the other at a distance 2 10 m h T = = µ .we have w T  This allows us to neglect both thick coasts of conductive structure.
Figure 3 presents the variation of skin depth of copper conductor as function of frequency.The value of the skin thickness can be done by applying the following analytical equation: with ρ = (1/σ).

Formulation of the Method
The theoretical formulation for the iterative method is based on determining the relationship between the incident waves " " defined in the spectral domain.The passage of the spatial domain to the spectral domain is using modal Fourier transform (FMT).The passage of the spectral domain to the space domain is using the transform inverse Fourier modal (MTF −1 ).These operations are done with repetitions until the convergence of the method.FMT and FMT −1 are used to speed up the computation time of the method.This evolution is initiated by the waves emitted by the excitation source on either side of the plane Ω 1 .Figure 4 summarizes the iterative method for planar structures with three layers incorporating different mediums.placed at plans (Ω 1 and Ω 2 ).ˆk Γ : Operator reflection ensuring the link between the incident waves and the reflected waves.It is defined in the spectral domain.It contains information on the housing walls and the relative permittivity of the different mediums of the structure, k ∈ {medium 1, medium 3}.ˆQ Γ : Diffraction Operator at each interface (Ω 1 and Ω 2 ).The evolution of iterations through the spectral domain to the space domain is done using the Fourier transform modal "FMT" which considerably reduces the calculation time.Modal Fourier transform requires the fragmentation of discontinuity planes (Ω 1 and Ω 2 ) in pixels and this so that the electromagnetic behavior of the overall circuit will be summarized by writing the boundary conditions and continuity of the tangential fields on each pixel.The iterative process stops when it reaches the convergence of results.
The terms below link the incident waves " The operators of diffraction  The flowchart in Figure 5 summarizes the evolution of the iterative method for a planar structure with three layers of different mediums.
• Diffraction Operator: ˆ1 Ω Γ For a source of bilateral excitation polarized in (oy), the overall diffraction operator is written from the diffraction operators in different regions of Ω 1 plane (metal region, source region of excitement, dielectric region): ) where: H s1 = 1 on the source and 0 elsewhere.H m1 = 1 on the metal and 0 elsewhere.H i1 = 1 on the dielectric and 0 elsewhere.where: H m2 = 1 on the metal and 0 elsewhere.H i2 = 1 on the dielectric and 0 elsewhere.

• Expression of the reflection operator: ˆk Γ
It is defined in the spectral domain and contains information about the nature of the housing and the relative permittivity of the medium 1 and 3 of the structure.It is expressed by the following relationship: For a top cover (or lower) placed at a distance h from Ω plan.
For an open circuit without top cover (or lower).

• Expression of the FMT
The Fourier transform in cosine and sine is defined by: The Fourier mode transform (FMT) is defined by: , , T : Passing modal operator in the area expressed by: ( ) • The reflection operator of the Quadruple: ˆQ Γ The reflection operator of the Quadruple is defined in layer 2 of the structure to be studied.It links the incident waves " 21 A  , 22 A  " the reflected waves " 21 B  , 22 B  " they the space domain to the spectral domain.The quadruple Q ensures the passage of plan Ω 1 to plan Ω 2 and inversely (see Figure 2).
According to the diagram in Figure 6 we can write: • The symmetry of the structure, allows us to write: After some mathematically manipulation, it is possible to determine the matrix: Plan

Convergence of the Method
This study begins by checking the convergence of results based on iterations.This is to optimize the calculation time and improve the accuracy of the method.Figure 7 shows that the real part of the structure of the input impedance, converges from 3500 iterations and the imaginary part converges from 2000 iterations, for a frequency f = 1 GHz.

Validation of the Method
The simulation result (Figure 8) shows a comparison between a micro-strip line with a thick conductor and between a microstrip line with a conductor without thickness.These results show that our model simulation results find are identical with those calculated by Ansoft HFSS software, for a thick conductor.These results show the effectiveness of the correction made by our Model to the iterative method (FWCIP).

Studied Structure
In this structure (Figure 9(a)) we have modeled the entire thick conductor.The upper and lower films of the conductor are modeled by two metal strips without thicknesses (δ = 0), respectively placed in Ω 1 and Ω 2 plans.The interior of the thick film conductor and both of the two vertical sides are modeled by effective permittivity ε eff characterizing the layer 2 of the structure (Figure 9(b)).It is calculated from the complex relative permittivity of the thick conductor ε ′′  and the relative permittivity ε r2 = 1 of the air filling the remainder of the layer 2 of the study structure.
• Calculus of the effective permittivity of the second layer: Figure 10 shows the equivalent model of the medium 2 located between the two planes Ω 1 and Ω 2 of Figure 9(b).
The permittivity of a metallic conductor is given by the following relationship:   The effective permittivity "ε eff " is calculated from the relative permittivity of the medium 2 (second layer located between the planes Ω 1 and Ω 2 of the study structure (Figure 10).This medium consists of air ε r2 = 1 and copper ε ′′  calculated by the relation (21).The effective permittivity "ε eff " is calculated by the relation with: V 1 : Volume of air which occupies the second layer.
V 2 : Volume of copper which occupies the second layer.S 1 : Surface of air which occupies the plan Ω 1 .S 2 : Surface of the copper occupies the plan Ω 1 .
( )  The analysis of the two models is analogous, simply changer ε r2 = 1 by ε eff who becomes in complex form, so we obtain: At the planes Ω 1 and Ω 2 diffractions operators are given by the following matrix: ) )

Convergence of the Method
This study begins by checking the convergence of results based on iterations.This is to optimize the calculation time and improve the accuracy of the method.Figure 11 shows that the real part of the structure of the input impedance, converges from 700 iterations and the imaginary part converges from 1500 iterations, for a frequency f = 1 GHz.

Validation of the Analytical Method
The simulation results (Figure 12) for a line with a thick conductor are in good agreement with those calculated by Ansoft HFSS software.For a micro-strip line with a driver without thickness the results obtained by the FWCIP method are shifted.This shows Indeed the use of exact model taking into account the thickness of the conductor seems necessary to improve the efficiency of the iterative method FWCIP.These models tested showed adapting iterative method considered in the calculation of complex structures.Indeed, through the new formulation of the method, we showed a net correction of the resonance frequency observed in the case of structures without thickness.These models provide almost the same results (see Figure 13), requiring multilayer modeling structures.
We have shown the efficiency of the correction allocated to the iterative method by the two modals proposed for the modeling of the planar structures integrating the thick conductors.The results which were found and compared to those which were calculated by the HFSS software well demonstrate the improvement allocated to the iterative method.

Conclusions
This work allowed taking stock of two electromagnetic models with which we had characterized a planar structure including a flat, thick copper conductor.In fact, the first model based on the phenomenon of skin effect allowed us to model the latter with two metallic ribbons without thicknesses, placed one above the other with a distance h 2 equal to the thickness of the thick conductor.Both sides of the conductor have summers neglected   since the width of the metal is strictly greater than its thickness.This is a simplifying assumption which has no effect on the results of the problem.This approach has been implemented and tested by the iterative method.Simulations results found were compared with those calculated by the software Ansoft HFFS, they were in good agreement, validating the method of analysis used.The second model is based on the calculation of the effective permittivity of the medium containing the thick conductor.This medium consists of a metallic region of complex relative permittivity, and the rest of this medium is filled with air ε r2 = 1.The effective permittivity ε eff calculated from these two relative permittivity ε r2 and ε ′′  .Comparing the results of simulations of this new for- mulation of the iterative method with those calculated by the software Ansoft HFFS shows that they are in good agreement which validates the model.Both models tested in the iterative method have improved accuracy.These models generate almost the same results and require multilayer structures to be implemented.
The different cases of the structures studied in this article allowed to highlight the potential of iterative FWCIP and suggest that it would be an essential tool in the global modeling of planar structures with thick con-

Figure 1 .
Figure 1.Principle of the skin effect.

Figure 3 .
Figure 3. Variation of skin thickness (according in terms as a function) to the frequency.

Figure 4 .
Figure 4. Synoptic diagram summarizing the iterative method for planar structures with three layers.
of circuit that being in the Ω 1 and Ω 2 plans.

Figure 2 (
b) defines these two planes of discontinuities.

Figure 5 .
Figure 5. Flowchart summarizes the evolution of the iterative method for a planar structure with three layers of different mediums.
) mn f α : Basic functions.It depends on the nature of the box.Y α : Mode admittance reduced to the level of Ω plan.

Figure 6 .
Figure 6.Equivalent electrical circuit of the planar structure including three layers of different medium.

Figure 8 .
Figure 8.Comparison between a micro-strip line without thickness, micro-strip line with thickness (HFSS) and equivalent first model (FWCIP).

Figure 9 .
Figure 9. Studied structure and equivalent model.

Figure 11 .
Figure 11.Convergence of Ze impedance based on iterations.

Figure 12 .
Figure 12.Structure study and simulation results.

Figure 13 .
Figure 13.Comparison between the two proposed models.