pdd_plast_flow_mixed/include/pdh_plast_flow_mixed_problem.h

Go to the documentation of this file.

/************************************************************************
File pdh_plast_flow_problem.h - problem module's type

Contains problem module defines (see below)

Contains definition of types:
  pdt_plast_flow_ctrls
  pdt_plast_flow_times
  pdt_plast_flow_nonls
  pdt_plast_flow_linss
  pdt_plast_flow_adpts
  pdt_plast_flow_problem - aggregates above ones
  
Procedures:
  pdr_plast_flow_problem_clear - clear problem data
  pdr_plast_flow_problem_read - read problem data

------------------------------
History:
    initial version - Krzysztof Banas (pobanas@cyf-kr.edu.pl)
*************************************************************************/

#ifndef PDH_PLAST_FLOW_PROBLEM
#define PDH_PLAST_FLOW_PROBLEM

#ifndef ANALYTIC_EXPRESSIONS
#define ANALYTIC_EXPRESSIONS
#endif

#ifndef CONVERGENCE_CURVE
#define CONVERGENCE_CURVE
#endif

#include <stdio.h>

/* problem dependent interface with the PDEs  */
#include "pdh_intf.h"

/* types and functions related to boundary conditions handling */
#include "pdh_plast_flow_mixed_bc.h" 

/* types and functions related to materials handling */
#include "pdh_plast_flow_mixed_materials.h"

