1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
#include "WeakEquationPoisson.h"
#include "Material.h"
#include <iostream>
#include <fstream>
namespace ATC {
//==============================================================
// Class WeakEquationPoisson
//==============================================================
// R = B^T flux + N source = 0
// 0 = B^T perm B phi - N rho_+
//--------------------------------------------------------------
// Constructor
//--------------------------------------------------------------
WeakEquationPoisson::WeakEquationPoisson()
: WeakEquation(STATIC_PDE,ELECTRIC_POTENTIAL,1)
{}
//---------------------------------------------------------------------
void WeakEquationPoisson::B_integrand(
const FIELD_MATS &fields,
const GRAD_FIELD_MATS &grad_fields,
const Material * material,
DENS_MAT_VEC &flux) const
{
// flux = - perm D_x phi
material->electric_displacement(fields, grad_fields, flux);
flux[0] *= -1;
flux[1] *= -1;
flux[2] *= -1;
}
//---------------------------------------------------------------------
void WeakEquationPoisson::BB_tangent_coefficients(
const FieldName /* field */,
const FIELD_MATS &fields,
const Material* material,
DENS_MAT &coefs) const
{
material->permittivity(fields,coefs); // convention evalues (B^T perm B) > 0
}
//---------------------------------------------------------------------
bool WeakEquationPoisson::N_integrand(
const FIELD_MATS &fields,
const GRAD_FIELD_MATS & /* grad_fields */,
const Material * material,
DENS_MAT &flux) const
{
return material->electron_charge_density(fields,flux); // - N rho_+ = N n
}
//---------------------------------------------------------------------
void WeakEquationPoisson::NN_tangent_coefficients(
const FieldName /* field */,
const FIELD_MATS &fields,
const Material* material,
DENS_MAT &coefs) const
{
material->D_electron_charge_density(ELECTRIC_POTENTIAL,fields,coefs);
}
//==============================================================
// Class WeakEquationPoissonConstantRHS
//==============================================================
//--------------------------------------------------------------
// Constructor
//--------------------------------------------------------------
WeakEquationPoissonConstantRHS::WeakEquationPoissonConstantRHS()
: WeakEquationPoisson()
{}
//---------------------------------------------------------------------
bool WeakEquationPoissonConstantRHS::N_integrand(
const FIELD_MATS &fields,
const GRAD_FIELD_MATS & /* grad_fields */,
const Material * /* material */,
DENS_MAT &flux) const
{
FIELD_MATS::const_iterator nIter = fields.find(ELECTRON_DENSITY);
if (nIter != fields.end()) {
const DENS_MAT & n = nIter->second;
flux = -1.0*n;
return true;
}
return false;
}
};
|
