FD-TD Mesh Generation

FD-TD Mesh Generation

Note: This document is under development.

Table Of Contents

Introduction
Mapping The Model Into The Mesh
Material Properties
Block Classification
Convert To 3D Cells
Convert To 2D Cells
Convert To 1D Fields
Active Blocks

Introduction

When simulation is run using an LC model, the first step is to convert data structures representing the ideal model geometry into a form suitable for the FD-TD electromagnetic solver. This process is referred to as mesh generation.

This document describes the mesh generation algorithm in LC.

LC Model

An LC model contains any number of adjacent, overlapping, or disjoint geometry blocks which define the material electromagnetic properties in that space. The blocks may be of just two primative shapes: rectangular or cylindrical. More complex shapes can be created by combining several primative blocks. For example, a conductor with a right-angle turn can be created with two adjacent rectangular blocks.

In addition to the geometry blocks which define the material properties of the model, active blocks are used for source excitation, loads, and calculating probe output values. These active blocks are mapped into the FD-TD mesh and interact with it during the simulation.

FD-TD Mesh

The FD-TD mesh is a rectangular space containing two cartesian grids, one for the electric field and one for the magnetic field. Each grid is composed of uniformly distributed nodes spaced by the cell width. The nodes define the corners of 3D cells, the faces of the 3D cells define 2D cells, and the edges define the electromagnetic vector field components. The two grids overlap, but are displaced by half of a cell width in each dimension. The nodes are indexed by (i,j,k) coordinates with the half cell displacement assumed. The electric field grid contains nodes from (0,0,0) to (imax,jmax,kmax), while the magnetic field grid contains nodes from (0,0,0) to (imax-1,jmax-1,kmax-1). Thus, the electric field components define the faces of the simulation domain, while the magnetic field components are all internal to the mesh.

Mapping The Model Into The Mesh

Finding The Model Bounds

Since the FD-TD mesh is a rectangular domain, while the region defined by the LC model can be of arbitrary shape, one of the first steps in mesh generation is to determine a bounding box for the LC model. This box will define the FD-TD domain. We choose a tight-fitting box that encloses all of the model geometry. To configure this box, the minimum and maximum coordinates within the model must be determined.

Since each model block is defined by the minimum and maximum coordinates of a bounding box, the simulation bounding box can be found by scanning the list of model blocks.

Converting To Standard Units

Since the model is defined in terms of arbitrary measurement units, it is convenient to convert the model coordinates into standard units of meters. It is easier to use this conversion factor the few times it is required than to modify the FD-TD update equations.

Model Origin

Since the mesh begins at coordinate (0,0,0) while the model bounding box begins at coordinate (xmin,ymin,zmin), then (-xmin,-ymin,-zmin) defines an offset from the model coordinate system origin and the mesh origin. This offset, along with the cell width, defines a transformation from model coordinates into mesh coordinates.

Mesh Size

Once the model bounds and the coordinate system transformation are known, then the mesh size can be determined by calculating the mesh maximum coordinates (imax,jmax,kmax) from the model maximum coordinates (xmax,ymax,zmax).

Material Properties

Block Materials

In an LC model, the geometry blocks themselves are assigned materials identified by a name. These materials are then assigned properties which are shared by all geometry blocks of the same material. This list of materials used by the geometry blocks is usually short, under ten materials.

Interface Materials

When any two geometry blocks are adjacent, the material properties of the interface surface must be determined. The list of block materials is augmented by an additional list of automatically generated interface materials.

The properties of an interface material is determined by the properties of the two base materials. The rules are as follows.

  1. If the interface is between a block and an undefined region, then the block material is used for the interface.

  2. If the interface is between two blocks of differing conductivities, then the material with the higher conductivity is used for the interface.

  3. If the conductivities are equal, then the permittivities are averaged, and a new interface material is created.

Coefficient List

The combined list of block and interface materials is converted into a list of FD-TD update equation coefficients. These values are used later during the mesh generation process to assign coefficients to the electric and magnetic fields.

Block Classification

Convert To 3D Cells

Convert To 2D Cells

Convert To 1D Fields

Active Blocks

Sources

Loads

Probes


LC Home
Copyright © Cray Inc.
Maintained by Kevin Thomas (kjt@cray.com).
Last modified Mon Aug 11 15:59:36 CDT 1997