Modeling Metallic Objects in Wave Electromagnetics Problems

May 14, 2015

Metals are materials that are highly conductive and reflect an incident electromagnetic wave — light, microwaves, and radio waves — very well. When using the RF Module or the Wave Optics Module to simulate electromagnetics problems in the frequency domain, there are several options for modeling metallic objects. Here, we will look at the Impedance and Transition boundary conditions as well as the Perfect Electric Conductor boundary condition, offering guidance on when to use each one.

What Is a Metal?

When approaching the question of what a metal is, we can do so from the point of view of the governing Maxwell’s equations that are solved for electromagnetic wave problems. Consider the frequency-domain form of Maxwell’s equations:

\nabla \times \left( \mu_r^{-1} \nabla \times \mathbf{E} \right) – {-\frac{\omega^2}{c_0^2}} \left( \epsilon_r -\frac{i \sigma}{\omega \epsilon_0} \right) \mathbf{E}= 0

The above equation is solved in the Electromagnetic Waves, Frequency Domain interface available in the RF Module and the Wave Optics Module. This equation solves for the electric field, \mathbf{E}, at the operating (angular) frequency \omega = 2 \pi f. The other inputs are the material properties: \mu_r is the relative permeability, \epsilon_r is the relative permittivity, and \sigma is the electrical conductivity.

For the purposes of this discussion, we will say that a metal is any material that is both lossy and has a relatively small skin depth. A lossy material is any material that has a complex-valued permittivity or permeability or a non-zero conductivity. That is, a lossy material introduces an imaginary-valued term into the governing equation. This will lead to electric currents within the material, and the skin depth is a measure of the distance into the material over which this current flows.

At any non-zero operating frequency, inductive effects will drive any current flowing in a lossy material towards the boundary. The skin depth is the distance into the material within which approximately 63% of the current flows. It is given by:

\delta=\left[ \operatorname{Re} \left( \sqrt{i \omega \mu_0 \mu_r (\sigma + i \omega \epsilon_0 \epsilon_r)} \right) \right] ^{-1}

where both \mu_r and \epsilon_r can be complex-valued.

At very high frequencies, approaching the optical regime, we are near the material plasma resonance and do in fact represent metals via a complex-valued permittivity. But when modeling metals below these frequencies, we can say that the permittivity is unity, the permeability is real-valued, and the electrical conductivity is very high. So the above equation reduces to:

\delta=\sqrt{\frac{2}{\omega \mu_0 \mu_r \sigma }}

Before you even begin your modeling in COMSOL Multiphysics, you should compute or have some rough estimate of the skin depth of all of the materials you are modeling. The skin depth, along with your knowledge of the dimensions of the part, will determine if it is possible to use the Impedance boundary condition or the Transition boundary condition.

The Impedance Boundary Condition

Now that we have the skin depth, we will want to compare this to the characteristic size, L_c, of the object we are simulating. There are different ways of defining L_c. Depending on the situation, the characteristic size can be defined as the ratio of volume to surface area or as the thickness of the thinnest part of the object being simulated.

Let’s consider an object in which L_c \gg \delta. That is, the object is much larger than the skin depth. Although there are currents flowing inside of the object, the skin effect drives these currents to the surface. So, from a modeling point of view, we can treat the currents as flowing on the surface. In this situation, it is appropriate to use the Impedance boundary condition, which treats any material “behind” the boundary as being infinitely large. From the point of view of the electromagnetic wave, this is true, since L_c \gg \delta means that the wave does not penetrate through the object.

Example of a metallic object with a very small skin depth.
The Impedance boundary condition is appropriate if the skin depth is much smaller than the object.

With the Impedance boundary condition (IBC), we are able to avoid modeling Maxwell’s equations in the interior of any of the model’s metal domains by assuming that the currents flow entirely on the surface. Thus, we can avoid meshing the interior of these domains and save significant computational effort. Additionally, the IBC computes losses due to the finite conductivity. For an example of the appropriate usage of the IBC and a comparison with analytic results, please see the Computing Q-Factors and Resonant Frequencies of Cavity Resonators tutorial.

The IBC becomes increasingly accurate as L_c / \delta \rightarrow \infty; however, it is accurate even when L_c / \delta \gt > 10 for smooth objects like spheres. Sharp-edged objects such as wedges will have some inaccuracy at the corners, but this is a local effect and also an inherent issue whenever a sharp corner is introduced into the model, as discussed in this previous blog post.

