test_nona.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This file has been generated by NS2DDV routine generate_setup_file.m %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This part has been automatically generated: please do not modify it !
% The parameters you can modify are located at the end of this file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
PARAMETERS.PARALLELIZATION.SERVER = 'LOCAL';
PARAMETERS.MODEL = 'NS';
PARAMETERS.TESTCASE = 'NONA';
PARAMETERS.CV_STUDY = 'NONE';
PARAMETERS.DOMAIN.GEOMETRY = 'DIHEDRON';
PARAMETERS.PARALLELIZATION.TOOLBOX = 'SERIAL';
PARAMETERS.LINSYS.DIRECT_INVERSION = 'YES';
PARAMETERS.MESH.ADAPTATIVE = 'NO';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% From this line, you can modify the values of the parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Domain parameters
%%%%%%%%%%%%%%%%%%%
% Parameters for a DIHEDRON domain
% Bounding rectangle (default = [-1,5]x[0,2])
PARAMETERS.DOMAIN.XMIN = -1;
PARAMETERS.DOMAIN.XMAX = 7;
PARAMETERS.DOMAIN.YMIN = 0;
PARAMETERS.DOMAIN.YMAX = 3;
% Coordinates of the arete point (default = (0,1))
PARAMETERS.DOMAIN.XAXIS_ANGLE1 = 0;
PARAMETERS.DOMAIN.YAXIS_ANGLE1 = 1;
% x-coordinate of the corner point (default = 3)
PARAMETERS.DOMAIN.XAXIS_ANGLE2 = 3;
% Physical parameters
%%%%%%%%%%%%%%%%%%%%%
% Viscosity kind
% 'CONST' Constant viscosity
PARAMETERS.PHYSICAL.VISCOSITY = 'CONST';
% Final time
PARAMETERS.FINAL_TIME = 1.;
% Reynolds number
PARAMETERS.PHYSICAL.RE = 100.;
% Gravity amplitude (Default = 1)
PARAMETERS.PHYSICAL.GRAVITY = 1.;
%% x-coordinate of the boundary layer origin
PARAMETERS.PHYSICAL.XSTART = -2;
% Signals points
% Warning: this list must have the format {[x1,y1], [x2,y2], ...}
PARAMETERS.PHYSICAL.SIGNALS = {[3.,0.1], [3.,1.], [3.,2.], [4.,0.1], [4.,1.], [4.,2.], [4.5,0.1], [4.5,1.], [4.5,2.], [4.9,0.1], [4.9,1.], [4.9,2.]};
%% x-limits of the injected stream zone
%% Format must be {ximin,ximax} with xA1 < ximin < ximax < xA2, or {} if there is no injection
PARAMETERS.PHYSICAL.INJECTION_XLIM = {};
%% Type of injection stream (allowed values: 'CONSTANT')
PARAMETERS.PHYSICAL.INJECTION_TYPE = 'CONSTANT';
%% Velocity magnitude of the injected stream
PARAMETERS.PHYSICAL.INJECTION_UMAGNI = 0.2;
%% Angle (in radians) between the injected stream velocity and the vector (1,0) in the clockwise order
PARAMETERS.PHYSICAL.INJECTION_UANGLE = 0.;
% Initialization with external data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Use a third-party file to initialize the simulation
% Allowed formats: '', 'xxx.h5', 'xxx.mat'
% If '' (default value) is selected, the simulation will be initialized with the test case parameters
% If a valid file is selected, some domain and physical parameters above will be ignored.
PARAMETERS.INIT_FILE = '';
% Build a new mesh at initialization (ignored if PARAMETERS.INIT_FILE = '')
% Allowed values: 'YES' or 'NO' (default)
PARAMETERS.REMESH = 'NO';
% Use domain parameters from the third-party initialization file or PARAMETERS.DOMAIN.XXX
% (ignored if PARAMETERS.INIT_FILE = '' or PARAMETERS.REMESH = 'NO')
% 'EXT_FILE' The domain geometry parameters are imposed by the contents of PARAMETERS.INIT_FILE
% 'CURRENT_FILE' Use the variables PARAMETERS.DOMAIN.XXX below
PARAMETERS.WHICH_GEOMETRY = 'EXT_FILE';
% Reset time variable if a third-party file is used to initialize the simulation (ignored if PARAMETERS.INIT_FILE = '')
% Allowed values: 'YES' (default) or 'NO'
PARAMETERS.RESET_TIME = 'YES';
% Mesh parameters
%%%%%%%%%%%%%%%%%
% Method for generating the space mesh
% 'PDET' Use the Matlab PDE toolbox
% 'NS2DDV' Use the routines embedded in NS2DDV code
% 'FROM_FILE' Load data from a mesh file (.msh or .p1/.t1 format allowed)
PARAMETERS.MESH.GENERATION = 'NS2DDV';
% Node renumbering method
% 'NONE' No renumbering
% 'AMD' Renumbering with Minimum Degree Algorithm Matlab built-in routines
% 'CMK' Renumbering with Cuthill-McKee Matlab built-in routines
PARAMETERS.MESH.RENUMBERING = 'NONE';
% Parameters for a NS2DDV mesh
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Height of the boundary layer
PARAMETERS.MESH.HEIGHT_BL = 0.2;
% Number of edges by segment of length 1 in x-direction (even value recommended)
PARAMETERS.MESH.NBSEG_PERUNIT_X = 6;
% Number of nodes in the boundary layer in y-direction
PARAMETERS.MESH.NBNODES_INBL_Y = 2;
% Geometric reason of nodes repartition in the boundary layer in y-direction
PARAMETERS.MESH.PROG_GEOM_BL_Y = 1.1;
% Number of edges out of the boundary layer in y-direction
PARAMETERS.MESH.NBSEG_OUTBL_Y = 9;
% Design of structured mesh ('DIAGONAL' or 'CROSS')
PARAMETERS.MESH.TRIANGLES_ORIENTATION = 'CROSS';
% Parameters for a PDET mesh
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Characteristic edge length in the mesh
PARAMETERS.MESH.H0 = 0.2;
% Parameters for a FROM_FILE mesh
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Mesh source file(s) : .msh or .p1/.t1 format allowed
% {'xxxx.msh'} cell array of size 1 for .msh format
% {'xxxx.p1', 'xxxx.t1'} cell array of size 2 for .p1/.t1 format
PARAMETERS.MESH.SRC_FILES = {};
% Rescale the mesh to fit with the domain geometry above ('YES' or 'NO')
PARAMETERS.MESH.RESCALE = 'NO';
% Boundary condition parameters for velocity
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 'DIRICHLET' Use Dirichlet boundary condition on the specified part of the boundary
% 'NEUMANN' Use Neumann boundary condition on the specified part of the boundary
% 'NATURAL' Use transparent boundary condition on the specified part of the boundary
% Boundary condition on the bottom bound for ux (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_BOTTOM_UX = 'DIRICHLET';
% Boundary condition on the bottom bound for uy (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_BOTTOM_UY = 'DIRICHLET';
% Boundary condition on the top bound for ux (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_TOP_UX = 'DIRICHLET';
% Boundary condition on the top bound for uy (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_TOP_UY = 'DIRICHLET';
% Boundary condition on the left bound for ux (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_LEFT_UX = 'DIRICHLET';
% Boundary condition on the left bound for uy (allowed values: 'DIRICHLET')
PARAMETERS.FE.BC_LEFT_UY = 'DIRICHLET';
% Boundary condition on the right bound for ux (allowed values: 'NEUMANN', 'NATURAL', 'STRESS_BALANCE')
PARAMETERS.FE.BC_RIGHT_UX = 'NEUMANN';
% Boundary condition on the right bound for uy (allowed values: 'NEUMANN', 'NATURAL', 'STRESS_BALANCE')
PARAMETERS.FE.BC_RIGHT_UY = 'NEUMANN';
% Non-linear term splitting parameters for 'NATURAL' boundary conditions
PARAMETERS.FE.ALPHA1 = 0.5;
PARAMETERS.FE.ALPHA2 = 1.;
PARAMETERS.FE.ALPHA3 = 1.;
% Amplitude of the correction in 'NATURAL' boundary conditions (recommended values: 0, 1, or 2; must be positive)
PARAMETERS.FE.BETA = 1.;
% Extension of Dirichlet boundary conditions to the open boundary
% Allowed values: 'BLASIUS' (default), 'COUETTE'
% 'BLASIUS' Use a barotropic Blasius profile
% 'COUETTE' Use a barotropic Couette profile
PARAMETERS.FE.BC_EXTENSIONREF = 'BLASIUS';
% Sharpness of the smoothed step function - should be close to 0 (only for 'STRESS_BALANCE' boundary conditions)
PARAMETERS.FE.BC_DELTA = 0.05;
% Finite Element parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Time semi-discretization of the Stokes equation
% 'BDF2' Implicit 2-steps Backward Differentiation Formula
% 'BDF2_PROJ' Implicit 2-steps Backward Differentiation Formula + Projection method
% Warning: 'STRESS_BALANCE' and 'NATURAL' boundary conditions are compatible only with 'BDF2'
PARAMETERS.FE.SCHEME = 'BDF2';
% Amplitude of the rotational term in the pressure correction for 'BDF2_PROJ*' methods
PARAMETERS.FE.ROT_CORRECTION = 0.5/PARAMETERS.PHYSICAL.RE;
% Space discretization of velocity
% 'P1B' P1-bubble finite elements approximation
% 'P2' P2 finite elements approximation
PARAMETERS.FE.TYPE = 'P2';
% Space discretization of density in Stokes equation
% 'P1' P1 finite elements approximation
PARAMETERS.FE.TYPE_DENSITY = 'P1';
% Constrain for closing the pressure problem
% 'ZERO_FIXED_POINT' p = 0 is imposed on a specific point in the space domain
% 'ZERO_AVERAGE' the average of p is set to 0 within the discrete system
% 'ZERO_AVERAGE_APOSTERIORI' the average of p is set to 0 as an a posteriori constrain (only works with PARAMETERS.FE.SCHEME = 'BDF2')
PARAMETERS.FE.PRESSURE_CONSTRAIN = 'ZERO_FIXED_POINT';
% Coordinates of the null pressure point
% WARNING: if the domain geometry parameters are obtained from an external file, it is not recommended to use
% the variables PARAMETERS.DOMAIN.XXX to define the coordinates of the null pressure point. Use numerical values instead
PARAMETERS.FE.NFX_X = 0.1*PARAMETERS.DOMAIN.XMIN+0.9*PARAMETERS.DOMAIN.XMAX;
PARAMETERS.FE.NFX_Y = 0.1*PARAMETERS.DOMAIN.YMIN+0.9*PARAMETERS.DOMAIN.YMAX;
% Time step for solving Stokes equation
% WARNING : this time step is of the form Cmax*hmax^alphamax + Cmin*hmin^alphamin
% where hmax (hmin) is the max (min) length of an edge
% Value of alphamax
PARAMETERS.FE.ALPHAMAX_STEP_TIME = 1;
% Value of alphamin
PARAMETERS.FE.ALPHAMIN_STEP_TIME = 1;
% Value of Cmax
PARAMETERS.FE.CMAX_STEP_TIME = 0;
% Value of Cmin
PARAMETERS.FE.CMIN_STEP_TIME = 1;
% Reduce the time step to reach the final time T exactly ('YES' or 'NO')
% Warning: answering 'NO' will probably make the simulation stop a little before the final time
PARAMETERS.FE.FIT_STEP_TIME = 'YES';
% Visualization parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%
% Display 2D results in a Matlab graphical window :
% 'SINGLE_FIGURE' All graphical outputs in a unique figure
% 'MULTIPLE_FIGURE' One graphical output per figure, stable
% 'NONE' No graphical output
PARAMETERS.OUTPUT.XDISPLAY = 'SINGLE_FIGURE';
% Refreshing frequency for 2D graphical display (default = 1)
PARAMETERS.OUTPUT.XDISPLAY_FREQUENCY = 1;
% List of 2D results to be plotted :
% 'MESH' The P1 mesh
% 'SUBTRI_U' The subtriangulation for velocity (either P2 or P1B)
% 'U_VECTOR' The velocity vector field
% 'UX' The isovalues of the x-component of velocity
% 'UY' The isovalues of the y-component of velocity
% 'UX_DX' The isovalues of the x-derivative of the x-component of velocity
% 'UX_DY' The isovalues of the y-derivative of the x-component of velocity
% 'UY_DX' The isovalues of the x-derivative of the y-component of velocity
% 'UY_DY' The isovalues of the y-derivative of the y-component of velocity
% 'SHEAR' The isovalues of the shear rate
% 'U_MODULE' The velocity magnitude isovalues
% 'VORTIC' The vorticity isovalues
% 'STREAM' The velocity streamlines
% 'P' The pressure isovalues
% 'GRAD_P' The pressure gradient
% 'P_DX' The isovalues of the x-derivative of pressure
% 'P_DY' The isovalues of the y-derivative of pressure
% 'W_VECTOR' The auxiliary Stokes velocity vector field (only for NATURAL boundary conditions)
% 'WX' The isovalues of the x-component of the auxiliary Stokes velocity (only for NATURAL boundary conditions)
% 'WY' The isovalues of the y-component of the auxiliary Stokes velocity (only for NATURAL boundary conditions)
% 'R' The auxiliary Stokes pressure isovalues (only for NATURAL boundary conditions)
PARAMETERS.OUTPUT.XDISPLAY_LIST = {'MESH', 'VORTIC', 'STREAM', 'U_MODULE', 'P_DX', 'P_DY'};
% Output parameters
%%%%%%%%%%%%%%%%%%%
% Directory for numerical results
PARAMETERS.OUTPUT.DIRECTORY_NAME = './RESULTS/EXAMPLES/test_nona';
% Name of output file(s)
PARAMETERS.OUTPUT.FILE_NAME = 'diags';
% Save method of these results :
% 'LAST_FRAME' The solution is only saved at final time
% 'ANIMATION' If time dynamics is considered, an animated result is saved
% 'NONE' Only the log file is saved
PARAMETERS.OUTPUT.XDISPLAY_SAVE = 'ANIMATION';
% Refreshing frequency for save 2D results (default = 1)
PARAMETERS.OUTPUT.XDISPLAY_SAVE_FREQUENCY = 1;
% List of optional diagnostics
% 'UX_DX' x-derivative of x-component of velocity
% 'UX_DY' y-derivative of x-component of velocity
% 'UY_DX' x-derivative of y-component of velocity
% 'UY_DY' y-derivative of y-component of velocity
% 'SHEAR' shear rate
% 'GRAD_P' Pressure gradient
% 'P_DX' x-derivative of the pressure
% 'P_DY' y-derivative of the pressure
% 'W_VECTOR' Auxiliary Stokes velocity (only for NATURAL boundary conditions)
% 'R' Auxiliary Stokes pressure (only for NATURAL boundary conditions)
PARAMETERS.OUTPUT.XDISPLAY_SAVE_OPTDIAGS = {'GRAD_P', 'UX_DX'};
% Backup parameters
%%%%%%%%%%%%%%%%%%%
% Number of backup files (default = 1)
% Each backup file contains the required data to restart the code after a critical crash
PARAMETERS.BACKUP.NB_FILES = 1;
% Backup frequency (default = 10)
PARAMETERS.BACKUP.FREQUENCY = 10;
% Backup file location
PARAMETERS.BACKUP.DIRECTORY_NAME = './BACKUP/EXAMPLES/test_nona';