-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmicro.i
134 lines (118 loc) · 2.06 KB
/
micro.i
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
[Mesh]
type = GeneratedMesh
dim = 1
xmax = 0.5
nx = 25
coord_type = RSPHERICAL
[]
[Problem]
type = FEProblem
[]
[Variables]
[./Cs]
family = LAGRANGE
order = FIRST
initial_condition = 0.5
[../]
[]
[Kernels]
[./dcsdt]
type = TimeDerivative
variable = Cs
[../]
[./csdiff]
type = MatDiffusion
variable = Cs
diffusivity = 1.0e-5
[../]
[]
[GlobalParams]
# Kappa = 1.61358
#Kappa = 0.0
# Chi = 2.5
Cm = 0.1714785651793526
# Sigma = 4426.43
# K = 174.728
K2 = 2.286
#######################################
############################
#Omega = 0.0837492
[]
[BCs]
[./bv_flux]
execute_on = 'NONE'
type = ParticleBVPostBCKernel
variable = Cs
boundary = right
pps_c2 = c2_from_macro
pps_phi1 = phi1_from_macro
pps_phi2 = phi2_from_macro
MateChoice = 4
[../]
[]
[Preconditioning]
[./smp]
type = SMP
full = true
[../]
[]
[Executioner]
type = Transient
solve_type = PJFNK
#solve_type = NEWTON
petsc_options_iname = '-pc_type -ksp_gmres_restart -pc_factor_mat_solver_type'
petsc_options_value = ' lu 1601 mumps'
nl_rel_tol = 8.5e-08
nl_abs_tol = 1.5e-07
[./TimeStepper]
type = IterationAdaptiveDT
dt = 5.0e-5
optimal_iterations = 5
growth_factor = 1.1
cutback_factor = 0.5
[../]
dtmax = 1.0
end_time = 36000.0
#automatic_scaling=true
[]
[Outputs]
execute_on = 'INITIAL TIMESTEP_END'
print_linear_residuals = false
console = false
csv = false
exodus = false
#interval = 5
[]
[Postprocessors]
[./socp]
type = ElementAverageValue
variable = Cs
execute_on = 'TIMESTEP_END'
[../]
[./cs_surface]
type = SideAverageValue
variable = Cs
boundary = right
[../]
# [./J_from_macro]
# type = Receiver
# [../]
[./c2_from_macro]
type = Receiver
[../]
[./phi1_from_macro]
type = Receiver
[../]
[./phi2_from_macro]
type = Receiver
[../]
# [./mem]
# type= MemoryUsage
# [../]
# [./dofs]
# type= NumDOFs
# [../]
# [./elements]
# type= NumElems
# [../]
[]