cahnHilliard_dd.dat.2d 3.14 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
dimension of world: 2

% =================== MESH ================================
mesh_name: mesh
mesh->macro file name:      ./macro/macro.square.fine.2d
%mesh->periodic file:        ./macro/square2d.per
mesh->global refinements:   0
mesh->check:                0
parallel->pre refine: 0

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

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

% =========== INPUT ==============================================
ch->phasefield->filename: ./phase.arh
ch->phasefield->perturbation: 1.e-3

% =========== OUTPUT ==============================================
ch->space->output->filename: ./output/ch2d_dd_

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

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

% ================== SOLVER ======================================
ch->space->solver:        petsc-ch
%ch->space->solver->backend:            mtl
ch->space->solver->petsc prefix:       ch_
%ch->space->solver->ksp: -ch_ksp_type fgmres
%ch->space->solver->ksp: -mass_ksp_max_it 10 -mass_ksp_type fgmres    -laplace_ksp_max_it 10 -laplace_ksp_type fgmres
%ch->space->solver->ksp: -ch_ksp_type=fgmres
ch->space->solver->use old initial guess: 1
%ch->space->solver->ksp: -mass_ksp_max_it 15 -laplace_ksp_max_it 1 -laplace_ksp_type richardson -laplace_pc_type lu -laplace_pc_factor_mat_solver_package mumps
%ch->space->solver->ksp: -mass_ksp_max_it 1 -mass_ksp_type richardson -mass_pc_type lu -mass_pc_factor_mat_solver_package mumps -laplace_ksp_max_it 1 -laplace_ksp_type richardson -laplace_pc_type lu -laplace_pc_factor_mat_solver_package mumps
%ch->space->solver->orthogonalization:  1
%ch->space->solver->left precon: bjacobi
%ch->space->solver->right precon:       bjacobi
ch->space->solver->restart:            50
ch->space->solver->max iteration:      2000
%ch->space->solver->relative tolerance: 0
ch->space->solver->tolerance:          1.e-6
ch->space->solver->print cycle:        10
ch->space->solver->info:               10

parallel->debug->print options info: 1

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

ch->phasefield->output->filename: ./output/phasefield_dd_
ch->phasefield->output->ParaView animation:        1
ch->phasefield->output->ParaView format:           1
ch->phasefield->output->write every i-th timestep: 1
ch->phasefield->output->append index:              1
ch->phasefield->output->index length:              9
ch->phasefield->output->index decimals:            7