SFEMaNS
version 4.1 (work in progress)
Reference documentation for SFEMaNS

It is used to initialize the velocity field, the pressure and the pressure increment.
The inputs of this function are the following:
mesh_f
is the finite element mesh used to approximate the velocity field. mesh_c
is the finite element mesh used to approximate the pressure. dt
is the time step. list_mode
is a list of integers which contains the Fourier modes approximated. As the mesh can be subdivised, we note that mesh_f and mesh_c depend of the processor considered when doing parallel computing. Same goes for the list of Fourier mode list_mode
.
The outputs of this function are the following:
time
is the time when the computations starts. un_m1
is the velocity field at the time timedt
. un
is the velocity field at the time time
. pn_m1
is the pressure at the time timedt
. pn
is the pressure at the time time
. phin_m1
is the increment pressure at the time timedt
. phin
is the increment pressure at the time time
. Remarks:
mesh_f
. mesh_p
. Here is an exemple where we use the function vv_exact and pp_exact to initialize the velocity field and the pressire. Thus, the initial conditions satisfy the boundary conditions.
We note that the integers mode, i and j have to be declared. Same for the real valued tabular pn_m2 whose dimension is mesh_c%np. It is done by adding the two following lines in the declaration of the function init_velocity_pressure
:
We refer to the sections Examples with manufactured solutions and Examples on physical problems for more examples.