コード例 #1
0
ファイル: pml.py プロジェクト: ixLeTh3E/pyfdtd
def pml(size, delta, thickness=20.0, mode='TMz'):
    """creates a perfectly matched layer as surounding boundary conditions"""
    # get patameter
    sizeX, sizeY = size
    deltaX, deltaY = delta

    # crate material
    shapeX, shapeY = int(sizeX / deltaX), int(sizeY / deltaY)
    sigma = {
        'electricX': numpy.zeros((shapeX, shapeY)),
        'electricY': numpy.zeros((shapeX, shapeY)),
        'magneticX': numpy.zeros((shapeX, shapeY)),
        'magneticY': numpy.zeros((shapeX, shapeY))
    }
    mask = numpy.zeros((shapeX, shapeY))

    # set constant
    c = constants.mu_0 / constants.epsilon_0

    # init PML
    sigmaMax = -(3.0 + 1.0) * constants.epsilon_0 * constants.c \
            * math.log(1.0e-5) / (2.0 * deltaX * thickness)

    for n in range(0, int(thickness + 1.0), 1):
        for j in range(0, int(shapeY), 1):
            sigma['electricY'][n, j] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticY'][n, j] = sigmaMax \
                    * math.pow(float(thickness - n - 0.5) / thickness, 3.0) * c
            mask[n, j] = 1.0

            sigma['electricY'][shapeX - 1 - n, j] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticY'][shapeX - 1 - n, j] = sigmaMax \
                    * math.pow(float(thickness - n + 0.5) / thickness, 3.0) * c
            mask[shapeX - 1 - n, j] = 1.0

        for i in range(0, int(shapeX), 1):
            sigma['electricX'][i, n] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticX'][i, n] = sigmaMax \
                    * math.pow(float(thickness - n - 0.5) / thickness, 3.0) * c
            mask[i, n] = 1.0

            sigma['electricX'][i, shapeY - 1 - n] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticX'][i, shapeY - 1 - n] = sigmaMax \
                    * math.pow(float(thickness - n + 0.5) / thickness, 3.0) * c
            mask[i, shapeY - 1 - n] = 1.0

    # create layer
    electric = (Material.epsilon(1.0, sigma['electricX']),
                Material.epsilon(1.0, sigma['electricY']))
    magnetic = (Material.mu(1.0, sigma['magneticX']),
                Material.mu(1.0, sigma['magneticY']))

    return electric, magnetic, mask
コード例 #2
0
ファイル: pml.py プロジェクト: 011120227/pyfdtd
def pml(size, delta, thickness=20.0, mode='TMz'):
    """creates a perfectly matched layer as surounding boundary conditions"""
    # get patameter
    sizeX, sizeY = size
    deltaX, deltaY = delta

    # crate material
    shapeX, shapeY = int(sizeX / deltaX), int(sizeY / deltaY)
    sigma = {
            'electricX': numpy.zeros((shapeX, shapeY)),
            'electricY': numpy.zeros((shapeX, shapeY)),
            'magneticX': numpy.zeros((shapeX, shapeY)),
            'magneticY': numpy.zeros((shapeX, shapeY))}
    mask = numpy.zeros((shapeX, shapeY))

    # set constant
    c = constants.mu_0 / constants.epsilon_0

    # init PML
    sigmaMax = -(3.0 + 1.0) * constants.epsilon_0 * constants.c \
            * math.log(1.0e-5) / (2.0 * deltaX * thickness)

    for n in range(0, int(thickness + 1.0), 1):
        for j in range(0, int(shapeY), 1):
            sigma['electricY'][n, j] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticY'][n, j] = sigmaMax \
                    * math.pow(float(thickness - n - 0.5) / thickness, 3.0) * c
            mask[n, j] = 1.0

            sigma['electricY'][shapeX - 1 - n, j] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticY'][shapeX - 1 - n, j] = sigmaMax \
                    * math.pow(float(thickness - n + 0.5) / thickness, 3.0) * c
            mask[shapeX - 1 - n, j] = 1.0

        for i in range(0, int(shapeX), 1):
            sigma['electricX'][i, n] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticX'][i, n] = sigmaMax \
                    * math.pow(float(thickness - n - 0.5) / thickness, 3.0) * c
            mask[i, n] = 1.0

            sigma['electricX'][i, shapeY - 1 - n] = sigmaMax \
                    * math.pow(float(thickness - n) / thickness, 3.0)
            sigma['magneticX'][i, shapeY - 1 - n] = sigmaMax \
                    * math.pow(float(thickness - n + 0.5) / thickness, 3.0) * c
            mask[i, shapeY - 1 - n] = 1.0

    # create layer
    electric = (Material.epsilon(1.0, sigma['electricX']),
            Material.epsilon(1.0, sigma['electricY']))
    magnetic = (Material.mu(1.0, sigma['magneticX']),
            Material.mu(1.0, sigma['magneticY']))

    return electric, magnetic, mask
コード例 #3
0
    def __init__(self, field, mode='TMz'):
        # save arguments
        self.field = field
        self.mode = mode

        # create sources
        self.source = Material(field.size, field.delta)

        # create listeners
        self.listener = []

        # create materials
        self.material = {}
        self.material['electric'] = Material(field.size, field.delta)
        self.material['magnetic'] = Material(field.size, field.delta)

        # add free space layer
        self.material['electric'][:, :] = Material.epsilon()
        self.material['magnetic'][:, :] = Material.mu()

        # add pml layer
        electric, magnetic, mask = pml(field.size, field.delta, mode=mode)
        self.material['electric'][mask] = electric
        self.material['magnetic'][mask] = magnetic
コード例 #4
0
ファイル: solver.py プロジェクト: psycow/pyfdtd
    def __init__(self, field, mode='TMz'):
        # save arguments
        self.field = field
        self.mode = mode

        # create sources
        self.source = Material(field.size, field.delta)

        # create listeners
        self.listener = []

        # create materials
        self.material = {}
        self.material['electric'] = Material(field.size, field.delta)
        self.material['magnetic'] = Material(field.size, field.delta)

        # add free space layer
        self.material['electric'][:, :] = Material.epsilon()
        self.material['magnetic'][:, :] = Material.mu()

        # add pml layer
        electric, magnetic, mask = pml(field.size, field.delta, mode=mode)
        self.material['electric'][mask] = electric
        self.material['magnetic'][mask] = magnetic