sibc1stdisc.h File Reference

---

Detailed Description

This is sibc1st. The numproc solves the electric field vector wave equation in 3-dimensional space which is truncated by a cavity wall of finite thickness, thereby causing ohmic loss in the cavity boundary; the surface impedance is approximated as being constant over frequency, instead of varying as determined by the skin effect.

Copyright by Benedikt Oswald, 2002-2007, all rights reserved.

Objective: implementation file

Author:

Benedikt Oswald

Date:

2007 mar 20 ~ 13:40:00 by benedikt oswald

Warning:

none

Attention:

none required

Bug:

this is a research code.

Todo:

still big important things to do!

#include "compiler.h"
#include "gm.h"
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include "ugstruct.h"
#include "misc.h"
#include "ugdevices.h"
#include "commands.h"
#include "cmdint.h"
#include "cmdline.h"
#include "general.h"
#include "np.h"
#include "assemble.h"
#include "quadrature.h"
#include "shapes.h"
#include "namespace.h"
#include "aqhdMaterials.h"
#include "disconst.h"
#include "whitney.h"
#include "sibc1st.h"

Go to the source code of this file.

Functions
INT	AssembleSIBC1stTetrahedron (NP_BASE *base, ELEMENT *t, INT argc, char **argv)

---

Function Documentation

INT AssembleSIBC1stTetrahedron	(	NP_BASE *	base,
		ELEMENT *	t,
		INT	argc,
		char **	argv
	)

AssembleSIBC1stTetrahedron - compute local matrices and accumulate them

SYNOPSIS: static int AssembleABC1stTetrahedron(...)

PARAMETERS:

DESCRIPTION:

This function constructs the elemental matrices that discretize the curl-curl vector wave equation terminated by a 1st order absorbing boundary condition, cf. Jin, pp. 530ff.

[T_ij] e_tt + ([R_ij] + [Q_ij]) e_t + [S_ij] e = -F_i

/ where [T_ij] = epsilon | w1_i * w1_j dV /

/ where [R_ij] = sigma | w1_i * w1_j dV /

/ where [Q_ij] = Y_0 | [n x w1_i] * [n x w1_j] dS /

/ where [S_ij] = (1/mu) | ( Nabla x w1_i ) * ( Nabla x w1_j ) dV /

/ where [F_i] = | w1_i ) * d/dt J_0 dV /

with J_0 defined as the impressed or forced current with [T_ij] and [R_ij] being scaled versions of the same integral expression, with [Q_ij] being the matrix resulting from the inclusion of the 1st order ABC

SEE ALSO: numproc implemenetations in itdi.h and itdi.c

RETURN VALUE: INT

FE_OK if ok FE_NOT_TETRAHEDRON if not tetrahedron FE_BAD_TRAFO if bad trafo FE_NO_CONNECTION if no connection

prepare Gaussian quadrature & get quadrature function for evaluation of surface integrals for the 1st order ABC

The minus signs before val1 and val2 are only testwise, to find out if we have made a mistake in the analytical derivation of the SIBC integration into the electric field vector wave equation

---