www.mooseframework.org
PFCRFFMaterial.C
Go to the documentation of this file.
1 /****************************************************************/
2 /* MOOSE - Multiphysics Object Oriented Simulation Environment */
3 /* */
4 /* All contents are licensed under LGPL V2.1 */
5 /* See LICENSE for full restrictions */
6 /****************************************************************/
7 #include "PFCRFFMaterial.h"
8 
9 template <>
10 InputParameters
12 {
13  InputParameters params = validParams<Material>();
14  params.addRequiredParam<unsigned int>(
15  "num_L", "specifies the number of complex L variables will be solved for");
16  return params;
17 }
18 
19 PFCRFFMaterial::PFCRFFMaterial(const InputParameters & parameters)
20  : Material(parameters),
21  _M(declareProperty<Real>("M")),
22  _alpha_R_0(declareProperty<Real>("alpha_R_0")),
23  _alpha_I_0(declareProperty<Real>("alpha_I_0")),
24  _A_R_0(declareProperty<Real>("A_R_0")),
25  _A_I_0(declareProperty<Real>("A_I_0")),
26  _alpha_R_1(declareProperty<Real>("alpha_R_1")),
27  _alpha_I_1(declareProperty<Real>("alpha_I_1")),
28  _A_R_1(declareProperty<Real>("A_R_1")),
29  _A_I_1(declareProperty<Real>("A_I_1")),
30  _alpha_R_2(declareProperty<Real>("alpha_R_2")),
31  _alpha_I_2(declareProperty<Real>("alpha_I_2")),
32  _A_R_2(declareProperty<Real>("A_R_2")),
33  _A_I_2(declareProperty<Real>("A_I_2")),
34  _alpha_R_3(declareProperty<Real>("alpha_R_3")),
35  _alpha_I_3(declareProperty<Real>("alpha_I_3")),
36  _A_R_3(declareProperty<Real>("A_R_3")),
37  _A_I_3(declareProperty<Real>("A_I_3")),
38  _alpha_R_4(declareProperty<Real>("alpha_R_4")),
39  _alpha_I_4(declareProperty<Real>("alpha_I_4")),
40  _A_R_4(declareProperty<Real>("A_R_4")),
41  _A_I_4(declareProperty<Real>("A_I_4")),
42  _num_L(getParam<unsigned int>("num_L"))
43 {
44 }
45 
46 void
48 {
49  // Mobility
50  _M[_qp] = 1.0;
51 
52  // Alpha and A constants
53  if (_num_L == 3)
54  {
55  // alpha constants
56  _alpha_R_0[_qp] = 2.352788316033853;
57  _alpha_I_0[_qp] = 0.0;
58  _alpha_R_1[_qp] = -4.371217046300305;
59  _alpha_I_1[_qp] = 6.153993830413678;
60  _alpha_R_2[_qp] = -4.371217046300305;
61  _alpha_I_2[_qp] = -6.153993830413678;
62 
63  // A constants
64  _A_R_0[_qp] = -1.254832460194660e2;
65  _A_I_0[_qp] = 4.141043034348927e-15;
66  _A_R_1[_qp] = 24.798843718179786;
67  _A_I_1[_qp] = 37.678064436502760;
68  _A_R_2[_qp] = 24.798843718179786;
69  _A_I_2[_qp] = -37.678064436502760;
70  }
71  else if (_num_L == 5)
72  {
73  // alpha constants
74  _alpha_R_0[_qp] = 2.429134088464706;
75  _alpha_I_0[_qp] = 0.0;
76  _alpha_R_1[_qp] = -3.972333899872749;
77  _alpha_I_1[_qp] = 6.499130135847140;
78  _alpha_R_2[_qp] = -3.972333899872749;
79  _alpha_I_2[_qp] = -6.499130135847140;
80  _alpha_R_3[_qp] = -18.943264072194637;
81  _alpha_I_3[_qp] = 9.349446845430961;
82  _alpha_R_4[_qp] = -18.943264072194637;
83  _alpha_I_4[_qp] = -9.349446845430961;
84 
85  // A constants
86  _A_R_0[_qp] = -1.282478656880326e02;
87  _A_I_0[_qp] = 9.910190130869531e-15;
88  _A_R_1[_qp] = 34.212475550662354;
89  _A_I_1[_qp] = 42.274652746493430;
90  _A_R_2[_qp] = 34.212475550666770;
91  _A_I_2[_qp] = -42.274652746496530;
92  _A_R_3[_qp] = 10.501019149011636;
93  _A_I_3[_qp] = -2.363585468012575;
94  _A_R_4[_qp] = 10.501019149026910;
95  _A_I_4[_qp] = 2.363585467971611;
96  }
97 }
98 
99 /*
100  _alpha_R_0[_qp] = 2.412;
101  _alpha_I_0[_qp] = 0.0;
102  _alpha_R_1[_qp] = -18.62;
103  _alpha_I_1[_qp] = 9.968;
104  _alpha_R_2[_qp] = -51.8;
105  _alpha_I_2[_qp] = -18.58;
106  _alpha_R_3[_qp] = -104.6;
107  _alpha_I_3[_qp] = -19.21;
108  _alpha_R_4[_qp] = -3.88;
109  _alpha_I_4[_qp] = 6.545;
110 
111  // A constants
112  _A_R_0[_qp] = -63.1;
113  _A_I_0[_qp] = 0.0;
114  _A_R_1[_qp] = 12.52;
115  _A_I_1[_qp] = -3.607;
116  _A_R_2[_qp] = 3.88;
117  _A_I_2[_qp] = 0.7762;
118  _A_R_3[_qp] = 0.9984;
119  _A_I_3[_qp] = 0.1591;
120  _A_R_4[_qp] = 36.7;
121  _A_I_4[_qp] = 42.66;
122 
123  else if (_num_L == 5)
124  {
125  // alpha constants
126  _alpha_R_0[_qp] = 2.4887266073084095552303551812656223773956298828125;
127  _alpha_I_0[_qp] = 0.0;
128  _alpha_R_1[_qp] = -19.2733746470461682065433706156909465789794921875;
129  _alpha_I_1[_qp] = -9.18277447910810451503493823111057281494140625;
130  _alpha_R_2[_qp] = -19.2733746470461682065433706156909465789794921875;
131  _alpha_I_2[_qp] = 9.18277447910810451503493823111057281494140625;
132  _alpha_R_3[_qp] = -3.8695517424123173633176975272363051772117614746094;
133  _alpha_I_3[_qp] = 6.7955256217773678528715208813082426786422729492188;
134  _alpha_R_4[_qp] = -3.8695517424123173633176975272363051772117614746094;
135  _alpha_I_4[_qp] = -6.7955256217773678528715208813082426786422729492188;
136 
137  // A constants
138  _A_R_0[_qp] = -133.2927098034036816898151300847530364990234375;
139  _A_I_0[_qp] = 0.0;
140  _A_R_1[_qp] = 10.728194854990965367846911249216645956039428710938;
141  _A_I_1[_qp] = 2.5027746492227604946378960448782891035079956054688;
142  _A_R_2[_qp] = 10.728194854990965367846911249216645956039428710938;
143  _A_I_2[_qp] = -2.5027746492227604946378960448782891035079956054688;
144  _A_R_3[_qp] = 35.9742594324134188354946672916412353515625;
145  _A_I_3[_qp] = 45.6070815722133602321264334022998809814453125;
146  _A_R_4[_qp] = 35.9742594324134188354946672916412353515625;
147  _A_I_4[_qp] = -45.6070815722133602321264334022998809814453125;
148  }
149 
150 */
MaterialProperty< Real > & _alpha_R_1
MaterialProperty< Real > & _A_R_4
MaterialProperty< Real > & _alpha_R_3
MaterialProperty< Real > & _A_I_1
MaterialProperty< Real > & _A_R_0
MaterialProperty< Real > & _A_I_3
virtual void computeQpProperties()
PFCRFFMaterial(const InputParameters &parameters)
MaterialProperty< Real > & _alpha_I_4
MaterialProperty< Real > & _alpha_R_2
MaterialProperty< Real > & _A_I_4
MaterialProperty< Real > & _A_R_1
unsigned int _num_L
MaterialProperty< Real > & _A_R_3
MaterialProperty< Real > & _alpha_R_0
MaterialProperty< Real > & _M
MaterialProperty< Real > & _A_R_2
InputParameters validParams< PFCRFFMaterial >()
MaterialProperty< Real > & _A_I_2
MaterialProperty< Real > & _alpha_I_3
MaterialProperty< Real > & _alpha_R_4
MaterialProperty< Real > & _alpha_I_0
MaterialProperty< Real > & _A_I_0
MaterialProperty< Real > & _alpha_I_1
MaterialProperty< Real > & _alpha_I_2