cahnHilliard_shell.dat.3d 2.35 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
dimension of world: 3

% =================== MESH ================================
mesh_name: mesh
mesh->macro file name:      ./macro/fliese_20x20x2.3d
mesh->global refinements:   3
mesh->check:                0

lattice: 4*M_PI/sqrt(3)
mesh->scale mesh: 1
N: 10
x: ${N}*${lattice}
y: ${x}
z: ${x} / 10.0

mesh->dimension: [${x}, ${y}, ${z}]

#include "init/reinit.inc.2d"

% ============== USER-PARAMETER ==========================
ch->use reinit:            0
ch->epsilon:               0.5
ch->gamma:                 1
ch->double-well type:      1

% =========== OUTPUT ==============================================
ch->output->filename: ./output/ch3d_shell_

% ==================== TIMESTEPS ===============================
adapt->timestep:      0.1
adapt->max timestep:  1e+10
adapt->min timestep:  1e-10
adapt->start time:    0.0
adapt->end time:      10000.0

% ============= PROBLEM-SPACES ==================================
ch->components:           2
ch->polynomial degree[0]: 2
ch->polynomial degree[1]: 2
ch->dim:                  3
ch->mesh:                 mesh

% ================== SOLVER ======================================
ch->solver:             fgmres
ch->solver->backend:            cahn_hilliard
%ch->solver->petsc prefix:       ch
ch->solver->right precon:    ch
ch->solver->orthogonalization:  1
ch->solver->restart:            30
ch->solver->max iteration:      2000
ch->solver->tolerance:          1.e-8
ch->solver->info:               10

ch->solver->right precon->subsolver M->solver: cg
ch->solver->right precon->subsolver M->solver->backend: mtl
ch->solver->right precon->subsolver M->solver->tolerance: 1.e-3
ch->solver->right precon->subsolver M->solver->max iteration: 5
ch->solver->right precon->subsolver M->solver->left precon: diag

ch->solver->right precon->subsolver MpL->solver: cg
ch->solver->right precon->subsolver MpL->solver->backend: mtl
ch->solver->right precon->subsolver MpL->solver->tolerance: 1.e-3
ch->solver->right precon->subsolver MpL->solver->max iteration: 20
ch->solver->right precon->subsolver MpL->solver->left precon: diag

% =================== OUTPUT =========================================
ch->output->ParaView animation:        1
ch->output->ParaView format:           1
ch->output->write every i-th timestep: 10
ch->output->append index:              1
ch->output->index length:              9
ch->output->index decimals:            7