navierStokes_TaylorHood.inc.2d 3.16 KB
Newer Older
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

ns->space->mesh: ${mesh_name}

% ============ PHYSICAL PARAMETERS - NS =================================

ns->viscosity1: 1  	% 	1   	0.1
ns->viscosity2: 10   
ns->density1: 100	%	100	1
ns->density2: 1000
main->sigma: 24.5		% 	24.5	1.96

ns->theta: 1
ns->force: [0.0, -0.98]		% gravitational force [m/s^2]

ns->force dirichlet bc: 0

ns->initial velocity: 0
ns->initial velocity value: 0.0
ns->non-linear term: 2 		% 0... u^old*grad(u_i^old), 1... u'*grad(u_i^old), 2... u^old*grad(u'_i)


ns->beta: 1

ns->exponent: 2


ns->poisson problem pertubation: 0 % 1..applySingularPertubation, 2..applySingularDBC, 0..nothing
ns->simple algorithm: 1
ns->calculate pressure: 0



% =========== OUTPUT ==============================================

ns->velocity->output->filename: ${output_folder}/velocity${output_postfix}_
ns->space->output[2]->filename: ${output_folder}/pressure${output_postfix}_
ns->space->output[0]->filename: ${output_folder}/u0_${output_postfix}_
ns->space->output[1]->filename: ${output_folder}/u1_${output_postfix}_

% ============= PROBLEM-SPACES ==================================
ns->space->components: 3

ns->space->polynomial degree[0]: ${polynomial-degree}+1
ns->space->polynomial degree[1]: ${polynomial-degree}+1
ns->space->polynomial degree[2]: ${polynomial-degree}

ns->space->dim: 2

% ================== SOLVER ======================================
ns->space->solver: petsc-navierstokes %umfpack % petsc-navierstokes
ns->space->solver->navierstokes->use old initial guess: 1
parallel->use zero start vector: 0
ns->space->solver->navierstokes->pressure component: 2
ns->space->solver->tolerance: 1.e-8
ns->space->solver->info: 8
ns->space->solver->left precon: no



% =================== OUTPUT =========================================
ns->velocity->output->ParaView animation: 1
ns->velocity->output->ParaView format: 1
ns->velocity->output->ParaView vector format: 1
ns->velocity->output->write vector as 3d vector: 1
ns->velocity->output->write every i-th timestep: 10
%ns->velocity->output->compression:  gzip
ns->velocity->output->append index: 1
ns->velocity->output->index length: 9
ns->velocity->output->index decimals: 7

ns->space->output[2]->ParaView animation: 1
ns->space->output[2]->ParaView format: 1
ns->space->output[2]->write every i-th timestep: 1
%ns->space->output[2]->compression:  gzip
ns->space->output[2]->append index: 1
ns->space->output[2]->index length: 9
ns->space->output[2]->index decimals: 7

ns->space->output[0]->ParaView animation: 1
ns->space->output[0]->ParaView format: 1
ns->space->output[0]->write every i-th timestep: 1
%ns->space->output[0]->compression:  gzip
ns->space->output[0]->append index: 1
ns->space->output[0]->index length: 9
ns->space->output[0]->index decimals: 7

ns->space->output[1]->ParaView animation: 1
ns->space->output[1]->ParaView format: 1
ns->space->output[1]->write every i-th timestep: 1
%ns->space->output[1]->compression:  gzip
ns->space->output[1]->append index: 1
ns->space->output[1]->index length: 9
ns->space->output[1]->index decimals: 7

ns->space->output->write serialization: 0
ns->space->output->serialization filename:	serial_ns
ns->space->input->read serialization:		0
ns->space->input->serialization filename:	serial_ns