Now, what if we are dealing with an object that has one dimension that is much smaller than the others, perhaps a thin film of material like aluminum foil? In that case, the skin depth in one direction may actually be comparable to the thickness, so the electromagnetic fields will partially penetrate through the material. Here, the IBC is not appropriate. We will instead want to use the Transition boundary condition.

The Transition Boundary Condition

The Transition boundary condition (TBC) is appropriate for a layer of conductive material with a thickness relatively smaller than the characteristic size, and curvature, of the objects being modeled. The TBC can be used even if the thickness is many times greater than the skin depth.

The TBC takes the material properties as well as the thickness of the film as inputs, computing an impedance through the thickness of the film as well as a tangential impedance. These are used to relate the current flowing on the surface of either side of the film. That is, the TBC will lead to a drop in the transmitted electric field.

From a computational point of view, the number of degrees of freedom on the boundary is doubled to compute the electric field on either surface of the TBC, as shown below. Additionally, the total losses through the thickness of the film are computed. For an example of using this boundary condition, see the Beam Splitter tutorial, which models a thin layer of silver via a complex-valued permittivity.

Diagram of a metallic object with surface currents.
The Transition boundary condition computes a surface current on either side of the boundary.

Adding Surface Roughness

So far, with both the TBC and the IBC, we have assumed that the surfaces are perfect. A planar boundary is assumed to be geometrically perfect. Curved boundaries will be resolved to within the accuracy of the finite element mesh used, the geometric discretization error, as discussed here.

An illustration of the surface currents on both a rough and smooth surface.
Rough surfaces impede current flow compared to smooth surfaces.

All real surfaces, however, have some roughness, which may be significant. Imperfections in the surface prevent the current from flowing purely tangentially and effectively reduce the conductivity of the surface (illustrated in the figure above). With COMSOL Multiphysics version 5.1, this effect can be accounted for with the Surface Roughness feature that can be added to the IBC and TBC conditions.

For the IBC, the input is the Root Mean Square (RMS) roughness of the surface height. For the TBC, the input is instead given in terms of the RMS of the thickness variation of the film. The magnitude of this roughness should be greater than the skin depth, but much smaller than the characteristic size of the part. The effective conductivity of the surface is decreased as the roughness increases, as described in “Accurate Models for Microstrip Computer-Aided Design” by E. Hammerstad and O. Jensen. There is a second roughness model available, known as the Snowball model, which uses the relationships described in The Foundation of Signal Integrity by P. G. Huray.

The Perfect Electric Conductor Boundary Condition

It is also worth looking at the idealized situation — the Perfect Electric Conductor (PEC) boundary condition. For many applications in the radio and microwave regime, the losses at metallic boundaries are quite small relative to the other losses within the system. In microwave circuits, for example, the losses in the dielectric substrate typically far exceed the losses at any metallization.

The PEC boundary condition is a surface without loss; it will reflect 100% of any incident wave. This boundary condition is good enough for many modeling purposes and can be used early in your model-building process. It is also sometimes interesting to see how well your device would perform without any material losses.

Additionally, the PEC boundary condition can be used as a symmetry condition to simplify your modeling. Depending on your foreknowledge of the fields, you can use the PEC boundary condition, as well as its complement — the Perfect Magnetic Conductor (PMC) boundary condition — to enforce symmetry of the electric fields. The Computing the Radar Cross Section of a Perfectly Conducting Sphere tutorial illustrates the use of the PEC and PMC boundary conditions as symmetry conditions.

Lastly, COMSOL Multiphysics also includes Surface Current, Magnetic Field, and Electric Field boundary conditions. These conditions are provided primarily for mathematical completeness, since the currents and fields at a surface are almost never known ahead of time.

Summary

In this blog post, we have highlighted how the Impedance, Transition, and Perfect Electric Conductor boundary conditions can be used for modeling metallic surfaces, helping to identify situations in which each should be used. But, what if you cannot use any of these boundary conditions? What if the characteristic size of the parts you are simulating are similar to the skin depth? In that case, you cannot use a boundary condition. You will have to model the metal domain explicitly, just as you would for any other material. This will be the next topic we focus on in this series, so stay tuned.


Comments (14)

Leave a Comment
Log In | Registration
Loading...
Ted Cook
Ted Cook
May 14, 2015

Is the j in Maxwell’s equation above the imaginary number? It would seem to be from this sentence, “lossy material introduces an imaginary-valued term into the governing equation,” and the statement that conductivity creates loss (and no mention of j as in input parameter). But then the skin depth equation uses imaginary i. Also, J is often used for currents in electrodynamics, so it is confusing.

Walter Frei
Walter Frei
May 14, 2015

Thank you Ted, The equations have been adjusted for consistency.

