Data processing: measuring – calibrating – or testing – Measurement system in a specific environment – Electrical signal parameter measurement system
Reexamination Certificate
2001-12-20
2003-12-09
Hoff, Marc S. (Department: 2857)
Data processing: measuring, calibrating, or testing
Measurement system in a specific environment
Electrical signal parameter measurement system
C702S066000, C702S057000, C702S073000, C702S077000, C702S189000, C703S002000, C703S013000, C708S840000, C708S832000, C708S679000, C708S514000, C324S076210, C324S076120, C324S076140
Reexamination Certificate
active
06662125
ABSTRACT:
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to an electromagnetic wave analyzer and program which analyze electromagnetic waves by solving Maxwell's equations in the time and spatial domains with finite difference methods. More particularly, the present invention relates to an electromagnetic wave analyzer and program which divide a computational domain into non-uniform cells to solve Maxwell's equations.
2. Description of the Related Art
Today's computational electromagnetics exploits the finite-difference time-domain (FD-TD) method as a technique to analyze the transitional behavior of electromagnetic waves, using computers for numerical calculation. The FD-TD algorithm solves Maxwell's equations in time and spatial domains with finite difference methods. Because of its wide scope of applications and improved accuracy and efficiency, FD-TD solvers are used in various situations in these years. For details of FD-TD methods, see the following literatures: K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media,” IEEE Trans. Antennas and Propagation, Vol. 14, pp. 302-307, 1966; and A. Taflove, “Computational Electrodynamics,” MA, Artech House, 1995.
As mentioned above, the FD-TD algorithm solves Maxwell's equations including two rotations, with difference methods in both the time domain (t) and spatial domain (x, y, z). The next part will present specific Maxwell's equations for some different numbers of spatial dimensions. The following notation will be used to describe each problem:
Ex, Ey, Ez: electric field components in the x-, y-, and z-axis directions
Hx, Hy, Hz: magnetic field components in the x-, y-, and z-axis directions
&mgr;: permeability
∈: dielectric constant
&sgr;: electrical conductivity
First, in a one-dimensional space, the following two equations express how the electromagnetic waves propagate in the x-axis direction.
μ
∂
H
y
∂
t
=
∂
E
z
∂
x
(
1
)
ϵ
∂
E
z
∂
t
+
σ
E
z
=
∂
H
y
∂
x
(
2
)
In a two-dimensional space, the transverse magnetic (TM) waves propagating in the x-axis and y-axis directions are expressed as follows.
μ
∂
H
x
∂
t
=
-
∂
E
z
∂
y
(
3
)
μ
∂
H
y
∂
t
=
∂
E
z
∂
x
(
4
)
ϵ
∂
E
z
∂
t
+
σ
E
z
=
∂
H
y
∂
x
-
∂
H
x
∂
y
(
5
)
Similarly, the following partial differential equations give the behavior of transverse electric (TE) waves in a two-dimensional space.
ϵ
∂
E
x
∂
t
+
σ
E
x
=
∂
H
z
∂
y
(
6
)
ϵ
∂
E
y
∂
t
+
σ
E
y
=
-
∂
H
z
∂
x
(
7
)
μ
∂
H
z
∂
t
=
∂
E
x
∂
y
-
∂
E
y
∂
x
(
8
)
Lastly, electromagnetic waves in a three-dimensional space are expressed in the systems of six coupled partial differential equations as follows.
ϵ
∂
E
x
∂
t
+
σ
E
x
=
∂
H
z
∂
y
-
∂
H
y
∂
z
(
9
)
ϵ
∂
E
y
∂
t
+
σ
E
y
=
∂
H
x
∂
z
-
∂
H
z
∂
x
(
10
)
ϵ
∂
E
z
∂
t
+
σ
E
z
=
∂
H
y
∂
x
-
∂
H
x
∂
y
(
11
)
μ
∂
H
x
∂
t
=
∂
E
y
∂
z
-
∂
E
z
∂
y
(
12
)
μ
∂
H
y
∂
t
=
∂
E
z
∂
x
-
∂
E
x
∂
z
(
13
)
μ
∂
H
z
∂
t
=
∂
E
x
∂
y
-
∂
E
y
∂
x
(
14
)
The FD-TD method is applied to a given computational domain, which is normally discretized into small meshes or boxes, called “cells” or “differencing grids.” The size of those cells is one of the major factors that govern computational errors in FD-TD analysis. Electromagnetic waves may exhibit a sudden change in their transitional behavior in some particular region of the computational domain, often resulting in increased numerical errors. Such errors, however, can be reduced to an acceptable level if the selected cell size is sufficiently small.
Obviously, the simplest way of spatial discretization is to divide a given computational domain into evenly spaced grids. This uniform discretization, however, is likely to produce a large number of cells, thus consuming more computational resources including CPU time and memory space. To optimize the analytical model of interest, researchers often divide a computational domain into differently sized cells. This non-uniform discretization technique reduces the total number of cells, effectively alleviating the processing burden on the computer.
However, such a non-uniform computational domain has a known problem that the numerical error would increase at a boundary between two regions having different cell sizes. Such boundaries are referred to herein as the “cell-size boundaries.” Since the spatial step size changes at a cell-size boundary, the midpoint between two adjacent cells does not lie on that cell-to-cell boundary, unlike the cases of uniform discretization. The consequent displacement of assumed calculation points aggravates the accuracy of finite difference approximation, resulting in increased truncation errors.
The central differencing technique is known to be second-order accurate for equally sized cells. However, the accuracy of this technique is reduced to the first order at cell-size boundaries because the assumption of spatial uniformity fails in such particular space points. The resulting computational errors could sometimes reach a intolerable level for certain types of objects to be analyzed.
SUMMARY OF THE INVENTION
In view of the foregoing, it is an object of the present invention to provide an electromagnetic wave analyzer capable of handling non-uniform cells with smaller computational errors.
Another object of the present invention is to provide an electromagnetic wave analyzing program which divides a given computational domain into a plurality of cells and numerically analyzes the behavior of transitional electromagnetic fields on an individual cell basis, without introducing any additional errors.
To accomplish the first object, according to the present invention, there is provided an electromagnetic wave analyzer which divides a given computational domain into a plurality of cells and numerically analyzes transitional behavior of electromagnetic fields on an individual cell basis. This analyzer comprises the following elements: a cell size identification unit which identifies the uniformity of cells that surround each space point defined in the computational domain; a first calculation unit which calculates electromagnetic field components at a space point with a first calculation method when the cell size identification unit has identified the surrounding cells as being uniform in size; a second calculation unit which calculates electromagnetic field components at a space point with a second calculation method when the cell size identification unit has identified the surrounding cells as being non-uniform in size, where the second calculation method has smaller computational errors than those of the first calculation method; and a data output unit which outputs the electromagnetic field components at every space point which have been calculated by the first and second calculation units.
To accomplish the second object, the present invention provides a computer program product which enables a computer to divide a given computational domain into a plurality of cells and numerically analyze transitional behavior of electromagnetic fields on an individual cell basis. This program causes the com
LandOfFree
Electromagnetic wave analyzer and program for same does not yet have a rating. At this time, there are no reviews or comments for this patent.
If you have personal experience with Electromagnetic wave analyzer and program for same, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Electromagnetic wave analyzer and program for same will most certainly appreciate the feedback.
Profile ID: LFUS-PAI-O-3154516