#ifdef __cplusplus
extern "C" 
{
#endif

/**************************************/
/* DEFINES                            */
/**************************************/
/* Rules:
/* - always uppercase */
/* - name starts with PDC_ */

// maximal number of equations (solution components) in component modules
#define PDC_PLAST_FLOW_MAXEQ 4
#define PDC_PLAST_FLOW_NREQ 4
#define PDC_PLAST_FLOW_PRIMARY_FIELD_NREQ 3
#define PDC_PLAST_FLOW_SECONDARY_FIELD_NREQ 1
  
#define PDC_PLAST_FLOW_STRAIN_MAXEQ 30  // derived field of:    strain_intensity, flow_stress, strain_rate_int, 
#define PDC_PLAST_FLOW_STRAIN_NREQ 30   //                      9 strains, 9 strain rates, 9 stresses
#define PDC_PLAST_FLOW_STRAIN_MAXELSD APC_DOUBLE_MAXELSD // maximum vector size for dofs data from strain field

/**************************************/
/* TYPES                              */
/**************************************/
/* Rules:
/* - type name starts always witn pdt_ */

/* structure with control parameters */
typedef struct {
  /*** GENERIC - DO NOT CHANGE !!! ***/
  char    name[300];    /* name (identifier for problem dependent routines) */
  int     mesh_id;      /* ID of the associated mesh */
  int     primary_field_id;     /* ID of the associated PRIMARY approximation field */
  int     secondary_field_id; /* ID of the associated SECONDARY approximation field */
  int     nr_sol;       /* number of solution vectors stored by approximation */
                        /* module (for time dependent/nonlinear problems) */
  int     nreq[3];         /* number of equations (solution components) */
  int     solver_id;    /* ID of the associated solver */
// for continuous basis functions problem->ctrl.base is not used
  //int     base;         /* parameter specifying the type of basis functions */
                        /* interpreted by particular approximation modules */
  int     pdeg[4];         /* PDEG value for approximation field 
                              pdeg[0] - PDEG value
                              pdeg[1] - PDEG value of PRIMARY field
                              pdeg[2] - PDEG value of SECONDARY field
                              pdeg[3] - PDEG value of STRAIN FIELD field */
  char    mesh_type[2];
  
  char    work_dir[300];
  FILE*   interactive_input;
  FILE*   interactive_output;
  char    mesh_filename[300];
  char    field_filename[300];
  char    material_filename[300];
  char    bc_filename[300];
  char    solver_filename[300];
  char    field_dmp_filepattern[50];
  char    mesh_dmp_filepattern[50];
  /*** PROBLEM SPECIFIC - CAN BE MODIFIED ***/ 
  int     mesh_id_initial;      /* ID of the associated initial mesh */
  int     strain_field_id;     /* ID of the associated approximation field */
  int     strain_nr_sol;       /* number of solution vectors stored by approximation */
                        /* module (for time dependent/nonlinear problems) */
  int     strain_nreq;         /* number of equations (solution components) */
  char    strain_field_filename[300];
  char    strain_field_dmp_filepattern[50];
  double  penalty; // for boundary conditions
  double  density; // for constant density calculations
  double  ref_temperature; // for constant temperature calculations
  int     fields_nb;    // fields number
} pdt_plast_flow_ctrls;

/* structure with time integration parameters */
typedef struct {
  /*** GENERIC - DO NOT CHANGE !!! ***/
  int    type;          /* type of time integration scheme */
                        /*      0 - no time integration */
                        /*      1 - alpha (a.k.a. theta) scheme  */
  double alpha; /* implicitnes parameter alpha (a.k.a. theta)*/ 

  int    cur_step;      /* current time-step number */
  double cur_time;      /* current time */
  double cur_dtime;     /* current time-step length */
  double prev_dtime;    /* previous time-step length */

  int    final_step;    /* time-step number to stop simulation */
  double final_time;    /* time to stop simulation */
  
  int    conv_type;     /* convergence criterion number */
  double conv_meas;     /* convergence measure */
  int    monitor;       /* monitoring level: */
                        /*   PDC_SILENT      0  */
                        /*   PDC_ERRORS      1  */
                        /*   PDC_INFO        2  */
                        /*   PDC_ALLINFO     3  */
  int    intv_dumpout;  /* interval (in time steps) for dumping out data */
  int    intv_graph;    /* interval (in time steps) for graphics output */
  int    graph_accu;    /* auto graphics dumpout accuracy */

  int    time_step_length_nonl_control; // indicator (0/1) of time step adaptivity
                                // based on the convergence of non-linear solution
  int    time_step_length_nonl_iter_max; // the number of non-linear iterations above which
                         //  time step length is decreased (0 for non-decreasing time-step)
  int    time_step_length_nonl_iter_min; // the number of non-linear iterations below which
                         //  time step length is increased (0 for non-increasing time-step)
  double time_step_length_nonl_iter_increase_mult;
  double time_step_length_nonl_iter_decrease_mult;

  /*** PROBLEM SPECIFIC - CAN BE MODIFIED ***/
} pdt_plast_flow_times;

/* structure with nonlinear solver control parameters */
typedef struct {
  /*** GENERIC - DO NOT CHANGE !!! ***/
  int    type;          /* method identifier */
                        /*      0 - problem is linear */
  int    max_iter;      /* maximal iteration number */
  int    conv_type;     /* convergence criterion number */
  double conv_meas;     /* convergence measure */
  int    monitor;       /* monitoring level: */
                        /*   PDC_SILENT      0  */
                        /*   PDC_ERRORS      1  */
                        /*   PDC_INFO        2  */
                        /*   PDC_ALLINFO     3  */
  /*** PROBLEM SPECIFIC - CAN BE MODIFIED ***/
} pdt_plast_flow_nonls;

/* structure with linear solver control parameters */
typedef struct {
  /*** GENERIC - DO NOT CHANGE !!! ***/
  int    type;          /* method identifier */
                        /*       0 - direct solver (?) */
                        /*       1 - precondittioned GMRES */
                        /*       2 - multigrid preconditioned  GMRES */
                        /*      10 - standard iterations */
                        /*      20 - V-cycle multigrid */
  int    max_iter;      /* maximal iteration number */
  int    conv_type;     /* convergence criterion number */
                        /*      0 - relative to initial residual */
                        /*      1 - absolute residual */
                        /*      2 - relative to rhs */
  double conv_meas;     /* convergence measure */
  int    monitor;       /* monitoring level: */
                        /*   PDC_SILENT      0  */
                        /*   PDC_ERRORS      1  */
                        /*   PDC_INFO        2  */
                        /*   PDC_ALLINFO     3  */
  /*** PROBLEM SPECIFIC - CAN BE MODIFIED ***/
} pdt_plast_flow_linss;

/* structure with adaptation parameters */
typedef struct {
  /*** GENERIC - DO NOT CHANGE !!! ***/
  int    type;          /* strategy number for adaptation */
  int    interval;      /* number of time steps between adaptations */
  int    maxgen;        /* maximum generation level for elements */
  double eps;           /* coefficient for choosing elements to adapt */
  double ratio;         /* ratio of errors for derefinements */
  int    monitor;       /* monitoring level: */
                        /*   PDC_SILENT      0  */
                        /*   PDC_ERRORS      1  */
                        /*   PDC_INFO        2  */
                        /*   PDC_ALLINFO     3  */
  /*** PROBLEM SPECIFIC - CAN BE MODIFIED ***/
} pdt_plast_flow_adpts;

/* problem definition data structure */
typedef struct {
    pdt_plast_flow_ctrls  ctrl;    /* structure with control parameters */
    pdt_plast_flow_times  time;   /* structure with time integration parameters */
    pdt_plast_flow_nonls  nonl;   /* structure with nonlinear solver parameters */
    pdt_plast_flow_linss  lins;   /* structure with linear solver parameters */
    pdt_plast_flow_adpts  adpt;    /* structure with adaptation parameters */
    //pdt_plast_flow_materials materials; /* structure containing materials data */
    pdt_plast_flow_bc bc; /* structure containing bc data */
} pdt_plast_flow_problem;


extern  int pdr_plast_flow_problem_clear(pdt_plast_flow_problem *Problem);

extern int pdr_plast_flow_problem_read(
  char *Work_dir,
  char *Filename,
  FILE *Interactive_output,
  pdt_plast_flow_problem *Problem,
  int Nr_sol // nr_sol is time integration dependent - specified by main
                             );
  
#ifdef __cplusplus
}
#endif

#endif

---