HIVE/input/AForm.i at main · farscape-project/HIVE · GitHub

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
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
!include Parameters.i

[Mesh]
[]

[Variables]
  [A]
    family = NEDELEC_ONE
    order = FIRST
  []
[]

[AuxVariables]
  [V]
    family = LAGRANGE
    order = FIRST
  []
  [P]
    family = MONOMIAL
    order = CONSTANT
  []
  [Vt]
    family = LAGRANGE
    order = FIRST
  []
  [E]
    family = NEDELEC_ONE
    order = FIRST
  []
  [B]
    family = NEDELEC_ONE
    order = FIRST
  []
[]

[Kernels]
  [curlcurlA_coil]
    type = CurlCurlField
    variable = A
    coeff = ${copper_reluctivity}
    block = coil
  []
  [curlcurlA_target]
    type = CurlCurlField
    variable = A
    coeff = ${steel_reluctivity}
    block = target
  []
  [curlcurlA_vacuum]
    type = CurlCurlField
    variable = A
    coeff = ${vacuum_reluctivity}
    block = vacuum_region
  []
  [dAdt_target]
    type = CoefVectorTimeDerivative
    variable = A
    coeff = ${steel_econductivity}
    block = target
  []
  [dAdt_coil_vacuum]
    type = CoefVectorTimeDerivative
    variable = A
    coeff = ${vacuum_econductivity}
    block = 'coil vacuum_region'
  []
  [gradV]
    type = CoupledGrad
    variable = A
    coupled_scalar_variable = V
    function = ${copper_econductivity}*sin(${voltage_wfrequency}*t)
    block = coil
  []
[]

[AuxKernels]
  [P]
    type = JouleHeatingAux
    variable = P
    vector_potential = A
    sigma = ${steel_econductivity}
    skip = ${skip_t_af}
    block = target
    execute_on = timestep_end
  []
  [Vt]
    type = ParsedAux
    variable = Vt
    coupled_variables = V
    use_xyzt = true
    expression = sin(${voltage_wfrequency}*t)*V
    block = coil
    execute_on = timestep_end
    enable = ${visualization}
  []
  [E]
    type = VectorTimeDerivativeAux
    variable = E
    coupled_vector_variable = A
    coeff = -1
    block = target
    execute_on = timestep_end
    enable = ${visualization}
  []
  [B]
    type = CurlAux
    variable = B
    coupled_vector_variable = A
    execute_on = timestep_end
    enable = ${visualization}
  []
[]

[BCs]
  [plane]
    type = VectorCurlPenaltyDirichletBC
    variable = A
    boundary = 'coil_in coil_out terminal_plane'
    penalty = 1e14
  []
[]

[Executioner]
  type = Transient
  solve_type = LINEAR
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_rtol -ksp_norm_type'
  petsc_options_value = 'hypre    ams            1e-9      preconditioned'
  start_time = 0.0
  end_time = ${end_t_af}
  dt = ${delta_t_af}
[]

[Outputs]
  exodus = ${visualization}
[]

[MultiApps]
  [VLaplace]
    type = FullSolveMultiApp
    input_files = VLaplace.i
    execute_on = initial
    clone_parent_mesh = true
  []
[]

[Transfers]
  [pull_potential]
    type = MultiAppCopyTransfer
    from_multi_app = VLaplace
    source_variable = V
    variable = V
  []
[]