-
Notifications
You must be signed in to change notification settings - Fork 0
/
example_2_1.i
121 lines (111 loc) · 1.99 KB
/
example_2_1.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
# "Applied multiphase flow in pipes and flow assurance oil and gas production"
# Al-Safran, E., Brill, J. P., 2017
# Example 2.1: Determine the pipeline inlet pressure?
[Mesh]
type = GeneratedMesh
dim = 1
xmin = 0
xmax = 9000
nx = 900
[]
[MeshModifiers]
[./rotate]
type = Transform
transform = ROTATE
vector_value = '7 0 0'
[../]
[]
[UserObjects]
[./eos]
type = MoskitoEOS1P_IdealFluid
bulk_modulus = 2e+012
reference_temperature = 293.15
reference_enthalpy = 0
reference_density = 1000
[../]
[./viscosity]
type = MoskitoViscosityConst
viscosity = 0.001
[../]
[]
[Materials]
[./area0]
type = MoskitoFluidWell_1p1c
pressure = p
temperature = T
flowrate = q
well_direction = x
well_type = 1
eos_uo = eos
viscosity_uo = viscosity
well_diameter = 0.152
roughness_type = smooth
gravity = '0 -9.8 0'
output_properties = 'density'
outputs = exodus
[../]
[]
[BCs]
[./pbc]
type = DirichletBC
variable = p
boundary = right
value = 2.07e6
[../]
[./qbc]
type = DirichletBC
variable = q
boundary = right
value = 0.046296296
[../]
[]
[Variables]
[./T]
initial_condition = 293.15
[../]
[./p]
initial_condition = 2.07e6
[../]
[./q]
scaling = 1e-2
initial_condition = 0.046296296
[../]
[]
[Kernels]
[./Tkernel]
type = NullKernel
variable = T
[../]
[./pkernel]
type = MoskitoMass_1p1c
variable = p
flowrate = q
temperature = T
[../]
[./qkernel]
type = MoskitoMomentum_1p1c
variable = q
pressure = p
temperature = T
[../]
[]
[Preconditioning]
[./p2]
type = SMP
full = true
petsc_options_iname = '-pc_type -sub_pc_type -sub_pc_factor_shift_type -ksp_gmres_restart'
petsc_options_value = 'asm lu NONZERO 51'
[../]
[]
[Executioner]
type = Steady
l_tol = 1e-10
l_max_its = 50
nl_rel_tol = 1e-8
nl_abs_tol = 1e-9
nl_max_its = 50
solve_type = NEWTON
[]
[Outputs]
exodus = true
[]