Göran Eriksson
Göran Eriksson
May 19, 2015

Thanks for a nice overview of how to model conducting bodies and layers. Just to clarify the usability of the Transition Boundary Condition, I would like to point out that it is in fact accurate for any ratio between the skin depth and the thickness of the layer. The only two conditions that have to be fulfilled are that (i) the thickness should be smaller than the typical radius of curvature and (ii) the material should be a good conductor. For those interested in more details about the TBC used in COMSOL MP I can recommend reading a conference paper, “Efficient 3D Simulation of Thin Conducting Layers of Arbitrary Thickness”, that I presented at the 20007 IEEE International Conference on Electromagnetic Compatability. There one can find some text on the theoretical and historical background, as well as a couple of validation and demonstration cases. I have a feeling than many COMSOL users are not so sure about when the TBC is applicable. It may therefore come as a positive surprise to see that it can in fact be used for both electrically thin and thick layers, as well as all intermediate cases.

Walter Frei
Walter Frei
May 20, 2015

Thank you Goran, I’ve revised the article with your comments.
Best, Walter

Stephan Savarese
Stephan Savarese
September 10, 2015

Thanks Walter fot this great overview ! How about modeling multiple thin layers using a TBC ? Or else is another method available with the RF module ? For instance: imagine two layers, copper and a less conductive material, say nylon, 10 micron each, on a PCB board at 2.4 GHz.

Walter Frei
Walter Frei
September 22, 2015

Dear Stephan,

Thank you for the comment. Although it is a bit early for me to say anything specific, please stay tuned to application example updates and product updates that are coming in the future.

Best Regards,
Walter

Omar Mohamed
Omar Mohamed
August 31, 2016

Hello Walter Frei,

Thank you for your post. I have a question about modeling metal nano-particle, which have size of 300nm.

Which boundary condition is suitable if I perform simulation at optical frequencies?

Thank you.
Omar

Walter Frei
Walter Frei
September 2, 2016

Dear Omar,

You will want to follow this example:
https://www.comsol.com/model/optical-scattering-by-gold-nanospheres-14697
Note that these modeling instructions can be followed using either the RF or the Wave Optics Modules. You may also find this blog useful:
https://www.comsol.com/blogs/modeling-of-materials-in-wave-electromagnetics-problems/

Anisuzzaman Boni
Anisuzzaman Boni
August 2, 2017

Dear sir,
Thanks for your post. I would like to ask for your suggestion for my work.

I am trying to find band diagram of a photonic crystal.My unit cell has periodic boundary condtion in x an dy direction .But z direction is open and also kz=0( no propagation).wave is only propagating in x and y direction.

PML and sacttering boundary condition has been used in top and bottom(open boundary).But the problem is the wave is reflecting back to the computational domain from the Sio2 substrate side.I have tried to increase the thickness and absorption coefficient(scaling factor).But no improvement has been found.

Could i use Impedence boundary condition at Sio2 boundary as u mention which treats any material “behind” the boundary as being infinitely large(is it only appropriate fro metal) .

Best Regards
Anisuzzaman Boni

Caty Fairclough
Caty Fairclough
August 4, 2017

Hi Anisuzzaman,

Thanks for your comment!

For questions related to your modeling, please contact our Support team.

Online Support Center: https://www.comsol.com/support
Email: support@comsol.com

kemerov kreshna
kemerov kreshna
June 28, 2021

Hi
Do you have any idea on simulating a selective frequency surface using split ring resonators but instead of pec,we use a metal(ie copper)?

Brianne Christopher
Brianne Christopher
June 28, 2021

Hello Kemerov,

Thank you for your comment. Please contact support@comsol.com for help with your modeling problem.

Best regards,
Brianne

Umid Jamolov
Umid Jamolov
September 8, 2023

Thank you Walter for this blog post. I was wondering if there is any simple way to model laminate iron (0.5mm thickness) stack using any boundary conditions in MHz. Is layered boundary condition can be used fo this purpose? Since anisotropic material models are not possible to simulate using impedance boundary condition, i am searching ways to simplify the problem.

Walter Frei
Walter Frei
September 8, 2023 COMSOL Employee

Hello Umid,
Thank you for raising this point. In the time since this blog was written, we have also included a Layered Impedance and Layered Transition boundary condition. There features were introduced in version 6.0 and 6.1. That being said, at MHz frequencies, the skin depth will be sufficiently less than 0.5mm, so I think that you could model solely the outside envelope of this stack. Of course, more detailed discussion is likely best to address via Support.
Best Regards,

EXPLORE COMSOL BLOG