.. Created using 4C version (git SHA1): 94d307e8b81d05a5cec0d1fd47ca58f4be040deb .. _headerparameters: Header parameters ================= .. _SECaledynamic: ALE DYNAMIC ----------- ALE DYNAMIC :: // ALE DYNAMIC: // TIMESTEP (default: 0.1) "time step size" // NUMSTEP (default: 41) "max number of time steps" // MAXTIME (default: 4) "max simulation time" // ALE_TYPE (default: solid) "ale mesh movement algorithm" // ASSESSMESHQUALITY (default: false) "Evaluate element quality measure according to [Oddy et al. 1988]" // UPDATEMATRIX (default: false) "Update stiffness matrix in every time step (only for linear/material strategies)" // MAXITER (default: 1) "Maximum number of newton iterations." // TOLRES (default: 1e-06) "Absolute tolerance for length scaled L2 residual norm" // TOLDISP (default: 1e-06) "Absolute tolerance for length scaled L2 increment norm" // NUM_INITSTEP (default: 0) // RESTARTEVERY (default: 1) "write restart data every RESTARTEVERY steps" // RESULTSEVERY (default: 0) "write results every RESULTSTEVERY steps" // DIVERCONT (default: continue) (choices: continue|stop) "What to do if nonlinear solver does not converge?" // MESHTYING (default: no) (choices: meshsliding|meshtying|no) "Flag to (de)activate mesh tying and mesh sliding algorithm" // INITIALDISP (default: zero_displacement) (choices: displacement_by_function|zero_displacement) "Initial displacement for structure problem" // STARTFUNCNO (default: -1) "Function for Initial displacement" // LINEAR_SOLVER (default: -1) "number of linear solver used for ale problems..." .. _SECalegeometry: ALE GEOMETRY ------------ ALE GEOMETRY :: // ALE GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECarterialdynamic: ARTERIAL DYNAMIC ---------------- ARTERIAL DYNAMIC :: // ARTERIAL DYNAMIC: // DYNAMICTYPE (default: ExpTaylorGalerkin) (choices: ExpTaylorGalerkin|Stationary) "Explicit Taylor Galerkin Scheme" // TIMESTEP (default: 0.01) "Time increment dt" // NUMSTEP (default: 0) "Number of Time Steps" // MAXTIME (default: 1000) "total simulation time" // RESTARTEVERY (default: 1) "Increment for writing restart" // RESULTSEVERY (default: 1) "Increment for writing solution" // SOLVESCATRA (default: false) "Flag to (de)activate solving scalar transport in blood" // LINEAR_SOLVER (default: -1) "number of linear solver used for arterial dynamics" // INITFUNCNO (default: -1) "function number for artery initial field" // INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "Initial Field for artery problem" .. _SECarterygeometry: ARTERY GEOMETRY --------------- ARTERY GEOMETRY :: // ARTERY GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECbeamcontact: BEAM CONTACT ------------ BEAM CONTACT :: // BEAM CONTACT: // BEAMS_STRATEGY (default: None) (choices: None|Penalty|none|penalty) "Type of employed solving strategy" // MODELEVALUATOR (default: Old) (choices: Old|Standard|old|standard) "Type of model evaluator" // BEAMS_NEWGAP (default: false) "choose between original or enhanced gapfunction" // BEAMS_SEGCON (default: false) "choose between beam contact with and without subsegment generation" // BEAMS_DEBUG (default: false) "This flag can be used for testing purposes. When it is switched on, some sanity checks are not performed!" // BEAMS_INACTIVESTIFF (default: false) "Always apply contact stiffness in first Newton step for pairs which have active in last time step" // BEAMS_BTSOL (default: false) "decide, if also the contact between beams and solids is possible" // BEAMS_ENDPOINTPENALTY (default: false) "Additional consideration of endpoint-line and endpoint-endpoint contacts" // BEAMS_SMOOTHING (default: None) (choices: Cpp|None|cpp|none) "Application of smoothed tangent field" // BEAMS_DAMPING (default: false) "Application of a contact damping force" // BEAMS_BTBPENALTYPARAM (default: 0) "Penalty parameter for beam-to-beam point contact" // BEAMS_BTBLINEPENALTYPARAM (default: -1) "Penalty parameter per unit length for beam-to-beam line contact" // BEAMS_BTSPENALTYPARAM (default: 0) "Penalty parameter for beam-to-solid contact" // BEAMS_DAMPINGPARAM (default: -1000) "Damping parameter for contact damping force" // BEAMS_DAMPREGPARAM1 (default: -1000) "First (at gap1, with gap1>gap2) regularization parameter for contact damping force" // BEAMS_DAMPREGPARAM2 (default: -1000) "Second (at gap2, with gap1>gap2) regularization parameter for contact damping force" // BEAMS_MAXDISISCALEFAC (default: -1) "Scale factor in order to limit maximal iterative displacement increment (resiudal displacement)" // BEAMS_MAXDELTADISSCALEFAC (default: 1) "Scale factor in order to limit maximal displacement per time step" // BEAMS_PERPSHIFTANGLE1 (default: -1) "Lower shift angle (in degrees) for penalty scaling of large-angle-contact" // BEAMS_PERPSHIFTANGLE2 (default: -1) "Upper shift angle (in degrees) for penalty scaling of large-angle-contact" // BEAMS_PARSHIFTANGLE1 (default: -1) "Lower shift angle (in degrees) for penalty scaling of small-angle-contact" // BEAMS_PARSHIFTANGLE2 (default: -1) "Upper shift angle (in degrees) for penalty scaling of small-angle-contact" // BEAMS_SEGANGLE (default: -1) "Maximal angle deviation allowed for contact search segmentation" // BEAMS_NUMINTEGRATIONINTERVAL (default: 1) "Number of integration intervals per element" // BEAMS_PENALTYLAW (default: LinPen) (choices: LinNegQuadPen|LinPen|LinPosCubPen|LinPosDoubleQuadPen|LinPosExpPen|LinPosQuadPen|QuadPen) "Applied Penalty Law" // BEAMS_PENREGPARAM_G0 (default: -1) "First penalty regularization parameter G0 >=0: For gap (default: -1) "Second penalty regularization parameter F0 >=0: F0 represents the force at the transition point between regularized and linear force law!" // BEAMS_PENREGPARAM_C0 (default: -1) "Third penalty regularization parameter C0 >=0: C0 has different physical meanings for the different penalty laws!" // BEAMS_GAPSHIFTPARAM (default: 0) "Parameter to shift penalty law!" // BEAMS_BASICSTIFFGAP (default: -1) "For gaps > -BEAMS_BASICSTIFFGAP, only the basic part of the contact linearization is applied!" // BEAMS_OCTREE (default: None) (choices: None|none|octree_axisaligned|octree_cylorient|octree_spherical) "octree and bounding box type for octree search routine" // BEAMS_ADDITEXT (default: true) "Switch between No==multiplicative extrusion factor and Yes==additive extrusion factor" // BEAMS_EXTVAL (default: -1.0) "extrusion value(s) of the bounding box, Depending on BEAMS_ADDITIVEEXTFAC is either additive or multiplicative. Give one or two values." // BEAMS_TREEDEPTH (default: 6) "max. tree depth of the octree" // BEAMS_BOXESINOCT (default: 8) "max number of bounding boxes in any leaf octant" .. _SECbeamcontact_runtimevtkoutput: BEAM CONTACT/RUNTIME VTK OUTPUT ------------------------------- BEAM CONTACT/RUNTIME VTK OUTPUT :: // BEAM CONTACT/RUNTIME VTK OUTPUT: // VTK_OUTPUT_BEAM_CONTACT (default: false) "write visualization output for beam contact" // INTERVAL_STEPS (default: -1) "write visualization output at runtime every INTERVAL_STEPS steps" // EVERY_ITERATION (default: false) "write output in every iteration of the nonlinear solver" // CONTACT_FORCES (default: false) "write visualization output for contact forces" // GAPS (default: false) "write visualization output for gap, i.e. penetration" .. _SECbeaminteraction: BEAM INTERACTION ---------------- BEAM INTERACTION :: // BEAM INTERACTION: // REPARTITIONSTRATEGY (default: Adaptive) (choices: Adaptive|Everydt|adaptive|everydt) "Type of employed repartitioning strategy" // SEARCH_STRATEGY (default: bruteforce_with_binning) "Type of search strategy used for finding coupling pairs" .. _SECbeaminteraction_beamtobeamcontact: BEAM INTERACTION/BEAM TO BEAM CONTACT ------------------------------------- BEAM INTERACTION/BEAM TO BEAM CONTACT :: // BEAM INTERACTION/BEAM TO BEAM CONTACT: // STRATEGY (default: None) (choices: None|Penalty|none|penalty) "Type of employed solving strategy" .. _SECbeaminteraction_beamtosolidsurfacecontact: BEAM INTERACTION/BEAM TO SOLID SURFACE CONTACT ---------------------------------------------- BEAM INTERACTION/BEAM TO SOLID SURFACE CONTACT :: // BEAM INTERACTION/BEAM TO SOLID SURFACE CONTACT: // CONTACT_DISCRETIZATION (default: none) "Type of employed contact discretization" // CONSTRAINT_STRATEGY (default: none) "Type of employed constraint enforcement strategy" // PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-solid surface contact" // CONTACT_TYPE (default: none) "How the contact constraints are formulated" // PENALTY_LAW (default: none) "Type of penalty law" // PENALTY_PARAMETER_G0 (default: 0) "First penalty regularization parameter G0 >=0: For gap (default: none) "Configuration where the mortar contact is defined" // MORTAR_SHAPE_FUNCTION (default: none) "Shape function for the mortar Lagrange-multipliers" // GEOMETRY_PAIR_STRATEGY (default: segmentation) "Type of employed segmentation strategy" // GEOMETRY_PAIR_SEGMENTATION_SEARCH_POINTS (default: 6) "Number of search points for segmentation" // GEOMETRY_PAIR_SEGMENTATION_NOT_ALL_GAUSS_POINTS_PROJECT_VALID_ACTION (default: fail) "What to do if not all Gauss points of a segment project valid" // GAUSS_POINTS (default: 6) "Number of Gauss Points for the integral evaluations" // INTEGRATION_POINTS_CIRCUMFERENCE (default: 6) "Number of Integration points along the circumferential direction of the beam. This is parameter is only used in beam to cylinder meshtying. No gauss integration is used along the circumferential direction, equally spaced integration points are used." // GEOMETRY_PAIR_SURFACE_NORMALS (default: standard) "How the surface normals are evaluated" .. _SECbeaminteraction_beamtosolidsurfacemeshtying: BEAM INTERACTION/BEAM TO SOLID SURFACE MESHTYING ------------------------------------------------ BEAM INTERACTION/BEAM TO SOLID SURFACE MESHTYING :: // BEAM INTERACTION/BEAM TO SOLID SURFACE MESHTYING: // CONTACT_DISCRETIZATION (default: none) "Type of employed contact discretization" // CONSTRAINT_STRATEGY (default: none) "Type of employed constraint enforcement strategy" // COUPLING_TYPE (default: none) "How the coupling constraints are formulated/" // MORTAR_SHAPE_FUNCTION (default: none) "Shape function for the mortar Lagrange-multipliers" // PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-solid surface meshtying" // ROTATIONAL_COUPLING (default: false) "Enable / disable rotational coupling" // ROTATIONAL_COUPLING_PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-solid surface rotational meshtying" // ROTATIONAL_COUPLING_SURFACE_TRIAD (default: none) "Construction method for surface triad" // GEOMETRY_PAIR_STRATEGY (default: segmentation) "Type of employed segmentation strategy" // GEOMETRY_PAIR_SEGMENTATION_SEARCH_POINTS (default: 6) "Number of search points for segmentation" // GEOMETRY_PAIR_SEGMENTATION_NOT_ALL_GAUSS_POINTS_PROJECT_VALID_ACTION (default: fail) "What to do if not all Gauss points of a segment project valid" // GAUSS_POINTS (default: 6) "Number of Gauss Points for the integral evaluations" // INTEGRATION_POINTS_CIRCUMFERENCE (default: 6) "Number of Integration points along the circumferential direction of the beam. This is parameter is only used in beam to cylinder meshtying. No gauss integration is used along the circumferential direction, equally spaced integration points are used." // GEOMETRY_PAIR_SURFACE_NORMALS (default: standard) "How the surface normals are evaluated" .. _SECbeaminteraction_beamtosolidsurface_runtimevtkoutput: BEAM INTERACTION/BEAM TO SOLID SURFACE/RUNTIME VTK OUTPUT --------------------------------------------------------- BEAM INTERACTION/BEAM TO SOLID SURFACE/RUNTIME VTK OUTPUT :: // BEAM INTERACTION/BEAM TO SOLID SURFACE/RUNTIME VTK OUTPUT: // WRITE_OUTPUT (default: false) "Enable / disable beam-to-solid volume mesh tying output." // NODAL_FORCES (default: false) "Enable / disable output of the resulting nodal forces due to beam to solid interaction." // AVERAGED_NORMALS (default: false) "Enable / disable output of averaged nodal normals on the surface." // MORTAR_LAMBDA_DISCRET (default: false) "Enable / disable output of the discrete Lagrange multipliers at the node of the Lagrange multiplier shape functions." // MORTAR_LAMBDA_CONTINUOUS (default: false) "Enable / disable output of the continuous Lagrange multipliers function along the beam." // MORTAR_LAMBDA_CONTINUOUS_SEGMENTS (default: 5) "Number of segments for continuous mortar output" // SEGMENTATION (default: false) "Enable / disable output of segmentation points." // INTEGRATION_POINTS (default: false) "Enable / disable output of used integration points. If the contact method has 'forces' at the integration point, they will also be output." // UNIQUE_IDS (default: false) "Enable / disable output of unique IDs (mainly for testing of created VTK files)." .. _SECbeaminteraction_beamtosolidvolumemeshtying: BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING ----------------------------------------------- BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING :: // BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING: // CONTACT_DISCRETIZATION (default: none) "Type of employed contact discretization" // CONSTRAINT_STRATEGY (default: none) "Type of employed constraint enforcement strategy" // MORTAR_SHAPE_FUNCTION (default: none) "Shape function for the mortar Lagrange-multipliers" // MORTAR_FOURIER_MODES (default: -1) "Number of fourier modes to be used for cross-section mortar coupling" // PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-solid volume meshtying" // COUPLE_RESTART_STATE (default: false) "Enable / disable the coupling of the restart configuration." // ROTATION_COUPLING (default: none) "Type of rotational coupling" // ROTATION_COUPLING_MORTAR_SHAPE_FUNCTION (default: none) "Shape function for the mortar Lagrange-multipliers" // ROTATION_COUPLING_PENALTY_PARAMETER (default: 0) "Penalty parameter for rotational coupling in beam-to-solid volume mesh tying" // GEOMETRY_PAIR_STRATEGY (default: segmentation) "Type of employed segmentation strategy" // GEOMETRY_PAIR_SEGMENTATION_SEARCH_POINTS (default: 6) "Number of search points for segmentation" // GEOMETRY_PAIR_SEGMENTATION_NOT_ALL_GAUSS_POINTS_PROJECT_VALID_ACTION (default: fail) "What to do if not all Gauss points of a segment project valid" // GAUSS_POINTS (default: 6) "Number of Gauss Points for the integral evaluations" // INTEGRATION_POINTS_CIRCUMFERENCE (default: 6) "Number of Integration points along the circumferential direction of the beam. This is parameter is only used in beam to cylinder meshtying. No gauss integration is used along the circumferential direction, equally spaced integration points are used." .. _SECbeaminteraction_beamtosolidvolumemeshtying_runtimevtkoutput: BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING/RUNTIME VTK OUTPUT ------------------------------------------------------------------ BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING/RUNTIME VTK OUTPUT :: // BEAM INTERACTION/BEAM TO SOLID VOLUME MESHTYING/RUNTIME VTK OUTPUT: // WRITE_OUTPUT (default: false) "Enable / disable beam-to-solid volume mesh tying output." // NODAL_FORCES (default: false) "Enable / disable output of the resulting nodal forces due to beam to solid interaction." // MORTAR_LAMBDA_DISCRET (default: false) "Enable / disable output of the discrete Lagrange multipliers at the node of the Lagrange multiplier shape functions." // MORTAR_LAMBDA_CONTINUOUS (default: false) "Enable / disable output of the continuous Lagrange multipliers function along the beam." // MORTAR_LAMBDA_CONTINUOUS_SEGMENTS (default: 5) "Number of segments for continuous mortar output" // MORTAR_LAMBDA_CONTINUOUS_SEGMENTS_CIRCUMFERENCE (default: 8) "Number of segments for continuous mortar output along the beam cross-section circumference" // SEGMENTATION (default: false) "Enable / disable output of segmentation points." // INTEGRATION_POINTS (default: false) "Enable / disable output of used integration points. If the contact method has 'forces' at the integration point, they will also be output." // UNIQUE_IDS (default: false) "Enable / disable output of unique IDs (mainly for testing of created VTK files)." .. _SECbeaminteraction_beamtospherecontact: BEAM INTERACTION/BEAM TO SPHERE CONTACT --------------------------------------- BEAM INTERACTION/BEAM TO SPHERE CONTACT :: // BEAM INTERACTION/BEAM TO SPHERE CONTACT: // STRATEGY (default: None) (choices: None|Penalty|none|penalty) "Type of employed solving strategy" // PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-rigidsphere contact" .. _SECbeaminteraction_crosslinking: BEAM INTERACTION/CROSSLINKING ----------------------------- BEAM INTERACTION/CROSSLINKING :: // BEAM INTERACTION/CROSSLINKING: // CROSSLINKER (default: false) "Crosslinker in problem" // INIT_LINKER_BOUNDINGBOX (default: 1e12 1e12 1e12 1e12 1e12 1e12) "Linker are initially set randomly within this bounding box" // TIMESTEP (default: -1) "time step for stochastic events concerning crosslinking (e.g. diffusion, p_link, p_unlink)" // VISCOSITY (default: 0) "viscosity" // KT (default: 0) "thermal energy" // MAXNUMINITCROSSLINKERPERTYPE (default: 0) "number of initial crosslinker of certain type (additional to NUMCROSSLINKERPERTYPE)" // NUMCROSSLINKERPERTYPE (default: 0) "number of crosslinker of certain type" // MATCROSSLINKERPERTYPE (default: -1) "material number characterizing crosslinker type" // MAXNUMBONDSPERFILAMENTBSPOT (default: 1) "maximal number of bonds per filament binding spot" // FILAMENTBSPOTINTERVALGLOBAL (default: -1.0) "distance between two binding spots on all filaments" // FILAMENTBSPOTINTERVALLOCAL (default: -1.0) "distance between two binding spots on current filament" // FILAMENTBSPOTRANGEGLOBAL (default: -1.0 -1.0) "Lower and upper arc parameter bound for binding spots on a filament" // FILAMENTBSPOTRANGELOCAL (default: 0.0 1.0) "Lower and upper arc parameter bound for binding spots on a filament" .. _SECbeaminteraction_spherebeamlink: BEAM INTERACTION/SPHERE BEAM LINK --------------------------------- BEAM INTERACTION/SPHERE BEAM LINK :: // BEAM INTERACTION/SPHERE BEAM LINK: // SPHEREBEAMLINKING (default: false) "Integrins in problem" // CONTRACTIONRATE (default: 0) "contraction rate of cell (integrin linker) in [microm/s]" // TIMESTEP (default: -1) "time step for stochastic events concerning sphere beam linking (e.g. catch-slip-bond behavior)" // MAXNUMLINKERPERTYPE (default: 0) "number of crosslinker of certain type" // MATLINKERPERTYPE (default: -1) "material number characterizing crosslinker type" // FILAMENTBSPOTINTERVALGLOBAL (default: -1.0) "distance between two binding spots on all filaments" // FILAMENTBSPOTINTERVALLOCAL (default: -1.0) "distance between two binding spots on current filament" // FILAMENTBSPOTRANGEGLOBAL (default: -1.0 -1.0) "Lower and upper arc parameter bound for binding spots on a filament" // FILAMENTBSPOTRANGELOCAL (default: 0.0 1.0) "Lower and upper arc parameter bound for binding spots on a filament" .. _SECbeampotential: BEAM POTENTIAL -------------- BEAM POTENTIAL :: // BEAM POTENTIAL: // POT_LAW_EXPONENT (default: 1.0) "negative(!) exponent(s) $m_i$ of potential law $\\Phi(r) = \\sum_i (k_i * r^{-m_i}).$" // POT_LAW_PREFACTOR (default: 0.0) "prefactor(s) $k_i$ of potential law $\\Phi(r) = \\sum_i (k_i * r^{-m_i})$." // TYPE (default: surface) "Type of potential interaction: surface (default) or volume potential" // STRATEGY (default: double_length_specific_large_separations) "strategy to evaluate interaction potential: double/single length specific, small/large separation approximation, ..." // CUTOFF_RADIUS > (default: none) "Neglect all potential contributions at separation largerthan this cutoff radius" // REGULARIZATION_TYPE (default: none) "Type of regularization applied to the force law" // REGULARIZATION_SEPARATION (default: -1) "Use regularization of force law at separations smaller than this separation" // N_INTEGRATION_SEGMENTS (default: 1) "Number of integration segments used per beam element" // N_GAUSS_POINTS (default: 10) "Number of Gauss points used per integration segment" // AUTOMATIC_DIFFERENTIATION (default: false) "apply automatic differentiation via FAD?" // CHOICE_MASTER_SLAVE (default: smaller_eleGID_is_slave) "According to which rule shall the role of master and slave be assigned to beam elements?" // POTENTIAL_REDUCTION_LENGTH > (default: none) "Within this length of the master beam end point the potential is smoothly reduced to one half to account for infinitely long master beam surrogates." .. _SECbeampotential_runtimevtkoutput: BEAM POTENTIAL/RUNTIME VTK OUTPUT --------------------------------- BEAM POTENTIAL/RUNTIME VTK OUTPUT :: // BEAM POTENTIAL/RUNTIME VTK OUTPUT: // VTK_OUTPUT_BEAM_POTENTIAL (default: false) "write visualization output for potential-based beam interactions" // INTERVAL_STEPS (default: 1) "write output at runtime every INTERVAL_STEPS steps" // EVERY_ITERATION (default: false) "write output in every iteration of the nonlinear solver" // FORCES (default: false) "write visualization output for forces" // MOMENTS (default: false) "write visualization output for moments" // WRITE_FORCE_MOMENT_PER_ELEMENTPAIR (default: false) "write visualization output for forces/moments separately for each element pair" // WRITE_UIDS (default: false) "write out the unique ID's for each visualization point,i.e., master and slave beam element global ID (uid_0_beam_1_gid, uid_1_beam_2_gid) and local Gauss point ID (uid_2_gp_id)" .. _SECbinningstrategy: BINNING STRATEGY ---------------- BINNING STRATEGY :: // BINNING STRATEGY: // BIN_SIZE_LOWER_BOUND (default: -1) "Lower bound for bin size. Exact bin size is computed via (Domain edge length)/BIN_SIZE_LOWER_BOUND. This also determines the number of bins in each spatial direction" // BIN_PER_DIR (default: -1 -1 -1) "Number of bins per direction (x, y, z) in particle simulations. Either Define this value or BIN_SIZE_LOWER_BOUND" // PERIODICONOFF (default: 0 0 0) "Turn on/off periodic boundary conditions in each spatial direction" // DOMAINBOUNDINGBOX (default: 1e12 1e12 1e12 1e12 1e12 1e12) "Bounding box for computational domain using binning strategy. Specify diagonal corner points" // WRITEBINS (default: none) "Write none, row or column bins for visualization" .. _SECbiofilmcontrol: BIOFILM CONTROL --------------- BIOFILM CONTROL :: // BIOFILM CONTROL: // BIOFILMGROWTH (default: false) "Scatra algorithm for biofilm growth" // AVGROWTH (default: false) "The calculation of growth parameters is based on averaged values" // FLUXCOEF (default: 0) "Coefficient for growth due to scalar flux" // NORMFORCEPOSCOEF (default: 0) "Coefficient for erosion due to traction normal surface forces" // NORMFORCENEGCOEF (default: 0) "Coefficient for erosion due to compression normal surface forces" // TANGONEFORCECOEF (default: 0) "Coefficient for erosion due to the first tangential surface force" // TANGTWOFORCECOEF (default: 0) "Coefficient for erosion due to the second tangential surface force" // BIOTIMESTEP (default: 0.05) "Time step size for biofilm growth" // BIONUMSTEP (default: 0) "Maximum number of steps for biofilm growth" // OUTPUT_GMSH (default: false) "Do you want to write Gmsh postprocessing files?" .. _SECboundingvolumestrategy: BOUNDINGVOLUME STRATEGY ----------------------- BOUNDINGVOLUME STRATEGY :: // BOUNDINGVOLUME STRATEGY: // BEAM_RADIUS_EXTENSION_FACTOR (default: 2) "Beams radius is multiplied with the factor and then the bounding volume only depending on the beam centerline is extended in all directions (+ and -) by that value." // SPHERE_RADIUS_EXTENSION_FACTOR (default: 2) "Bounding volume of the sphere is the sphere center extended by this factor times the sphere radius in all directions (+ and -)." // WRITE_GEOMETRIC_SEARCH_VISUALIZATION (default: false) "If visualization output for the geometric search should be written" .. _SECbrowniandynamics: BROWNIAN DYNAMICS ----------------- BROWNIAN DYNAMICS :: // BROWNIAN DYNAMICS: // BROWNDYNPROB (default: false) "switch Brownian dynamics on/off" // VISCOSITY (default: 0) "viscosity" // KT (default: 0) "thermal energy" // MAXRANDFORCE (default: -1) "Any random force beyond MAXRANDFORCE*(standard dev.) will be omitted and redrawn. -1.0 means no bounds.'" // TIMESTEP (default: -1) "Within this time interval the random numbers remain constant. -1.0" // BEAMS_DAMPING_COEFF_SPECIFIED_VIA (default: Cylinder_geometry_approx) (choices: Cylinder_geometry_approx|Input_file|cylinder_geometry_approx|input_file) "In which way are damping coefficient values for beams specified?" // BEAMS_DAMPING_COEFF_PER_UNITLENGTH (default: 0.0 0.0 0.0) "values for beam damping coefficients (per unit length and NOT yet multiplied by fluid viscosity): translational perpendicular/parallel to beam axis, rotational around axis" .. _SECcardiacmonodomaincontrol: CARDIAC MONODOMAIN CONTROL -------------------------- CARDIAC MONODOMAIN CONTROL :: // CARDIAC MONODOMAIN CONTROL: // WRITEMAXINTSTATE (default: 0) "number of maximal internal state variables to be postprocessed" // WRITEMAXIONICCURRENTS (default: 0) "number of maximal ionic currents to be postprocessed" // ACTTHRES (default: 1) "threshold for the potential for computing and postprocessing activation time" .. _SECcardiovascular0d-structurecoupling: CARDIOVASCULAR 0D-STRUCTURE COUPLING ------------------------------------ CARDIOVASCULAR 0D-STRUCTURE COUPLING :: // CARDIOVASCULAR 0D-STRUCTURE COUPLING: // TOL_CARDVASC0D_RES (default: 1e-08) "tolerance in the cardiovascular0d error norm for the Newton iteration" // TOL_CARDVASC0D_DOFINCR (default: 1e-08) "tolerance in the cardiovascular0d dof increment error norm for the Newton iteration" // TIMINT_THETA (default: 0.5) "theta for one-step-theta time-integration scheme of Cardiovascular0D" // RESTART_WITH_CARDVASC0D (default: false) "Must be chosen if a non-cardiovascular0d simulation is to be restarted as cardiovascular0d-structural coupled problem." // ENHANCED_OUTPUT (default: false) "Set to yes for enhanced output (like e.g. derivative information)" // LINEAR_COUPLED_SOLVER (default: -1) "number of linear solver used for cardiovascular 0D-structural problems" // SOLALGORITHM (default: direct) (choices: block|direct) // T_PERIOD (default: -1) "periodic time" // EPS_PERIODIC (default: 1e-16) "tolerance for periodic state" // PTC_3D0D (default: false) "Set to yes for doing PTC 2x2 block system." // K_PTC (default: 0) "PTC parameter: 0 means normal Newton, ->infty means steepest desc" .. _SECcardiovascular0d-structurecoupling_respiratoryparameters: CARDIOVASCULAR 0D-STRUCTURE COUPLING/RESPIRATORY PARAMETERS ----------------------------------------------------------- CARDIOVASCULAR 0D-STRUCTURE COUPLING/RESPIRATORY PARAMETERS :: // CARDIOVASCULAR 0D-STRUCTURE COUPLING/RESPIRATORY PARAMETERS: // RESPIRATORY_MODEL (default: None) (choices: None|Standard) // L_alv (default: 0) "alveolar inertance" // R_alv (default: 0) "alveolar resistance" // E_alv (default: 0) "alveolar elastance" // V_lung_tidal (default: 400000) "tidal volume (the total volume of inspired air, in a single breath)" // V_lung_dead (default: 150000) "dead space volume" // V_lung_u (default: 0) "unstressed lung volume (volume of the lung when it is fully collapsed outside the body)" // U_t_curve (default: -1) "time-varying, prescribed pleural pressure curve driven by diaphragm" // U_m (default: 0) "in-breath pressure" // fCO2_ext (default: 0.03) "atmospheric CO2 gas fraction" // fO2_ext (default: 0.21) "atmospheric O2 gas fraction" // kappa_CO2 (default: 0) "diffusion coefficient for CO2 across the hemato-alveolar membrane, in molar value / (time * pressure)" // kappa_O2 (default: 0) "diffusion coefficient for O2 across the hemato-alveolar membrane, in molar value / (time * pressure)" // V_m_gas (default: 22400) "molar volume of an ideal gas" // p_vap_water_37 (default: 6.27949) "vapor pressure of water at 37 degrees celsius" // alpha_CO2 (default: 0) "CO2 solubility constant, in molar value / (volume * pressure)" // alpha_O2 (default: 0) "O2 solubility constant, in molar value / (volume * pressure)" // c_Hb (default: 0) "hemoglobin concentration of the blood, in molar value / volume" // M_CO2_arspl (default: 0) "splanchnic metabolic rate of CO2 production" // M_O2_arspl (default: 0) "splanchnic metabolic rate of O2 consumption" // M_CO2_arespl (default: 0) "extra-splanchnic metabolic rate of CO2 production" // M_O2_arespl (default: 0) "extra-splanchnic metabolic rate of O2 consumption" // M_CO2_armsc (default: 0) "muscular metabolic rate of CO2 production" // M_O2_armsc (default: 0) "muscular metabolic rate of O2 consumption" // M_CO2_arcer (default: 0) "cerebral metabolic rate of CO2 production" // M_O2_arcer (default: 0) "cerebral metabolic rate of O2 consumption" // M_CO2_arcor (default: 0) "coronary metabolic rate of CO2 production" // M_O2_arcor (default: 0) "coronary metabolic rate of O2 consumption" // V_tissspl (default: 1) "splanchnic tissue volume" // V_tissespl (default: 1) "extra-splanchnic tissue volume" // V_tissmsc (default: 1) "muscular tissue volume" // V_tisscer (default: 1) "cerebral tissue volume" // V_tisscor (default: 1) "coronary tissue volume" // V_alv_0 (default: -1) "initial alveolar volume" // q_alv_0 (default: 0) "initial alveolar flux" // p_alv_0 (default: -1) "initial alveolar pressure" // fCO2_alv_0 (default: 0.05263) "initial alveolar CO2 fraction" // fO2_alv_0 (default: 0.1368) "initial alveolar O2 fraction" // q_arspl_sys_in_0 (default: 0) "initial arterial splanchnic in-flux" // q_arsspl_sys_in_0 (default: 0) "initial arterial extra-splanchnic in-flux" // q_armsc_sys_in_0 (default: 0) "initial arterial muscular in-flux" // q_arcer_sys_in_0 (default: 0) "initial arterial cerebral in-flux" // q_arcor_sys_in_0 (default: 0) "initial arterial coronary in-flux" // ppCO2_at_r_0 (default: 1) "initial right atrial CO2 partial pressure" // ppO2_at_r_0 (default: 1) "initial right atrial O2 partial pressure" // ppCO2_v_r_0 (default: 1) "initial right ventricular CO2 partial pressure" // ppO2_v_r_0 (default: 1) "initial right ventricular O2 partial pressure" // ppCO2_ar_pul_0 (default: 1) "initial pulmonary arterial CO2 partial pressure" // ppO2_ar_pul_0 (default: 1) "initial pulmonary arterial O2 partial pressure" // ppCO2_cap_pul_0 (default: 1) "initial pulmonary capillary CO2 partial pressure" // ppO2_cap_pul_0 (default: 1) "initial pulmonary capillary O2 partial pressure" // ppCO2_ven_pul_0 (default: 1) "initial pulmonary venous CO2 partial pressure" // ppO2_ven_pul_0 (default: 1) "initial pulmonary venous O2 partial pressure" // ppCO2_at_l_0 (default: 1) "initial left atrial CO2 partial pressure" // ppO2_at_l_0 (default: 1) "initial left atrial O2 partial pressure" // ppCO2_v_l_0 (default: 1) "initial left ventricular CO2 partial pressure" // ppO2_v_l_0 (default: 1) "initial left ventricular O2 partial pressure" // ppCO2_ar_sys_0 (default: 1) "initial systemic arterial CO2 partial pressure" // ppO2_ar_sys_0 (default: 1) "initial systemic arterial O2 partial pressure" // ppCO2_arspl_sys_0 (default: 1) "initial systemic arterial splanchnic CO2 partial pressure" // ppO2_arspl_sys_0 (default: 1) "initial systemic arterial splanchnic O2 partial pressure" // ppCO2_arespl_sys_0 (default: 1) "initial systemic arterial extra-splanchnic CO2 partial pressure" // ppO2_arespl_sys_0 (default: 1) "initial systemic arterial extra-splanchnic O2 partial pressure" // ppCO2_armsc_sys_0 (default: 1) "initial systemic arterial muscular CO2 partial pressure" // ppO2_armsc_sys_0 (default: 1) "initial systemic arterial muscular O2 partial pressure" // ppCO2_arcer_sys_0 (default: 1) "initial systemic arterial cerebral CO2 partial pressure" // ppO2_arcer_sys_0 (default: 1) "initial systemic arterial cerebral O2 partial pressure" // ppCO2_arcor_sys_0 (default: 1) "initial systemic arterial coronary CO2 partial pressure" // ppO2_arcor_sys_0 (default: 1) "initial systemic arterial coronary O2 partial pressure" // ppCO2_venspl_sys_0 (default: 1) "initial systemic venous splanchnic CO2 partial pressure" // ppO2_venspl_sys_0 (default: 1) "initial systemic venous splanchnic O2 partial pressure" // ppCO2_venespl_sys_0 (default: 1) "initial systemic venous extra-splanchnic CO2 partial pressure" // ppO2_venespl_sys_0 (default: 1) "initial systemic venous extra-splanchnic O2 partial pressure" // ppCO2_venmsc_sys_0 (default: 1) "initial systemic venous muscular CO2 partial pressure" // ppO2_venmsc_sys_0 (default: 1) "initial systemic venous muscular O2 partial pressure" // ppCO2_vencer_sys_0 (default: 1) "initial systemic venous cerebral CO2 partial pressure" // ppO2_vencer_sys_0 (default: 1) "initial systemic venous cerebral O2 partial pressure" // ppCO2_vencor_sys_0 (default: 1) "initial systemic venous coronary CO2 partial pressure" // ppO2_vencor_sys_0 (default: 1) "initial systemic venous coronary O2 partial pressure" // ppCO2_ven_sys_0 (default: 1) "initial systemic venous CO2 partial pressure" // ppO2_ven_sys_0 (default: 1) "initial systemic venous O2 partial pressure" .. _SECcardiovascular0d-structurecoupling_sys-pulcirculationparameters: CARDIOVASCULAR 0D-STRUCTURE COUPLING/SYS-PUL CIRCULATION PARAMETERS ------------------------------------------------------------------- CARDIOVASCULAR 0D-STRUCTURE COUPLING/SYS-PUL CIRCULATION PARAMETERS :: // CARDIOVASCULAR 0D-STRUCTURE COUPLING/SYS-PUL CIRCULATION PARAMETERS: // R_arvalve_max_l (default: 0) "maximal left arterial (semilunar) valve resistance" // R_arvalve_min_l (default: 0) "minimal left arterial (semilunar) valve resistance" // R_atvalve_max_l (default: 0) "maximal left atrial (atrioventricular) valve resistance" // R_atvalve_min_l (default: 0) "minimal left atrial (atrioventricular) valve resistance" // R_arvalve_max_r (default: 0) "maximal right arterial (semilunar) valve resistance" // R_arvalve_min_r (default: 0) "minimal right arterial (semilunar) valve resistance" // R_atvalve_max_r (default: 0) "maximal right atrial (atrioventricular) valve resistance" // R_atvalve_min_r (default: 0) "minimal right atrial (atrioventricular) valve resistance" // ATRIUM_MODEL (default: 0D) (choices: 0D|3D|prescribed) // Atrium_act_curve_l (default: -1) "left atrial activation curve (ONLY for ATRIUM_MODEL '0D'!)" // Atrium_act_curve_r (default: -1) "right atrial activation curve (ONLY for ATRIUM_MODEL '0D'!)" // Atrium_prescr_E_curve_l (default: -1) "left atrial elastance prescription curve (ONLY for ATRIUM_MODEL 'prescribed'!)" // Atrium_prescr_E_curve_r (default: -1) "right atrial elastance prescription curve (ONLY for ATRIUM_MODEL 'prescribed'!)" // E_at_max_l (default: 0) "0D maximum left atrial elastance" // E_at_min_l (default: 0) "0D baseline left atrial elastance" // E_at_max_r (default: 0) "0D maximum right atrial elastance" // E_at_min_r (default: 0) "0D baseline right atrial elastance" // VENTRICLE_MODEL (default: 3D) (choices: 0D|3D|prescribed) // Ventricle_act_curve_l (default: -1) "left ventricular activation curve (ONLY for VENTRICLE_MODEL '0D'!)" // Ventricle_act_curve_r (default: -1) "right ventricular activation curve (ONLY for VENTRICLE_MODEL '0D'!)" // Ventricle_prescr_E_curve_l (default: -1) "left ventricular elastance prescription curve (ONLY for VENTRICLE_MODEL 'prescribed'!)" // Ventricle_prescr_E_curve_r (default: -1) "right ventricular elastance prescription curve (ONLY for VENTRICLE_MODEL 'prescribed'!)" // E_v_max_l (default: 0) "0D maximum left ventricular elastance" // E_v_min_l (default: 0) "0D baseline left ventricular elastance" // E_v_max_r (default: 0) "0D maximum right ventricular elastance" // E_v_min_r (default: 0) "0D baseline right ventricular elastance" // C_ar_sys (default: 0) "systemic arterial compliance" // R_ar_sys (default: 0) "systemic arterial resistance" // L_ar_sys (default: 0) "systemic arterial inertance" // Z_ar_sys (default: 0) "systemic arterial impedance" // C_ar_pul (default: 0) "pulmonary arterial compliance" // R_ar_pul (default: 0) "pulmonary arterial resistance" // L_ar_pul (default: 0) "pulmonary arterial inertance" // Z_ar_pul (default: 0) "pulmonary arterial impedance" // C_ven_sys (default: 0) "systemic venous compliance" // R_ven_sys (default: 0) "systemic venous resistance" // L_ven_sys (default: 0) "systemic venous inertance" // C_ven_pul (default: 0) "pulmonary venous compliance" // R_ven_pul (default: 0) "pulmonary venous resistance" // L_ven_pul (default: 0) "pulmonary venous inertance" // q_vin_l_0 (default: 0) "initial left ventricular in-flux" // p_at_l_0 (default: 0) "initial left atrial pressure" // q_vout_l_0 (default: 0) "initial left ventricular out-flux" // p_v_l_0 (default: 0) "initial left ventricular pressure" // p_ar_sys_0 (default: 0) "initial systemic arterial pressure" // q_ar_sys_0 (default: 0) "initial systemic arterial flux" // p_ven_sys_0 (default: 0) "initial systemic venous pressure" // q_ven_sys_0 (default: 0) "initial systemic venous flux" // q_vin_r_0 (default: 0) "initial right ventricular in-flux" // p_at_r_0 (default: 0) "initial right atrial pressure" // q_vout_r_0 (default: 0) "initial right ventricular out-flux" // p_v_r_0 (default: 0) "initial right ventricular pressure" // p_ar_pul_0 (default: 0) "initial pulmonary arterial pressure" // q_ar_pul_0 (default: 0) "initial pulmonary arterial flux" // p_ven_pul_0 (default: 0) "initial pulmonary venous pressure" // q_ven_pul_0 (default: 0) "initial pulmonary venous flux" // V_at_l_u (default: 0) "unstressed volume of left 0D atrium" // V_v_l_u (default: 0) "unstressed volume of left 0D ventricle" // V_ar_sys_u (default: 0) "unstressed volume of systemic arteries and capillaries" // V_ven_sys_u (default: 100000) "unstressed volume of systemic veins" // V_at_r_u (default: 0) "unstressed volume of right 0D atrium" // V_v_r_u (default: 0) "unstressed volume of right 0D ventricle" // V_ar_pul_u (default: 0) "unstressed volume of pulmonary arteries and capillaries" // V_ven_pul_u (default: 120000) "unstressed volume of pulmonary veins" // C_arspl_sys (default: 0) "systemic arterial splanchnic compliance" // R_arspl_sys (default: 0) "systemic arterial splanchnic resistance" // C_arespl_sys (default: 0) "systemic arterial extra-splanchnic compliance" // R_arespl_sys (default: 0) "systemic arterial extra-splanchnic resistance" // C_armsc_sys (default: 0) "systemic arterial muscular compliance" // R_armsc_sys (default: 0) "systemic arterial muscular resistance" // C_arcer_sys (default: 0) "systemic arterial cerebral compliance" // R_arcer_sys (default: 0) "systemic arterial cerebral resistance" // C_arcor_sys (default: 0) "systemic arterial coronary compliance" // R_arcor_sys (default: 0) "systemic arterial coronary resistance" // C_venspl_sys (default: 0) "systemic venous splanchnic compliance" // R_venspl_sys (default: 0) "systemic venous splanchnic resistance" // C_venespl_sys (default: 0) "systemic venous extra-splanchnic compliance" // R_venespl_sys (default: 0) "systemic venous extra-splanchnic resistance" // C_venmsc_sys (default: 0) "systemic venous muscular compliance" // R_venmsc_sys (default: 0) "systemic venous muscular resistance" // C_vencer_sys (default: 0) "systemic venous cerebral compliance" // R_vencer_sys (default: 0) "systemic venous cerebral resistance" // C_vencor_sys (default: 0) "systemic venous coronary compliance" // R_vencor_sys (default: 0) "systemic venous coronary resistance" // C_cap_pul (default: 0) "pulmonary capillary compliance" // R_cap_pul (default: 0) "pulmonary capillary resistance" // p_arperi_sys_0 (default: 0) "initial systemic peripheral arterial pressure" // q_arspl_sys_0 (default: 0) "initial systemic arterial splanchnic flux" // q_arespl_sys_0 (default: 0) "initial systemic arterial extra-splanchnic flux" // q_armsc_sys_0 (default: 0) "initial systemic arterial muscular flux" // q_arcer_sys_0 (default: 0) "initial systemic arterial cerebral flux" // q_arcor_sys_0 (default: 0) "initial systemic arterial coronary flux" // p_venspl_sys_0 (default: 0) "initial systemic venous splanchnic pressure" // q_venspl_sys_0 (default: 0) "initial systemic venous splanchnic flux" // p_venespl_sys_0 (default: 0) "initial systemic venous extra-splanchnic pressure" // q_venespl_sys_0 (default: 0) "initial systemic venous extra-splanchnic flux" // p_venmsc_sys_0 (default: 0) "initial systemic venous muscular pressure" // q_venmsc_sys_0 (default: 0) "initial systemic venous muscular flux" // p_vencer_sys_0 (default: 0) "initial systemic venous cerebral pressure" // q_vencer_sys_0 (default: 0) "initial systemic venous cerebral flux" // p_vencor_sys_0 (default: 0) "initial systemic venous coronary pressure" // q_vencor_sys_0 (default: 0) "initial systemic venous coronary flux" // p_cap_pul_0 (default: 0) "initial pulmonary capillary pressure" // q_cap_pul_0 (default: 0) "initial pulmonary capillary flux" // V_arspl_sys_u (default: 274400) "unstressed volume of systemic splanchnic arteries" // V_arespl_sys_u (default: 134640) "unstressed volume of systemic extra-splanchnic arteries" // V_armsc_sys_u (default: 105800) "unstressed volume of systemic muscular arteries" // V_arcer_sys_u (default: 72130) "unstressed volume of systemic cerebral arteries" // V_arcor_sys_u (default: 24000) "unstressed volume of systemic coronary arteries" // V_venspl_sys_u (default: 1.121e+06) "unstressed volume of systemic splanchnic veins" // V_venespl_sys_u (default: 550000) "unstressed volume of systemic extra-splanchnic veins" // V_venmsc_sys_u (default: 432140) "unstressed volume of systemic muscular veins" // V_vencer_sys_u (default: 294640) "unstressed volume of systemic cerebral veins" // V_vencor_sys_u (default: 98210) "unstressed volume of systemic coronary veins" // V_cap_pul_u (default: 123000) "unstressed volume of pulmonary capillaries" .. _SECcontactdynamic: CONTACT DYNAMIC --------------- CONTACT DYNAMIC :: // CONTACT DYNAMIC: // LINEAR_SOLVER (default: -1) "number of linear solver used for meshtying and contact" // RESTART_WITH_CONTACT (default: false) "Must be chosen if a non-contact simulation is to be restarted with contact" // ADHESION (default: none) "Type of adhesion law" // FRICTION (default: None) (choices: Coulomb|None|Stick|Tresca) "Type of friction law" // FRLESS_FIRST (default: false) "If chosen the first time step of a newly in contact slave node is regarded as frictionless" // GP_SLIP_INCR (default: false) "If chosen the slip increment is computed gp-wise which results to a non-objective quantity, but this would be consistent to wear and tsi calculations." // STRATEGY (default: Lagrange) (choices: Ehl|Lagrange|LagrangianMultipliers|MultiScale|Nitsche|Penalty|Uzawa|lagrange|penalty) "Type of employed solving strategy" // SYSTEM (default: Condensed) (choices: Condensed|Condensedlagmult|SaddlePoint|Saddlepoint|cond|condensed|condensedlagmult|condlm|none|saddlepoint|sp) "Type of linear system setup / solution" // PENALTYPARAM (default: 0) "Penalty parameter for penalty / Uzawa augmented solution strategy" // PENALTYPARAMTAN (default: 0) "Tangential penalty parameter for penalty / Uzawa augmented solution strategy" // UZAWAMAXSTEPS (default: 10) "Maximum no. of Uzawa steps for Uzawa solution strategy" // UZAWACONSTRTOL (default: 1e-08) "Tolerance of constraint norm for Uzawa solution strategy" // SEMI_SMOOTH_NEWTON (default: true) "If chosen semi-smooth Newton concept is applied" // SEMI_SMOOTH_CN (default: 1) "Weighting factor cn for semi-smooth PDASS" // SEMI_SMOOTH_CT (default: 1) "Weighting factor ct for semi-smooth PDASS" // CONTACTFORCE_ENDTIME (default: false) "If chosen, the contact force is not evaluated at the generalized midpoint, but at the end of the time step" // VELOCITY_UPDATE (default: false) "If chosen, velocity update method is applied" // INITCONTACTBYGAP (default: false) "Initialize init contact by weighted gap vector" // INITCONTACTGAPVALUE (default: 0) "Value for initialization of init contact set with gap vector" // NORMCOMBI_RESFCONTCONSTR (default: And) (choices: And|Or) "binary operator to combine contact constraints and residual force values" // NORMCOMBI_DISPLAGR (default: And) (choices: And|Or) "binary operator to combine displacement increments and Lagrange multiplier increment values" // TOLCONTCONSTR (default: 1e-06) "tolerance in the contact constraint norm for the newton iteration (saddlepoint formulation only)" // TOLLAGR (default: 1e-06) "tolerance in the LM norm for the newton iteration (saddlepoint formulation only)" // CONSTRAINT_DIRECTIONS (default: ntt) "formulation of constraints in normal/tangential or xyz-direction" // NONSMOOTH_GEOMETRIES (default: false) "If chosen the contact algorithm combines mortar and nts formulations. This is needed if contact between entities of different geometric dimension (such as contact between surfaces and lines, or lines and nodes) can occur" // NONSMOOTH_CONTACT_SURFACE (default: false) "This flag is used to alter the criterion for the evaluation of the so-called qualified vectors in the case of a self contact scenario. This is needed as the standard criterion is only valid for smooth surfaces and thus has to be altered, if the surface that is defined to be a self contact surface is non-smooth!" // HYBRID_ANGLE_MIN (default: -1) "Non-smooth contact: angle between cpp normal and element normal: begin transition (Mortar)" // HYBRID_ANGLE_MAX (default: -1) "Non-smooth contact: angle between cpp normal and element normal: end transition (NTS)" // CPP_NORMALS (default: false) "If chosen the nodal normal field is created as averaged CPP normal field." // TIMING_DETAILS (default: false) "Enable and print detailed contact timings to screen." // NITSCHE_THETA (default: 0) "+1: symmetric, 0: non-symmetric, -1: skew-symmetric" // NITSCHE_THETA_2 (default: 1) "+1: Chouly-type, 0: Burman penalty-free (only with theta=-1)" // NITSCHE_WEIGHTING (default: harmonic) "how to weight consistency terms in Nitsche contact formulation" // NITSCHE_PENALTY_ADAPTIVE (default: true) "adapt penalty parameter after each converged time step" // REGULARIZED_NORMAL_CONTACT (default: false) "add a regularized normal contact formulation" // REGULARIZATION_THICKNESS (default: -1) "maximum contact penetration" // REGULARIZATION_STIFFNESS (default: -1) "initial contact stiffness (i.e. initial \"penalty parameter\")" .. _SECcoupledreduced-dairwaysandtissuedynamic: COUPLED REDUCED-D AIRWAYS AND TISSUE DYNAMIC -------------------------------------------- COUPLED REDUCED-D AIRWAYS AND TISSUE DYNAMIC :: // COUPLED REDUCED-D AIRWAYS AND TISSUE DYNAMIC: // CONVTOL_P (default: 1e-06) "Coupled red_airway and tissue iteration convergence for pressure" // CONVTOL_Q (default: 1e-06) "Coupled red_airway and tissue iteration convergence for flux" // MAXITER (default: 5) "Maximum coupling iterations" // RELAXTYPE (default: norelaxation) "Dynamic Relaxation Type" // TIMESTEP (default: 0.01) "Time increment dt" // NUMSTEP (default: 1) "Number of Time Steps" // MAXTIME (default: 4) // NORMAL (default: 1) .. _SECcutgeneral: CUT GENERAL ----------- CUT GENERAL :: // CUT GENERAL: // KERNEL_INTERSECTION_FLOATTYPE (default: double) (choices: cln|double) "The floattype of the cut surface-edge intersection" // KERNEL_DISTANCE_FLOATTYPE (default: double) (choices: cln|double) "The floattype of the cut distance computation" // GENERAL_POSITION_DISTANCE_FLOATTYPE (default: none) (choices: cln|double|none) "A general floattype for Cut::Position for Embedded Elements (compute_distance)" // GENERAL_POSITION_POSITION_FLOATTYPE (default: none) (choices: cln|double|none) "A general floattype for Cut::Position Elements (ComputePosition)" // DIRECT_DIVERGENCE_REFPLANE (default: all) (choices: all|diagonal|diagonal_side|facet|none|side) "Specify which Referenceplanes are used in DirectDivergence" // SPLIT_CUTSIDES (default: true) "Split Quad4 CutSides into Tri3-Subtriangles?" // DO_SELFCUT (default: true) "Do the SelfCut?" // SELFCUT_DO_MESHCORRECTION (default: true) "Do meshcorrection in the SelfCut?" // SELFCUT_MESHCORRECTION_MULTIPLICATOR (default: 30) "ISLANDS with maximal size of the bounding box of h*multiplacator will be removed in the meshcorrection" // BOUNDARYCELL_CUBATURDEGREE (default: 20) "Cubaturedegree utilized for the numerical integration on the CUT BoundaryCells." // INTEGRATE_INSIDE_CELLS (default: true) "Should the integration be done on inside cells" .. _SECdiscretisation: DISCRETISATION -------------- DISCRETISATION :: // DISCRETISATION: // NUMFLUIDDIS (default: 1) "Number of meshes in fluid field" // NUMSTRUCDIS (default: 1) "Number of meshes in structural field" // NUMALEDIS (default: 1) "Number of meshes in ale field" // NUMARTNETDIS (default: 1) "Number of meshes in arterial network field" // NUMTHERMDIS (default: 1) "Number of meshes in thermal field" // NUMAIRWAYSDIS (default: 1) "Number of meshes in reduced dimensional airways network field" .. _SECelastohydrodynamic: ELASTO HYDRO DYNAMIC -------------------- ELASTO HYDRO DYNAMIC :: // ELASTO HYDRO DYNAMIC: // RESTARTEVERYTIME (default: 0) "write restart possibility every RESTARTEVERY steps" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: -1) "time step size dt" // DIFFTIMESTEPSIZE (default: false) "use different step size for lubrication and solid" // RESULTSEVERYTIME (default: 0) "increment for writing solution" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // ITEMIN (default: 1) "minimal number of iterations over fields" // FIELDCOUPLING (default: none) (choices: matching|none) "Type of coupling strategy between fields" // COUPALGO (default: ehl_Monolithic) "Coupling strategies for EHL solvers" // UNPROJ_ZERO_DBC (default: false) "set unprojectable nodes to zero pressure via Dirichlet condition" // DRY_CONTACT_MODEL (default: false) "set unprojectable nodes to zero pressure via Dirichlet condition" .. _SECelastohydrodynamic_monolithic: ELASTO HYDRO DYNAMIC/MONOLITHIC ------------------------------- ELASTO HYDRO DYNAMIC/MONOLITHIC :: // ELASTO HYDRO DYNAMIC/MONOLITHIC: // CONVTOL (default: 1e-06) "tolerance for convergence check of EHL" // TOLINC (default: 1e-06) "tolerance for convergence check of EHL-increment in monolithic EHL" // NORM_RESF (default: Abs) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // NORM_INC (default: Abs) (choices: Abs|Mix|Rel) "type of norm for convergence check of primary variables in EHL" // NORMCOMBI_RESFINC (default: Coupl_And_Single) (choices: And|And_Single|Coupl_And_Single|Coupl_Or_Single|Or|Or_Single) "binary operator to combine primary variables and residual force values" // ITERNORM (default: Rms) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // PTCDT (default: 0.1) "pseudo time step for pseudo-transient continuation (PTC) stabilised Newton procedure" // LINEAR_SOLVER (default: -1) "number of linear solver used for monolithic EHL problems" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // INFNORMSCALING (default: true) "Scale blocks of matrix with row infnorm?" .. _SECelastohydrodynamic_partitioned: ELASTO HYDRO DYNAMIC/PARTITIONED -------------------------------- ELASTO HYDRO DYNAMIC/PARTITIONED :: // ELASTO HYDRO DYNAMIC/PARTITIONED: // MAXOMEGA (default: 10) "largest omega allowed for Aitken relaxation" // MINOMEGA (default: 0.1) "smallest omega allowed for Aitken relaxation" // STARTOMEGA (default: 1) "fixed relaxation parameter" // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteration within partitioned EHL" .. _SECelchcontrol: ELCH CONTROL ------------ ELCH CONTROL :: // ELCH CONTROL: // MOVBOUNDARYITEMAX (default: 10) "Maximum number of outer iterations in electrode shape change computations" // MOVBOUNDARYCONVTOL (default: 1e-06) "Convergence check tolerance for outer loop in electrode shape change computations" // TEMPERATURE (default: 298) "Constant temperature (Kelvin)" // TEMPERATURE_FROM_FUNCT (default: -1) "Homogeneous temperature within electrochemistry field that can be time dependent according to function definition" // FARADAY_CONSTANT (default: 96485.3) "Faraday constant (in unit system as chosen in input file)" // GAS_CONSTANT (default: 8.31447) "(universal) gas constant (in unit system as chosen in input file)" // MOVINGBOUNDARY (default: No) (choices: No|fully-transient|pseudo-transient) "ELCH algorithm for deforming meshes" // MOLARVOLUME (default: 0) "Molar volume for electrode shape change computations" // MOVBOUNDARYTHETA (default: 0) "One-step-theta factor in electrode shape change computations" // GALVANOSTATIC (default: false) "flag for galvanostatic mode" // GSTAT_APPROX_ELECT_RESIST (default: relation_pot_cur) (choices: effective_length_with_initial_cond|effective_length_with_integrated_cond|relation_pot_cur) "relation of potential and current flow" // GSTATCONDID_CATHODE (default: 0) "condition id of electrode kinetics for cathode" // GSTATCONDID_ANODE (default: 1) "condition id of electrode kinetics for anode" // GSTATCONVTOL (default: 1e-05) "Convergence check tolerance for galvanostatic mode" // GSTATCURTOL (default: 1e-15) "Current Tolerance" // GSTATFUNCTNO (default: -1) "function number defining the imposed current curve" // GSTATITEMAX (default: 10) "maximum number of iterations for galvanostatic mode" // GSTAT_LENGTH_CURRENTPATH (default: 0) "average length of the current path" // EQUPOT (default: Undefined) (choices: ENC|ENC_PDE|ENC_PDE_ELIM|Laplace|Poisson|Undefined|divi) "type of closing equation for electric potential" // DIFFCOND_FORMULATION (default: false) "Activation of diffusion-conduction formulation" // INITPOTCALC (default: false) "Automatically calculate initial field for electric potential" // ONLYPOTENTIAL (default: false) "Coupling of general ion transport equation with Laplace equation" // COUPLE_BOUNDARY_FLUXES (default: true) "Coupling of lithium-ion flux density and electric current density at Dirichlet and Neumann boundaries" // CYCLING_TIMESTEP (default: -1) "modified time step size for CCCV cell cycling" // ELECTRODE_INFO_EVERY_STEP (default: false) "the cell voltage, SOC, and C-Rate will be written to the csv file every step, even if RESULTSEVERY is not 1" .. _SECelchcontrol_diffcond: ELCH CONTROL/DIFFCOND --------------------- ELCH CONTROL/DIFFCOND :: // ELCH CONTROL/DIFFCOND: // CURRENT_SOLUTION_VAR (default: false) "Current as a solution variable" // MAT_DIFFCOND_DIFFBASED (default: true) "Coupling terms of chemical diffusion for current equation are based on t and kappa" // MAT_NEWMAN_CONST_A (default: 2) "Constant A for the Newman model(term for the concentration overpotential)" // MAT_NEWMAN_CONST_B (default: -2) "Constant B for the Newman model(term for the concentration overpotential)" // MAT_NEWMAN_CONST_C (default: -1) "Constant C for the Newman model(term for the concentration overpotential)" // PERMITTIVITY_VACUUM (default: 8.85419e-12) "Vacuum permittivity" .. _SECelchcontrol_scl: ELCH CONTROL/SCL ---------------- ELCH CONTROL/SCL :: // ELCH CONTROL/SCL: // ADD_MICRO_MACRO_COUPLING (default: false) "flag for micro macro coupling with scls" // COUPLING_OUTPUT (default: false) "write coupled node gids and node coordinates to csv file" // INITPOTCALC (default: false) "calculate initial potential field?" // SOLVER (default: -1) "solver for coupled SCL problem" // MATRIXTYPE (default: undefined) (choices: block|sparse|undefined) "type of global system matrix in global system of equations" // ADAPT_TIME_STEP (default: -1) "time step when time step size should be updated to 'ADAPTED_TIME_STEP_SIZE'." // ADAPTED_TIME_STEP_SIZE (default: -1) "new time step size." // INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "Initial Field for scalar transport problem" // INITFUNCNO (default: -1) "function number for scalar transport initial field" .. _SECembeddedmeshcoupling: EMBEDDED MESH COUPLING ---------------------- EMBEDDED MESH COUPLING :: // EMBEDDED MESH COUPLING: // COUPLING_STRATEGY (default: none) "Strategy to couple background and overlapping mesh" // MORTAR_SHAPE_FUNCTION (default: none) "Shape functions that should be use in case of coupling using the Mortar/Lagrange Multiplier method" // CONSTRAINT_ENFORCEMENT (default: none) "Apply a constraint enforcement in the embedded mesh coupling strategy" // CONSTRAINT_ENFORCEMENT_PENALTYPARAM (default: 0) "Penalty parameter for the constraint enforcement in embedded mesh coupling" .. _SECfluidbeaminteraction: FLUID BEAM INTERACTION ---------------------- FLUID BEAM INTERACTION :: // FLUID BEAM INTERACTION: // COUPLING (default: two-way) (choices: fluid|solid|two-way) "Type of FBI coupling" // STARTSTEP (default: 0) "Time Step at which to begin the fluid beam coupling. Usually this will be the first step." // PRESORT_STRATEGY (default: bruteforce) "Presort strategy for the beam elements" .. _SECfluidbeaminteraction_beamtofluidmeshtying: FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING ---------------------------------------------- FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING :: // FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING: // MESHTYING_DISCRETIZATION (default: none) "Type of employed meshtying discretization" // CONSTRAINT_STRATEGY (default: none) "Type of employed constraint enforcement strategy" // PENALTY_PARAMETER (default: 0) "Penalty parameter for beam-to-Fluid volume meshtying" // SEARCH_RADIUS (default: 1000) "Absolute Search radius for beam-to-fluid volume meshtying. Choose carefully to not blow up memory demand but to still find all interaction pairs!" // MORTAR_SHAPE_FUNCTION (default: none) "Shape function for the mortar Lagrange-multipliers" // GEOMETRY_PAIR_STRATEGY (default: segmentation) "Type of employed segmentation strategy" // GEOMETRY_PAIR_SEGMENTATION_SEARCH_POINTS (default: 6) "Number of search points for segmentation" // GEOMETRY_PAIR_SEGMENTATION_NOT_ALL_GAUSS_POINTS_PROJECT_VALID_ACTION (default: fail) "What to do if not all Gauss points of a segment project valid" // GAUSS_POINTS (default: 6) "Number of Gauss Points for the integral evaluations" // INTEGRATION_POINTS_CIRCUMFERENCE (default: 6) "Number of Integration points along the circumferential direction of the beam. This is parameter is only used in beam to cylinder meshtying. No gauss integration is used along the circumferential direction, equally spaced integration points are used." .. _SECfluidbeaminteraction_beamtofluidmeshtying_runtimevtkoutput: FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING/RUNTIME VTK OUTPUT ----------------------------------------------------------------- FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING/RUNTIME VTK OUTPUT :: // FLUID BEAM INTERACTION/BEAM TO FLUID MESHTYING/RUNTIME VTK OUTPUT: // WRITE_OUTPUT (default: false) "Enable / disable beam-to-fluid mesh tying output." // NODAL_FORCES (default: false) "Enable / disable output of the resulting nodal forces due to beam to Fluid interaction." // SEGMENTATION (default: false) "Enable / disable output of segmentation points." // INTEGRATION_POINTS (default: false) "Enable / disable output of used integration points. If the meshtying method has 'forces' at the integration point, they will also be output." // CONSTRAINT_VIOLATION (default: false) "Enable / disable output of the constraint violation into a output_name.penalty csv file." // MORTAR_LAMBDA_DISCRET (default: false) "Enable / disable output of the discrete Lagrange multipliers at the node of the Lagrange multiplier shape functions." // MORTAR_LAMBDA_CONTINUOUS (default: false) "Enable / disable output of the continuous Lagrange multipliers function along the beam." // MORTAR_LAMBDA_CONTINUOUS_SEGMENTS (default: 5) "Number of segments for continuous mortar output" .. _SECfluiddynamic: FLUID DYNAMIC ------------- FLUID DYNAMIC :: // FLUID DYNAMIC: // PHYSICAL_TYPE (default: Incompressible) (choices: Artificial_compressibility|Boussinesq|Incompressible|Loma|Oseen|Stokes|Temp_dep_water|Varying_density|Weakly_compressible|Weakly_compressible_dens_mom|Weakly_compressible_stokes|Weakly_compressible_stokes_dens_mom) "Physical Type" // LINEAR_SOLVER (default: -1) "number of linear solver used for fluid dynamics" // SIMPLER_SOLVER (default: -1) "number of linear solver used for fluid dynamic (ONLY NECESSARY FOR BlockGaussSeidel solver block within fluid mehstying case any more!!!!)" // WSS_TYPE (default: Standard) (choices: Aggregation|Mean|Standard) "which type of stresses and wall shear stress" // WSS_ML_AGR_SOLVER (default: -1) "Set ML-solver number for smoothing of residual-based calculated wallshearstress via plain aggregation." // TIMEINTEGR (default: One_Step_Theta) (choices: Af_Gen_Alpha|BDF2|Np_Gen_Alpha|One_Step_Theta|Stationary) "Time Integration Scheme" // OST_CONT_PRESS (default: Cont_normal_Press_normal) "One step theta option for time discretization of continuity eq. and pressure" // NONLINITER (default: fixed_point_like) "Nonlinear iteration scheme" // PREDICTOR (default: steady_state) (choices: TangVel|constant_acceleration|constant_increment|explicit_second_order_midpoint|steady_state|zero_acceleration) "Predictor for first guess in nonlinear iteration" // CONVCHECK (default: L_2_norm) (choices: L_2_norm) "Norm for convergence check (relative increments and absolute residuals)" // INCONSISTENT_RESIDUAL (default: false) "do not evaluate residual after solution has converged (->faster)" // INITIALFIELD (default: zero_field) (choices: BELTRAMI-FLOW|FLAME_VORTEX_INTERACTION|KIM-MOIN-FLOW|channel_weakly_compressible|disturbed_field_from_function|field_by_function|forced_hit_numeric_spectrum|forced_hit_passive|forced_hit_simple_algebraic_spectrum|hit_comte_bellot_corrsin_initial_field|zero_field) "Initial Field for transport problem" // OSEENFIELDFUNCNO (default: -1) "function number of Oseen advective field" // LIFTDRAG (default: false) "Calculate lift and drag forces along specified boundary" // CONVFORM (default: convective) (choices: conservative|convective) "form of convective term" // NONLINEARBC (default: false) "Flag to activate check for potential nonlinear boundary conditions" // MESHTYING (default: no) (choices: Condensed_Bmat|Condensed_Bmat_merged|Condensed_Smat|no) "Flag to (de)activate mesh tying algorithm" // GRIDVEL (default: BE) "scheme for determination of gridvelocity from displacements" // ALLDOFCOUPLED (default: true) "all dof (incl. pressure) are coupled" // CALCERROR (default: no) "Flag to (de)activate error calculations" // CALCERRORFUNCNO (default: -1) "Function for Error Calculation" // CORRTERMFUNCNO (default: -1) "Function for calculation of the correction term for the weakly compressible problem" // BODYFORCEFUNCNO (default: -1) "Function for calculation of the body force for the weakly compressible problem" // STAB_DEN_REF (default: 0) "Reference stabilization parameter for the density for the HDG weakly compressible formulation" // STAB_MOM_REF (default: 0) "Reference stabilization parameter for the momentum for the HDG weakly compressible formulation" // VARVISCFUNCNO (default: -1) "Function for calculation of a variable viscosity for the weakly compressible problem" // PRESSAVGBC (default: false) "Flag to (de)activate imposition of boundary condition for the considered element average pressure" // REFMACH (default: 1) "Reference Mach number" // BLOCKMATRIX (default: false) "Indicates if system matrix should be assembled into a sparse block matrix type." // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // INFNORMSCALING (default: false) "Scale blocks of matrix with row infnorm?" // GMSH_OUTPUT (default: false) "write output to gmsh files" // COMPUTE_DIVU (default: false) "Compute divergence of velocity field at the element center" // COMPUTE_EKIN (default: false) "Compute kinetic energy at the end of each time step and write it to file." // NEW_OST (default: false) "Solve the Navier-Stokes equation with the new One Step Theta algorithm" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 20) "Increment for writing restart" // NUMSTEP (default: 1) "Total number of Timesteps" // STEADYSTEP (default: -1) "steady state check every step" // NUMSTASTEPS (default: 0) "Number of Steps for Starting Scheme" // STARTFUNCNO (default: -1) "Function for Initial Starting Field" // ITEMAX (default: 10) "max. number of nonlin. iterations" // INITSTATITEMAX (default: 5) "max number of nonlinear iterations for initial stationary solution" // TIMESTEP (default: 0.01) "Time increment dt" // MAXTIME (default: 1000) "Total simulation time" // ALPHA_M (default: 1) "Time integration factor" // ALPHA_F (default: 1) "Time integration factor" // GAMMA (default: 1) "Time integration factor" // THETA (default: 0.66) "Time integration factor" // START_THETA (default: 1) "Time integration factor for starting scheme" // STRONG_REDD_3D_COUPLING_TYPE (default: false) "Flag to (de)activate potential Strong 3D redD coupling" // VELGRAD_PROJ_SOLVER (default: -1) "Number of linear solver used for L2 projection" // VELGRAD_PROJ_METHOD (default: none) (choices: L2_projection|none|superconvergent_patch_recovery) "Flag to (de)activate gradient reconstruction." // OFF_PROC_ASSEMBLY (default: false) "Do not evaluate ghosted elements but communicate them --> faster if element call is expensive" .. _SECfluiddynamic_edge-basedstabilization: FLUID DYNAMIC/EDGE-BASED STABILIZATION -------------------------------------- FLUID DYNAMIC/EDGE-BASED STABILIZATION :: // FLUID DYNAMIC/EDGE-BASED STABILIZATION: // EOS_PRES (default: none) (choices: none|std_eos|xfem_gp) "Flag to (de)activate pressure edge-based stabilization. Options: do not use pressure edge-based stabilization, or, use pressure edge-based stabilization as standard edge-based stabilization on the entire domain, or, use pressure edge-based stabilization as xfem ghost-penalty stabilization just around cut elements" // EOS_CONV_STREAM (default: none) (choices: none|std_eos|xfem_gp) "Flag to (de)activate convective streamline edge-based stabilization. Options: do not use convective streamline edge-based stabilization, or, use convective streamline edge-based stabilization as standard edge-based stabilization on the entire domain, or, use convective streamline edge-based stabilization as xfem ghost-penalty stabilization just around cut elements" // EOS_CONV_CROSS (default: none) (choices: none|std_eos|xfem_gp) "Flag to (de)activate convective crosswind edge-based stabilization. Options:do not use convective crosswind edge-based stabilization, or, use convective crosswind edge-based stabilization as standard edge-based stabilization on the entire domain, or, use convective crosswind edge-based stabilization as xfem ghost-penalty stabilization just around cut elements" // EOS_DIV (default: none) (choices: div_jump_std_eos|div_jump_xfem_gp|none|vel_jump_std_eos|vel_jump_xfem_gp) "Flag to (de)activate divergence edge-based stabilization. Options: do not use divergence edge-based stabilization, or, divergence edge-based stabilization based on velocity jump on the entire domain, or, divergence edge-based stabilization based on divergence jump just around cut elements, or, divergence edge-based stabilization based on velocity jump on the entire domain, or, divergence edge-based stabilization based on divergence jump just around cut elements" // PRES_KRYLOV_2Dz (default: false) "residual based without second derivatives (i.e. only consistent for tau->0, but faster)" // EOS_DEFINITION_TAU (default: Burman_Hansbo_DAngelo_Zunino) (choices: Braack_Burman_John_Lube|Braack_Burman_John_Lube_wo_divjump|Burman|Burman_Fernandez|Burman_Fernandez_Hansbo|Burman_Fernandez_Hansbo_wo_dt|Burman_Hansbo_DAngelo_Zunino|Burman_Hansbo_DAngelo_Zunino_wo_dt|Franca_Barrenechea_Valentin_Wall|Schott_Massing_Burman_DAngelo_Zunino|Schott_Massing_Burman_DAngelo_Zunino_wo_dt|Taylor_Hughes_Zarins_Whiting_Jansen_Codina_scaling|tau_not_defined) "Definition of stabilization parameter for edge-based stabilization" // EOS_H_DEFINITION (default: EOS_he_max_diameter_to_opp_surf) "Definition of element length for edge-based stabilization. Options:take the maximal (nsd-1)D diameter of faces that connect the internal face to its opposite faces, or, take the maximal 1D distance along 1D edge to opposite surface for both parent elements, or, take the maximal (nsd-1)D face diameter of all faces for both parent elements, or, maximal nD diameter of the neighboring elements, or, maximal (n-1)D diameter of the internal face/edge, or, take the maximal volume equivalent diameter of adjacent elements" .. _SECfluiddynamic_multifractalsubgridscales: FLUID DYNAMIC/MULTIFRACTAL SUBGRID SCALES ----------------------------------------- FLUID DYNAMIC/MULTIFRACTAL SUBGRID SCALES :: // FLUID DYNAMIC/MULTIFRACTAL SUBGRID SCALES: // CSGS (default: 0) "Modelparameter of multifractal subgrid-scales." // SCALE_SEPARATION (default: no_scale_sep) (choices: algebraic_multigrid_operator|box_filter|no_scale_sep) "Specify the filter type for scale separation in LES. no_scale_sep: no scale separation. box_filter: classical box filter. algebraic_multigrid_operator: scale separation by algebraic multigrid operator." // ML_SOLVER (default: -1) "Set solver number for scale separation via level set transfer operators from plain aggregation." // CALC_N (default: false) "Flag to (de)activate calculation of N from the Reynolds number." // N (default: 1) "Set grid to viscous scale ratio." // REF_LENGTH (default: cube_edge) (choices: cube_edge|gradient_based|metric_tensor|sphere_diameter|streamlength) "Specify the reference length for Re-dependent N. cube_edge: edge length of volume equivalent cube. sphere_diameter: diameter of volume equivalent sphere. streamlength: streamlength taken from stabilization. gradient_based: gradient based length taken from stabilization. metric_tensor: metric tensor taken from stabilization." // REF_VELOCITY (default: strainrate) (choices: fine_scale|resolved|strainrate) "Specify the reference velocity for Re-dependent N. strainrate: norm of strain rate. resolved: resolved velocity. fine_scale: fine-scale velocity." // C_NU (default: 1) "Proportionality constant between Re and ratio viscous scale to element length." // NEAR_WALL_LIMIT (default: false) "Flag to (de)activate near-wall limit." // EVALUATION_B (default: element_center) (choices: element_center|integration_point) "Location where B is evaluated element_center: evaluate B at element center. integration_point: evaluate B at integration point." // BETA (default: 0) "Cross- and Reynolds-stress terms only on right-hand-side." // CONVFORM (default: convective) (choices: conservative|convective) "form of convective term convective: Use the convective form. conservative: Use the conservative form." // CSGS_PHI (default: 0) "Modelparameter of multifractal subgrid-scales for scalar transport." // ADAPT_CSGS_PHI (default: false) "Flag to (de)activate adaption of CsgsD to CsgsB." // NEAR_WALL_LIMIT_CSGS_PHI (default: false) "Flag to (de)activate near-wall limit for scalar field." // CONSISTENT_FLUID_RESIDUAL (default: false) "Flag to (de)activate the consistency term for residual-based stabilization." // C_DIFF (default: 1) "Proportionality constant between Re*Pr and ratio dissipative scale to element length. Usually equal cnu." // SET_FINE_SCALE_VEL (default: false) "Flag to set fine-scale velocity for parallel nightly tests." // LOMA_CONTI (default: false) "Flag to (de)activate cross- and Reynolds-stress terms in loma continuity equation." .. _SECfluiddynamic_nonlinearsolvertolerances: FLUID DYNAMIC/NONLINEAR SOLVER TOLERANCES ----------------------------------------- FLUID DYNAMIC/NONLINEAR SOLVER TOLERANCES :: // FLUID DYNAMIC/NONLINEAR SOLVER TOLERANCES: // TOL_VEL_RES (default: 1e-06) "Tolerance for convergence check of velocity residual" // TOL_VEL_INC (default: 1e-06) "Tolerance for convergence check of velocity increment" // TOL_PRES_RES (default: 1e-06) "Tolerance for convergence check of pressure residual" // TOL_PRES_INC (default: 1e-06) "Tolerance for convergence check of pressure increment" .. _SECfluiddynamic_porous-flowstabilization: FLUID DYNAMIC/POROUS-FLOW STABILIZATION --------------------------------------- FLUID DYNAMIC/POROUS-FLOW STABILIZATION :: // FLUID DYNAMIC/POROUS-FLOW STABILIZATION: // STAB_BIOT (default: false) "Flag to (de)activate BIOT stabilization." // STAB_BIOT_SCALING (default: 1) "Scaling factor for stabilization parameter for biot stabilization of porous flow." // STABTYPE (default: residual_based) (choices: edge_based|no_stabilization|residual_based) "Apply (un)stabilized fluid formulation. No stabilization is only possible for inf-sup stable elements. Use a residual-based stabilization or, more generally, a stabilization based on the concept of the residual-based variational multiscale method... Expecting additional inputUse an edge-based stabilization, especially for XFEM" // INCONSISTENT (default: false) "residual based without second derivatives (i.e. only consistent for tau->0, but faster)" // Reconstruct_Sec_Der (default: false) "residual computed with a reconstruction of the second derivatives via projection or superconvergent patch recovery" // TDS (default: quasistatic) (choices: quasistatic|time_dependent) "Flag to allow time dependency of subscales for residual-based stabilization." // TRANSIENT (default: no_transient) (choices: no_transient|transient_complete|yes_transient) "Specify how to treat the transient term." // PSPG (default: true) "Flag to (de)activate PSPG stabilization." // SUPG (default: true) "Flag to (de)activate SUPG stabilization." // GRAD_DIV (default: true) "Flag to (de)activate grad-div term." // VSTAB (default: no_vstab) (choices: no_vstab|vstab_gls|vstab_gls_rhs|vstab_usfem|vstab_usfem_rhs) "Flag to (de)activate viscous term in residual-based stabilization." // RSTAB (default: no_rstab) (choices: no_rstab|rstab_gls|rstab_usfem) "Flag to (de)activate reactive term in residual-based stabilization." // CROSS-STRESS (default: no_cross) (choices: cross_rhs|no_cross|yes_cross) "Flag to (de)activate cross-stress term -> residual-based VMM." // REYNOLDS-STRESS (default: no_reynolds) (choices: no_reynolds|reynolds_rhs|yes_reynolds) "Flag to (de)activate Reynolds-stress term -> residual-based VMM." // DEFINITION_TAU (default: Franca_Barrenechea_Valentin_Frey_Wall) (choices: Codina|Codina_wo_dt|Franca_Barrenechea_Valentin_Frey_Wall|Franca_Barrenechea_Valentin_Frey_Wall_wo_dt|Franca_Madureira_Valentin_Badia_Codina|Franca_Madureira_Valentin_Badia_Codina_wo_dt|Shakib_Hughes_Codina|Shakib_Hughes_Codina_wo_dt|Taylor_Hughes_Zarins|Taylor_Hughes_Zarins_Whiting_Jansen|Taylor_Hughes_Zarins_Whiting_Jansen_wo_dt|Taylor_Hughes_Zarins_scaled|Taylor_Hughes_Zarins_scaled_wo_dt|Taylor_Hughes_Zarins_wo_dt) "Definition of tau_M and Tau_C" // CHARELELENGTH_U (default: streamlength) (choices: root_of_volume|streamlength|volume_equivalent_diameter) "Characteristic element length for tau_Mu" // CHARELELENGTH_PC (default: volume_equivalent_diameter) (choices: root_of_volume|streamlength|volume_equivalent_diameter) "Characteristic element length for tau_Mp/tau_C" // EVALUATION_TAU (default: element_center) (choices: element_center|integration_point) "Location where tau is evaluated" // EVALUATION_MAT (default: element_center) (choices: element_center|integration_point) "Location where material law is evaluated" // LOMA_CONTI_SUPG (default: false) "Flag to (de)activate SUPG stabilization in loma continuity equation." // LOMA_CONTI_CROSS_STRESS (default: no_cross) (choices: cross_rhs|no_cross|yes_cross) "Flag to (de)activate cross-stress term loma continuity equation-> residual-based VMM." // LOMA_CONTI_REYNOLDS_STRESS (default: no_reynolds) (choices: no_reynolds|reynolds_rhs|yes_reynolds) "Flag to (de)activate Reynolds-stress term loma continuity equation-> residual-based VMM." .. _SECfluiddynamic_residual-basedstabilization: FLUID DYNAMIC/RESIDUAL-BASED STABILIZATION ------------------------------------------ FLUID DYNAMIC/RESIDUAL-BASED STABILIZATION :: // FLUID DYNAMIC/RESIDUAL-BASED STABILIZATION: // STABTYPE (default: residual_based) (choices: edge_based|no_stabilization|pressure_projection|residual_based) "Apply (un)stabilized fluid formulation. No stabilization is only possible for inf-sup stable elements.Use a residual-based stabilization or, more generally, a stabilization based on the concept of the residual-based variational multiscale method... Expecting additional input. Use an edge-based stabilization, especially for XFEM. Alternative: Element/cell based polynomial pressure projection, see Dohrmann/Bochev 2004, IJNMF" // INCONSISTENT (default: false) "residual based without second derivatives (i.e. only consistent for tau->0, but faster)" // Reconstruct_Sec_Der (default: false) "residual computed with a reconstruction of the second derivatives via projection or superconvergent patch recovery" // TDS (default: quasistatic) (choices: quasistatic|time_dependent) "Flag to allow time dependency of subscales for residual-based stabilization." // TRANSIENT (default: no_transient) (choices: no_transient|transient_complete|yes_transient) "Specify how to treat the transient term. Use transient term (recommended for time dependent subscales) or use transient term including a linearisation of 1/tau" // PSPG (default: true) "Flag to (de)activate PSPG stabilization." // SUPG (default: true) "Flag to (de)activate SUPG stabilization." // GRAD_DIV (default: true) "Flag to (de)activate grad-div term." // VSTAB (default: no_vstab) (choices: no_vstab|vstab_gls|vstab_gls_rhs|vstab_usfem|vstab_usfem_rhs) "Flag to (de)activate viscous term in residual-based stabilization. Options: No viscous term in stabilization, or, Viscous stabilization of GLS type, or, Viscous stabilization of GLS type, included only on the right hand side, or, Viscous stabilization of USFEM type, or, Viscous stabilization of USFEM type, included only on the right hand side" // RSTAB (default: no_rstab) (choices: no_rstab|rstab_gls|rstab_usfem) "Flag to (de)activate reactive term in residual-based stabilization." // CROSS-STRESS (default: no_cross) (choices: cross_rhs|no_cross|yes_cross) "Flag to (de)activate cross-stress term -> residual-based VMM. Options:No cross-stress term, or,Include the cross-stress term with a linearization of the convective part, or, Include cross-stress term, but only explicitly on right hand side" // REYNOLDS-STRESS (default: no_reynolds) (choices: no_reynolds|reynolds_rhs|yes_reynolds) "Flag to (de)activate Reynolds-stress term -> residual-based VMM." // DEFINITION_TAU (default: Franca_Barrenechea_Valentin_Frey_Wall) (choices: Codina|Codina_convscaled|Codina_wo_dt|Franca_Barrenechea_Valentin_Frey_Wall|Franca_Barrenechea_Valentin_Frey_Wall_wo_dt|Franca_Madureira_Valentin_Badia_Codina|Franca_Madureira_Valentin_Badia_Codina_wo_dt|Hughes_Franca_Balestra_wo_dt|Shakib_Hughes_Codina|Shakib_Hughes_Codina_wo_dt|Taylor_Hughes_Zarins|Taylor_Hughes_Zarins_Whiting_Jansen|Taylor_Hughes_Zarins_Whiting_Jansen_wo_dt|Taylor_Hughes_Zarins_scaled|Taylor_Hughes_Zarins_scaled_wo_dt|Taylor_Hughes_Zarins_wo_dt) "Definition of tau" // CHARELELENGTH_U (default: streamlength) (choices: root_of_volume|streamlength|volume_equivalent_diameter) "Characteristic element length for tau_Mu" // CHARELELENGTH_PC (default: volume_equivalent_diameter) (choices: root_of_volume|streamlength|volume_equivalent_diameter) "Characteristic element length for tau_Mp/tau_C" // EVALUATION_TAU (default: element_center) (choices: element_center|integration_point) "Location where tau is evaluated" // EVALUATION_MAT (default: element_center) (choices: element_center|integration_point) "Location where material law is evaluated" // LOMA_CONTI_SUPG (default: false) "Flag to (de)activate SUPG stabilization in loma continuity equation." // LOMA_CONTI_CROSS_STRESS (default: no_cross) (choices: cross_rhs|no_cross|yes_cross) "Flag to (de)activate cross-stress term loma continuity equation-> residual-based VMM." // LOMA_CONTI_REYNOLDS_STRESS (default: no_reynolds) (choices: no_reynolds|reynolds_rhs|yes_reynolds) "Flag to (de)activate Reynolds-stress term loma continuity equation-> residual-based VMM." .. _SECfluiddynamic_subgridviscosity: FLUID DYNAMIC/SUBGRID VISCOSITY ------------------------------- FLUID DYNAMIC/SUBGRID VISCOSITY :: // FLUID DYNAMIC/SUBGRID VISCOSITY: // C_SMAGORINSKY (default: 0) "Constant for the Smagorinsky model. Something between 0.1 to 0.24. Vreman constant if the constant vreman model is applied (something between 0.07 and 0.01)." // C_YOSHIZAWA (default: -1) "Constant for the compressible Smagorinsky model: isotropic part of subgrid-stress tensor. About 0.09 or 0.0066. Ci will not be squared!" // C_SMAGORINSKY_AVERAGED (default: false) "Flag to (de)activate averaged Smagorinksy constant" // C_INCLUDE_CI (default: false) "Flag to (de)inclusion of Yoshizawa model" // CHANNEL_L_TAU (default: 0) "Used for normalisation of the wall normal distance in the Van Driest Damping function. May be taken from the output of the apply_mesh_stretching.pl preprocessing script." // C_TURBPRANDTL (default: 1) "(Constant) turbulent Prandtl number for the Smagorinsky model in scalar transport." // FILTER_WIDTH (default: CubeRootVol) (choices: CubeRootVol|Direction_dependent|Minimum_length) "The Vreman model requires a filter width." .. _SECfluiddynamic_timeadaptivity: FLUID DYNAMIC/TIMEADAPTIVITY ---------------------------- FLUID DYNAMIC/TIMEADAPTIVITY :: // FLUID DYNAMIC/TIMEADAPTIVITY: // ADAPTIVE_TIME_STEP_ESTIMATOR (default: none) (choices: cfl_number|none|only_print_cfl_number) "Method used to determine adaptive time step size." // CFL_NUMBER (default: -1) "CFL number for adaptive time step" // FREEZE_ADAPTIVE_DT_AT (default: 1000000) "keep time step constant after this step, otherwise turbulence statistics sampling is not consistent" // ADAPTIVE_DT_INC (default: 0.8) "Increment of whole step for adaptive dt via CFL" .. _SECfluiddynamic_turbulencemodel: FLUID DYNAMIC/TURBULENCE MODEL ------------------------------ FLUID DYNAMIC/TURBULENCE MODEL :: // FLUID DYNAMIC/TURBULENCE MODEL: // TURBULENCE_APPROACH (default: DNS_OR_RESVMM_LES) (choices: CLASSICAL_LES|DNS_OR_RESVMM_LES) "Try to solve flow as an underresolved DNS. Mind that your stabilisation already acts as a kind of turbulence model! Perform a classical Large Eddy Simulation adding addititional turbulent viscosity. This may be based on various physical models.)" // PHYSICAL_MODEL (default: no_model) (choices: Dynamic_Smagorinsky|Dynamic_Vreman|Multifractal_Subgrid_Scales|Smagorinsky|Smagorinsky_with_van_Driest_damping|Vreman|no_model) "If classical LES is our turbulence approach, this is a contradiction and should cause a FOUR_C_THROW. Classical constant coefficient Smagorinsky model. Be careful if you have a wall bounded flow domain! Use an exponential damping function for the turbulent viscosity close to the wall. This is only implemented for a channel geometry of height 2 in y direction. The viscous lengthscale l_tau is required as additional input. The solution is filtered and by comparison of the filtered velocity field with the real solution, the Smagorinsky constant is estimated in each step --- mind that this procedure includes an averaging in the xz plane, hence this implementation will only work for a channel flow. Multifractal Subgrid-Scale Modeling based on the work of burton. Vremans constant model. Dynamic Vreman model according to You and Moin (2007)" // FSSUGRVISC (default: No) (choices: No|Smagorinsky_all|Smagorinsky_small) "fine-scale subgrid viscosity" // SAMPLING_START (default: 10000000) "Time step after when sampling shall be started" // SAMPLING_STOP (default: 1) "Time step when sampling shall be stopped" // DUMPING_PERIOD (default: 1) "Period of time steps after which statistical data shall be dumped" // SUBGRID_DISSIPATION (default: false) "Flag to (de)activate estimation of subgrid-scale dissipation (only for seclected flows)." // OUTMEAN (default: false) "Flag to (de)activate averaged paraview output" // TURBMODEL_LS (default: true) "Flag to (de)activate turbulence model in level-set equation" // CANONICAL_FLOW (default: no) (choices: backward_facing_step|backward_facing_step2|blood_fda_flow|bubbly_channel_flow|channel_flow_of_height_2|combust_oracles|decaying_homogeneous_isotropic_turbulence|forced_homogeneous_isotropic_turbulence|lid_driven_cavity|loma_backward_facing_step|loma_channel_flow_of_height_2|loma_lid_driven_cavity|no|periodic_hill|rotating_circular_cylinder_nurbs|rotating_circular_cylinder_nurbs_scatra|scatra_channel_flow_of_height_2|scatra_forced_homogeneous_isotropic_turbulence|square_cylinder|square_cylinder_nurbs|taylor_green_vortex|time_averaging) "Sampling is different for different canonical flows - so specify what kind of flow you've got no: The flow is not further specified, so spatial averaging and hence the standard sampling procedure is not possible time_averaging: The flow is not further specified, but time averaging of velocity and pressure field is performed channel_flow_of_height_2: For this flow, all statistical data could be averaged in the homogeneous planes - it is essentially a statistically one dimensional flow. lid_driven_cavity: For this flow, all statistical data are evaluated on the center lines of the xy-midplane, averaged only over time. backward_facing_step: For this flow, statistical data are evaluated on various lines, averaged over time and z. square_cylinder: For this flow, statistical data are evaluated on various lines of the xy-midplane, averaged only over time. square_cylinder_nurbs: For this flow, statistical data are evaluated on various lines of the xy-midplane, averaged over time and eventually in one home.direction. rotating_circular_cylinder_nurbs: For this flow, statistical data is computed in concentric surfaces and averaged. in time and in one home. direction rotating_circular_cylinder_nurbs_scatra: For this flow with mass transport, statistical data is computed in concentric surfaces and averaged. in time and in one home. direction loma_channel_flow_of_height_2: For this low-Mach-number flow, all statistical data could be averaged in the homogeneous planes - it is essentially a statistically one dimensional flow. loma_lid_driven_cavity: For this low-Mach-number flow, all statistical data are evaluated on the center lines of the xy-midplane, averaged only over time. loma_backward_facing_step: For this low-Mach-number flow, statistical data are evaluated on various lines, averaged over time and z. combust_oracles: ORACLES test rig for turbulent premixed combustion. bubbly_channel_flow: Turbulent two-phase flow: bubbly channel flow, statistical data are averaged in homogeneous planes and over time. scatra_channel_flow_of_height_2: For this flow, all statistical data could be averaged in the homogeneous planes - it is essentially a statistically one dimensional flow. decaying_homogeneous_isotropic_turbulence: For this flow, all statistical data could be averaged in the in all homogeneous directions - it is essentially a statistically zero dimensional flow. forced_homogeneous_isotropic_turbulence: For this flow, all statistical data could be averaged in the in all homogeneous directions - it is essentially a statistically zero dimensional flow. scatra_forced_homogeneous_isotropic_turbulence: For this flow, all statistical data could be averaged in the in all homogeneous directions - it is essentially a statistically zero dimensional flow. taylor_green_vortex: For this flow, dissipation rate could be averaged in the in all homogeneous directions - it is essentially a statistically zero dimensional flow. periodic_hill: For this flow, statistical data is evaluated on various lines, averaged over time and z. blood_fda_flow: For this flow, statistical data is evaluated on various planes. backward_facing_step2: For this flow, statistical data is evaluated on various planes." // HOMDIR (default: not_specified) (choices: not_specified|x|xy|xyz|xz|y|yz|z) "Specify the homogeneous direction(s) of a flow. not_specified: no homogeneous directions available, averaging is restricted to time averaging x: average along x-direction y: average along y-direction z: average along z-direction xy: Wall normal direction is z, average in x and y direction xz: Wall normal direction is y, average in x and z direction yz: Wall normal direction is x, average in y and z direction xyz: Averaging in all directions" // CHAN_AMPL_INIT_DIST (default: 0.1) "Max. amplitude of the random disturbance in percent of the initial value in mean flow direction." // FORCING_TYPE (default: linear_compensation_from_intermediate_spectrum) "forcing strategy" // CHA_NUMSUBDIVISIONS (default: 5) "Number of homogeneous sampling planes in element" // FORCING_TIME_STEPS (default: 0) "Number of time steps during which forcing is applied" // THRESHOLD_WAVENUMBER (default: 0) "Forcing is only applied to wave numbers lower or equal than the given threshold wave number." // POWER_INPUT (default: 0) "power of forcing" // SCALAR_FORCING (default: no) (choices: isotropic|mean_scalar_gradient|no) "no: Do not force the scalar field isotropic: Force scalar field isotropically such as the fluid field. mean_scalar_gradient: Force scalar field by imposed mean-scalar gradient." // MEAN_SCALAR_GRADIENT (default: 0) "Value of imposed mean-scalar gradient to force scalar field." // EXCLUDE_XFEM (default: false) "Flag to (de)activate XFEM dofs in calculation of fine-scale velocity." .. _SECfluiddynamic_turbulentinflow: FLUID DYNAMIC/TURBULENT INFLOW ------------------------------ FLUID DYNAMIC/TURBULENT INFLOW :: // FLUID DYNAMIC/TURBULENT INFLOW: // TURBULENTINFLOW (default: false) "Flag to (de)activate potential separate turbulent inflow section" // INITIALINFLOWFIELD (default: zero_field) (choices: disturbed_field_from_function|field_by_function|zero_field) "Initial field for inflow section" // INFLOWFUNC (default: -1) "Function number for initial flow field in inflow section" // INFLOW_INIT_DIST (default: 0.1) "Max. amplitude of the random disturbance in percent of the initial value in mean flow direction." // NUMINFLOWSTEP (default: 1) "Total number of time steps for development of turbulent flow" // CANONICAL_INFLOW (default: no) (choices: channel_flow_of_height_2|loma_channel_flow_of_height_2|no|scatra_channel_flow_of_height_2|time_averaging) "Sampling is different for different canonical flows --- so specify what kind of flow you've got no: The flow is not further specified, so spatial averaging and hence the standard sampling procedure is not possible. time_averaging: The flow is not further specified, but time averaging of velocity and pressure field is performed. channel_flow_of_height_2: For this flow, all statistical data could be averaged in the homogeneous planes --- it is essentially a statistically one dimensional flow. loma_channel_flow_of_height_2: For this low-Mach-number flow, all statistical data could be averaged in the homogeneous planes --- it is essentially a statistically one dimensional flow. scatra_channel_flow_of_height_2: For this flow, all statistical data could be averaged in the homogeneous planes --- it is essentially a statistically one dimensional flow." // INFLOW_CHA_SIDE (default: 0) "Most right side of inflow channel. Necessary to define sampling domain." // INFLOW_HOMDIR (default: not_specified) (choices: not_specified|x|xy|xz|y|yz|z) "Specify the homogeneous direction(s) of a flow not_specified: no homogeneous directions available, averaging is restricted to time averaging. x: average along x-direction. y: average along y-direction. z: average along z-direction. xy: Wall normal direction is z, average in x and y direction. xz: Wall normal direction is y, average in x and z direction (standard case). yz: Wall normal direction is x, average in y and z direction." // INFLOW_SAMPLING_START (default: 10000000) "Time step after when sampling shall be started" // INFLOW_SAMPLING_STOP (default: 1) "Time step when sampling shall be stopped" // INFLOW_DUMPING_PERIOD (default: 1) "Period of time steps after which statistical data shall be dumped" .. _SECfluiddynamic_wallmodel: FLUID DYNAMIC/WALL MODEL ------------------------ FLUID DYNAMIC/WALL MODEL :: // FLUID DYNAMIC/WALL MODEL: // X_WALL (default: false) "Flag to switch on the xwall model" // Tauw_Type (default: constant) (choices: between_steps|constant) "constant: Use the constant wall shear stress given in the input file for the whole simulation. between_steps: Calculate wall shear stress in between time steps." // Tauw_Calc_Type (default: residual) (choices: gradient|gradient_to_residual|residual) "residual: Residual (force) divided by area. gradient: Gradient via shape functions and nodal values. gradient_to_residual: First gradient, then residual." // Switch_Step (default: -1) "Switch from gradient to residual based tauw." // Projection (default: No) (choices: No|l2projectionwithcontinuityconstraint|onlyl2projection) "Flag to switch projection of the enriched dofs after updating tauw, alternatively with or without continuity constraint." // C_Tauw (default: 1) "Constant wall shear stress for Spalding's law, if applicable" // Min_Tauw (default: 2e-09) "Minimum wall shear stress preventing system to become singular" // Inc_Tauw (default: 1) "Increment of Tauw of full step, between 0.0 and 1.0" // Blending_Type (default: none) (choices: none|ramp_function) "Methods for blending the enrichment space. none: No ramp function, does not converge! ramp_function: Enrichment is multiplied with linear ramp function resulting in zero enrichment at the interface." // GP_Wall_Normal (default: 3) "Gauss points in wall normal direction" // GP_Wall_Normal_Off_Wall (default: 3) "Gauss points in wall normal direction, off-wall elements" // GP_Wall_Parallel (default: 3) "Gauss points in wall parallel direction" // Treat_Tauw_on_Dirichlet_Inflow (default: false) "Flag to treat residual on Dirichlet inflow nodes for calculation of wall shear stress" // PROJECTION_SOLVER (default: -1) "Set solver number for l2-projection." .. _SECfluidgeometry: FLUID GEOMETRY -------------- FLUID GEOMETRY :: // FLUID GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECfpsidynamic: FPSI DYNAMIC ------------ FPSI DYNAMIC :: // FPSI DYNAMIC: // COUPALGO (default: fpsi_monolithic_plain) "Iteration Scheme over the fields" // SHAPEDERIVATIVES (default: false) "Include linearization with respect to mesh movement in Navier Stokes equation. Supported in monolithic FPSI for now." // USESHAPEDERIVATIVES (default: false) "Add linearization with respect to mesh movement in Navier Stokes equation to stiffness matrix. Supported in monolithic FPSI for now." // PARTITIONED (default: RobinNeumann) "Coupling strategies for partitioned FPSI solvers." // SECONDORDER (default: false) "Second order coupling at the interface." // RESTOL (default: 1e-8 1e-8 1e-8 1e-8 1e-8 1e-8) "Tolerances for single fields in the residual norm for the Newton iteration. For NORM_RESF != Abs_sys_split only the first value is used for all fields. Order of fields: porofluidvelocity, porofluidpressure, porostructure, fluidvelocity, fluidpressure, ale" // INCTOL (default: 1e-8 1e-8 1e-8 1e-8 1e-8 1e-8) "Tolerance in the increment norm for the Newton iteration. For NORM_INC != \\*_split only the first value is used for all fields. Order of fields: porofluidvelocity, porofluidpressure, porostructure, fluidvelocity, fluidpressure, ale" // NORM_INC (default: Abs) (choices: Abs|Abs_sys_split|Rel_sys) "Type of norm for primary variables convergence check. Abs: absolute values, Abs_sys_split: absolute values with correction of systemsize for every field separate, Rel_sys: relative values with correction of systemsize." // NORM_RESF (default: Abs) (choices: Abs|Abs_sys_split|Rel_sys) "Type of norm for primary variables convergence check. Abs: absolute values, Abs_sys_split: absolute values with correction of systemsize for every field separate, Rel_sys: relative values with correction of systemsize." // NORMCOMBI_RESFINC (default: And) (choices: And|Or) "binary operator to combine primary variables and residual force values" // LineSearch (default: false) "adapt increment in case of non-monotonic residual convergence or residual oscillations" // FDCheck (default: false) "perform FPSIFDCheck() finite difference check" // LINEAR_SOLVER (default: -1) "number of linear solver used for FPSI problems" // ITEMIN (default: 1) "minimal number of iterations over fields" // NUMSTEP (default: 200) "Total number of Timesteps" // ITEMAX (default: 100) "Maximum number of iterations over fields" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // FDCheck_row (default: 0) "print row value during fd_check" // FDCheck_column (default: 0) "print column value during fd_check" // TIMESTEP (default: 0.1) "Time increment dt" // MAXTIME (default: 1000) "Total simulation time" // CONVTOL (default: 1e-06) "Tolerance for iteration over fields" // ALPHABJ (default: 1) "Beavers-Joseph-Coefficient for Slip-Boundary-Condition at Fluid-Porous-Interface (0.1-4)" .. _SECfs3idynamic: FS3I DYNAMIC ------------ FS3I DYNAMIC :: // FS3I DYNAMIC: // TIMESTEP (default: 0.1) "Time increment dt" // NUMSTEP (default: 20) "Total number of time steps" // MAXTIME (default: 1000) "Total simulation time" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // SCATRA_SOLVERTYPE (default: nonlinear) (choices: linear|nonlinear) "type of scalar transport solver" // INF_PERM (default: true) "Flag for infinite permeability" // CONSTHERMPRESS (default: Yes) (choices: No_energy|No_mass|Yes) "treatment of thermodynamic pressure in time" // COUPLED_LINEAR_SOLVER (default: -1) "number of linear solver used for fs3i problem" // LINEAR_SOLVER1 (default: -1) "number of linear solver used for fluid problem" // LINEAR_SOLVER2 (default: -1) "number of linear solver used for structural problem" // STRUCTSCAL_CONVFORM (default: conservative) (choices: conservative|convective) "form of convective term of structure scalar" // STRUCTSCAL_INITIALFIELD (default: zero_field) (choices: field_by_function|zero_field) "Initial Field for structure scalar transport problem" // STRUCTSCAL_INITFUNCNO (default: -1) "function number for structure scalar transport initial field" // STRUCTSCAL_FIELDCOUPLING (default: volume_matching) (choices: volume_matching|volume_nonmatching) "Type of coupling strategy between structure and structure-scalar field" // FLUIDSCAL_FIELDCOUPLING (default: volume_matching) (choices: volume_matching|volume_nonmatching) "Type of coupling strategy between fluid and fluid-scalar field" // FLUIDSCAL_SCATRATYPE (default: ConvectionDiffusion) (choices: Advanced_Reaction|Chemo_Reac|Chemotaxis|ConvectionDiffusion|Loma|Undefined) "Type of scalar transport problem" // RESTART_FROM_PART_FSI (default: false) "restart from partitioned fsi problem (e.g. from prestress calculations) instead of fs3i" .. _SECfs3idynamic_partitioned: FS3I DYNAMIC/PARTITIONED ------------------------ FS3I DYNAMIC/PARTITIONED :: // FS3I DYNAMIC/PARTITIONED: // COUPALGO (default: fs3i_IterStagg) "Coupling strategies for FS3I solvers" // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteration within partitioned FS3I" // ITEMAX (default: 10) "Maximum number of outer iterations" .. _SECfs3idynamic_structurescalarstabilization: FS3I DYNAMIC/STRUCTURE SCALAR STABILIZATION ------------------------------------------- FS3I DYNAMIC/STRUCTURE SCALAR STABILIZATION :: // FS3I DYNAMIC/STRUCTURE SCALAR STABILIZATION: // STABTYPE (default: SUPG) (choices: GLS|SUPG|USFEM|centered|no_stabilization|upwind) "Type of stabilization (if any). No stabilization is only reasonable for low-Peclet-number." // SUGRVEL (default: false) "potential incorporation of subgrid-scale velocity" // ASSUGRDIFF (default: false) "potential incorporation of all-scale subgrid diffusivity (a.k.a. discontinuity-capturing) term" // DEFINITION_TAU (default: Franca_Valentin) (choices: Codina|Codina_wo_dt|Exact_1D|Franca_Madureira_Valentin|Franca_Madureira_Valentin_wo_dt|Franca_Valentin|Franca_Valentin_wo_dt|Numerical_Value|Shakib_Hughes_Codina|Shakib_Hughes_Codina_wo_dt|Taylor_Hughes_Zarins|Taylor_Hughes_Zarins_wo_dt|Zero) "Definition of tau" // CHARELELENGTH (default: streamlength) "Characteristic element length for tau" // DEFINITION_ASSGD (default: artificial_linear) (choices: Almeida_Silva_97_nonlinear|Codina_nonlinear|Hughes_etal_86_nonlinear|Tezduyar_Park_86_nonlinear|Tezduyar_Park_86_nonlinear_wo_phizero|YZbeta_nonlinear|artificial_linear|artificial_linear_reinit|doCarmo_Galeao_91_nonlinear) "Definition of (all-scale) subgrid diffusivity" // EVALUATION_TAU (default: element_center) (choices: element_center|integration_point) "Location where tau is evaluated" // EVALUATION_MAT (default: element_center) (choices: element_center|integration_point) "Location where material law is evaluated" // CONSISTENCY (default: no) (choices: L2_projection_lumped|no) "improvement of consistency for stabilization" // TAU_VALUE (default: 0) "Numerical value for tau for stabilization" .. _SECfsidynamic: FSI DYNAMIC ----------- FSI DYNAMIC :: // FSI DYNAMIC: // COUPALGO (default: iter_stagg_AITKEN_rel_param) (choices: basic_sequ_stagg|iter_fluidfluid_monolithicfluidsplit|iter_fluidfluid_monolithicfluidsplit_nonox|iter_fluidfluid_monolithicstructuresplit|iter_fluidfluid_monolithicstructuresplit_nonox|iter_monolithicfluidsplit|iter_monolithicstructuresplit|iter_mortar_monolithicfluidsplit|iter_mortar_monolithicfluidsplit_saddlepoint|iter_mortar_monolithicstructuresplit|iter_sliding_monolithicfluidsplit|iter_sliding_monolithicstructuresplit|iter_stagg_AITKEN_rel_param|iter_stagg_MFNK_FD|iter_stagg_MFNK_FSI|iter_stagg_MPE|iter_stagg_NLCG|iter_stagg_RRE|iter_stagg_fixed_rel_param|iter_stagg_steep_desc|iter_xfem_monolithic) "Iteration Scheme over the fields" // DEBUGOUTPUT (default: false) "Output of unconverged interface values during FSI iteration. There will be a new control file for each time step. This might be helpful to understand the coupling iteration." // MATCHGRID_FLUIDALE (default: true) "is matching grid (fluid-ale)" // MATCHGRID_STRUCTALE (default: true) "is matching grid (structure-ale)" // MATCHALL (default: true) "is matching grid (fluid-ale) and is full fluid-ale (without euler part)" // MAXTIME (default: 1000) "Total simulation time" // NUMSTEP (default: 200) "Total number of Timesteps" // RESTARTEVERY (default: 1) "Increment for writing restart" // RESTART_FROM_PART_FSI (default: false) "restart from partitioned fsi (e.g. from prestress calculations) instead of monolithic fsi" // SECONDORDER (default: false) "Second order displacement-velocity conversion at the interface." // SLIDEALEPROJ (default: None) (choices: Curr|None|Ref|RotZ|RotZSphere) "Projection method to use for sliding FSI." // TIMESTEP (default: 0.1) "Time increment dt" // RESULTSEVERY (default: 1) "Increment for writing solution" // VERBOSITY (default: full) (choices: full|low|medium|subproblem) "Verbosity of the FSI problem." .. _SECfsidynamic_constraint: FSI DYNAMIC/CONSTRAINT ---------------------- FSI DYNAMIC/CONSTRAINT :: // FSI DYNAMIC/CONSTRAINT: // PRECONDITIONER (default: Simple) (choices: Simple|Simplec) "preconditioner to use" // SIMPLEITER (default: 2) "Number of iterations for simple pc" // ALPHA (default: 0.8) "alpha parameter for simple pc" .. _SECfsidynamic_monolithicsolver: FSI DYNAMIC/MONOLITHIC SOLVER ----------------------------- FSI DYNAMIC/MONOLITHIC SOLVER :: // FSI DYNAMIC/MONOLITHIC SOLVER: // ADAPTIVEDIST (default: 0) "Required distance for adaptive convergence check in Newton-type FSI. This is the improvement we want to achieve in the linear extrapolation of the adaptive convergence check. Set to zero to avoid the adaptive check altogether." // BASETOL (default: 0.001) "Basic tolerance for adaptive convergence check in monolithic FSI. This tolerance will be used for the linear solve of the FSI block system. The linear convergence test will always use the relative residual norm (AZ_r0). Not to be confused with the Newton tolerance (CONVTOL) that applies to the nonlinear convergence test using a absolute residual norm." // CONVTOL (default: 1e-06) "Nonlinear tolerance for lung/constraint/fluid-fluid FSI" // ENERGYFILE (default: false) "Write artificial interface energy due to temporal discretization to file" // FSIAMGANALYZE (default: false) "run analysis on fsiamg multigrid scheme" // INFNORMSCALING (default: true) "Scale Blocks with row infnorm?" // ITEMAX (default: 100) "Maximum allowed number of nonlinear iterations" // KRYLOV_ITEMAX (default: 1000) "Max Iterations for linear solver." // KRYLOV_SIZE (default: 50) "Size of Krylov Subspace." // LINEARBLOCKSOLVER (default: PreconditionedKrylov) (choices: LinalgSolver|PreconditionedKrylov) "Linear block preconditioner for block system in monolithic FSI." // LINEAR_SOLVER (default: -1) "Number of SOLVER block describing the linear solver and preconditioner" // NORM_INC (default: Rel) (choices: Abs|Mix|Rel) "type of norm for primary variables convergence check" // NORM_RESF (default: Rel) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // NORMCOMBI_RESFINC (default: And) (choices: And) "binary operator to combine primary variables and residual force values" // PRECONDREUSE (default: 0) "Number of iterations in one time step reusing the preconditioner before rebuilding it" // REBUILDPRECEVERYSTEP (default: true) "Enforce rebuilding the preconditioner at the beginning of every time step" // SHAPEDERIVATIVES (default: false) "Include linearization with respect to mesh movement in Navier Stokes equation." // SYMMETRICPRECOND (default: false) "Symmetric block GS preconditioner or ordinary GS" // ALEPCOMEGA (default: 1.0 1.0 1.0 1.0) "Relaxation factor for Richardson iteration on ale block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // ALEPCITER (default: 1 1 1 1) "Number of Richardson iterations on ale block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // FLUIDPCOMEGA (default: 1.0 1.0 1.0 1.0) "Relaxation factor for Richardson iteration on fluid block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // FLUIDPCITER (default: 1 1 1 1) "Number of Richardson iterations on fluid block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // STRUCTPCOMEGA (default: 1.0 1.0 1.0 1.0) "Relaxation factor for Richardson iteration on structural block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // STRUCTPCITER (default: 1 1 1 1) "Number of Richardson iterations on structural block in MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // PCOMEGA (default: 1.0 1.0 1.0) "Relaxation factor for Richardson iteration on whole MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // PCITER (default: 1 1 1) "Number of Richardson iterations on whole MFSI block preconditioner FSIAMG: each number belongs to a level PreconditiondKrylov: only first number is used for finest level" // BLOCKSMOOTHER (default: BGS BGS BGS) "Type of block smoother, can be BGS or Schur" // SCHUROMEGA (default: 0.001 0.01 0.1) "Damping factor for Schur complement construction" // TOL_DIS_RES_L2 (default: 1e-06) "Absolute tolerance for structure displacement residual in L2-norm" // TOL_DIS_RES_INF (default: 1e-06) "Absolute tolerance for structure displacement residual in Inf-norm" // TOL_DIS_INC_L2 (default: 1e-06) "Absolute tolerance for structure displacement increment in L2-norm" // TOL_DIS_INC_INF (default: 1e-06) "Absolute tolerance for structure displacement increment in Inf-norm" // TOL_FSI_RES_L2 (default: 1e-06) "Absolute tolerance for interface residual in L2-norm" // TOL_FSI_RES_INF (default: 1e-06) "Absolute tolerance for interface residual in Inf-norm" // TOL_FSI_INC_L2 (default: 1e-06) "Absolute tolerance for interface increment in L2-norm" // TOL_FSI_INC_INF (default: 1e-06) "Absolute tolerance for interface increment in Inf-norm" // TOL_PRE_RES_L2 (default: 1e-06) "Absolute tolerance for fluid pressure residual in L2-norm" // TOL_PRE_RES_INF (default: 1e-06) "Absolute tolerance for fluid pressure residual in Inf-norm" // TOL_PRE_INC_L2 (default: 1e-06) "Absolute tolerance for fluid pressure increment in L2-norm" // TOL_PRE_INC_INF (default: 1e-06) "Absolute tolerance for fluid pressure increment in Inf-norm" // TOL_VEL_RES_L2 (default: 1e-06) "Absolute tolerance for fluid velocity residual in L2-norm" // TOL_VEL_RES_INF (default: 1e-06) "Absolute tolerance for fluid velocity residual in Inf-norm" // TOL_VEL_INC_L2 (default: 1e-06) "Absolute tolerance for fluid velocity increment in L2-norm" // TOL_VEL_INC_INF (default: 1e-06) "Absolute tolerance for fluid velocity increment in Inf-norm" .. _SECfsidynamic_partitionedsolver: FSI DYNAMIC/PARTITIONED SOLVER ------------------------------ FSI DYNAMIC/PARTITIONED SOLVER :: // FSI DYNAMIC/PARTITIONED SOLVER: // BASETOL (default: 0.001) "Basic tolerance for adaptive convergence check in monolithic FSI. This tolerance will be used for the linear solve of the FSI block system. The linear convergence test will always use the relative residual norm (AZ_r0). Not to be confused with the Newton tolerance (CONVTOL) that applies to the nonlinear convergence test using a absolute residual norm." // CONVTOL (default: 1e-06) "Tolerance for iteration over fields in case of partitioned scheme" // COUPMETHOD (default: conforming) (choices: conforming|immersed|mortar) "Coupling Method mortar or conforming nodes at interface" // COUPVARIABLE (default: Displacement) (choices: Displacement|Force|Velocity) "Coupling variable at the interface" // DIVPROJECTION (default: false) "Project velocity into divergence-free subspace for partitioned fsi" // ITEMAX (default: 100) "Maximum number of iterations over fields" // MAXOMEGA (default: 0) "largest omega allowed for Aitken relaxation (0.0 means no constraint)" // MINOMEGA (default: -1) "smallest omega allowed for Aitken relaxation (default is -1.0)" // PARTITIONED (default: DirichletNeumann) (choices: DirichletNeumann|DirichletNeumannSlideALE|DirichletNeumannVolCoupl) "Coupling strategies for partitioned FSI solvers." // PREDICTOR (default: d(n)) (choices: d(n)|d(n)+dt*(1.5*v(n)-0.5*v(n-1))|d(n)+dt*v(n)|d(n)+dt*v(n)+0.5*dt^2*a(n)) "Predictor for interface displacements" // RELAX (default: 1) "fixed relaxation parameter for partitioned FSI solvers" .. _SECfsidynamic_timeadaptivity: FSI DYNAMIC/TIMEADAPTIVITY -------------------------- FSI DYNAMIC/TIMEADAPTIVITY :: // FSI DYNAMIC/TIMEADAPTIVITY: // ADAPTSTEPMAX (default: 5) "Maximum number of repetitions of one time step for adapting/reducing the time step size (>0)" // AUXINTEGRATORFLUID (default: AB2) (choices: AB2|ExplicitEuler|None) "Method for error estimation in the fluid field" // AVERAGINGDT (default: 0.3 0.7) "Averaging of time step sizes in case of increasing time step size. Parameters are ordered from most recent weight to the most historic one. Number of parameters determines the number of previous time steps that are involved in the averaging procedure." // DIVERCONT (default: stop) (choices: continue|halve_step|revert_dt|stop) "What to do if nonlinear solver does not converge?" // DTMAX (default: 0.1) "Limit maximally permitted time step size (>0)" // DTMIN (default: 0.0001) "Limit minimally allowed time step size (>0)" // LOCERRTOLFLUID (default: 0.001) "Tolerance for the norm of local velocity error" // NUMINCREASESTEPS (default: 0) "Number of consecutive steps that want to increase time step size before actually increasing it. Set 0 to deactivate this feature." // SAFETYFACTOR (default: 0.9) "This is a safety factor to scale theoretical optimal step size, should be lower than 1 and must be larger than 0" // SIZERATIOMAX (default: 2) "Limit maximally permitted change of time step size compared to previous size (>0)." // SIZERATIOMIN (default: 0.5) "Limit minimally permitted change of time step size compared to previous size (>0)." // TIMEADAPTON (default: false) "Activate or deactivate time step size adaptivity" .. _SECio: IO -- IO :: // IO: // OUTPUT_GMSH (default: false) // OUTPUT_ROT (default: false) // OUTPUT_SPRING (default: false) // OUTPUT_BIN (default: true) "Do you want to have binary output?" // OUTPUT_EVERY_ITER (default: false) "Do you desire structural displ. output every Newton iteration" // OEI_FILE_COUNTER (default: 0) "Add an output name affix by introducing a additional number" // ELEMENT_MAT_ID (default: false) "Output of the material id of each element" // STRUCT_ELE (default: true) "Output of element properties" // STRUCT_DISP (default: true) "Output of displacements" // STRUCT_STRESS (default: NO) (choices: 2PK|2pk|Cauchy|NO|No|YES|Yes|cauchy|no|yes) "Output of stress" // STRUCT_COUPLING_STRESS (default: NO) (choices: 2PK|2pk|Cauchy|NO|No|YES|Yes|cauchy|no|yes) // STRUCT_STRAIN (default: NO) (choices: EA|GL|LOG|NO|No|YES|Yes|ea|gl|log|no|yes) "Output of strains" // STRUCT_PLASTIC_STRAIN (default: NO) (choices: EA|GL|NO|No|YES|Yes|ea|gl|no|yes) // STRUCT_SURFACTANT (default: false) // STRUCT_JACOBIAN_MATLAB (default: false) // STRUCT_CONDITION_NUMBER (default: none) (choices: gmres_estimate|inf-norm|max_min_ev_ratio|none|one-norm) "Compute the condition number of the structural system matrix and write it to a text file." // FLUID_STRESS (default: false) // FLUID_WALL_SHEAR_STRESS (default: false) // FLUID_ELEDATA_EVERY_STEP (default: false) // FLUID_NODEDATA_FIRST_STEP (default: false) // THERM_TEMPERATURE (default: false) // THERM_HEATFLUX (default: NO) (choices: Current|Initial|NO|No|None|no) // THERM_TEMPGRAD (default: NO) (choices: Current|Initial|NO|No|None|no) // FILESTEPS (default: 1000) "Amount of timesteps written to a single result file" // STDOUTEVERY (default: 1) "Print to screen every n step" // WRITE_TO_SCREEN (default: true) "Write screen output" // WRITE_TO_FILE (default: false) "Write the output into a file" // WRITE_INITIAL_STATE (default: true) "Do you want to write output for initial state ?" // WRITE_FINAL_STATE (default: false) "Enforce to write output/restart data at the final state regardless of the other output/restart intervals" // PREFIX_GROUP_ID (default: false) "Put a : in front of every line" // LIMIT_OUTP_TO_PROC (default: -1) "Only the specified procs will write output" // VERBOSITY (default: Verbose) (choices: Debug|Minimal|Standard|Verbose|debug|minimal|standard|verbose) // RESTARTWALLTIMEINTERVAL (default: -1) "Enforce restart after this walltime interval (in seconds), smaller zero to disable" // RESTARTEVERY (default: -1) "write restart every RESTARTEVERY steps" .. _SECio_everyiteration: IO/EVERY ITERATION ------------------ IO/EVERY ITERATION :: // IO/EVERY ITERATION: // OUTPUT_EVERY_ITER (default: false) "Do you wish output every Newton iteration?" // RUN_NUMBER (default: -1) "Create a new folder for different runs of the same simulation. If equal -1, no folder is created." // STEP_NP_NUMBER (default: -1) "Give the number of the step (i.e. step_{n+1}) for which you want to write the debug output. If a negative step number is provided, all steps willbe written." // WRITE_OWNER_EACH_NEWTON_ITER (default: false) "If yes, the ownership of elements and nodes are written each Newton step, instead of only once per time/load step." .. _SECio_monitorstructuredbc: IO/MONITOR STRUCTURE DBC ------------------------ IO/MONITOR STRUCTURE DBC :: // IO/MONITOR STRUCTURE DBC: // INTERVAL_STEPS (default: -1) "write reaction force output every INTERVAL_STEPS steps" // PRECISION_FILE (default: 16) "precision for written file" // PRECISION_SCREEN (default: 5) "precision for written screen output" // FILE_TYPE (default: CSV) (choices: CSV|Csv|DATA|Data|csv|data) "type of written output file" // WRITE_HEADER (default: false) "write information about monitored boundary condition to output file" .. _SECio_runtimevtkoutput: IO/RUNTIME VTK OUTPUT --------------------- IO/RUNTIME VTK OUTPUT :: // IO/RUNTIME VTK OUTPUT: // INTERVAL_STEPS (default: -1) "write visualization output at runtime every INTERVAL_STEPS steps" // STEP_OFFSET (default: 0) "An offset added to the current step to shift the steps to be written." // OUTPUT_DATA_FORMAT (default: binary) "data format for written numeric data" // COMPRESSION_LEVEL (default: best_speed) "Specify the compression level of written vtk output." // TIMESTEP_RESERVE_DIGITS (default: 5) "Specify the maximum digits in the number of time steps that shall be written. This only affects the number of leading zeros in the output file names." // EVERY_ITERATION (default: false) "write output in every iteration of the nonlinear solver" // EVERY_ITERATION_VIRTUAL_TIME_INCREMENT (default: 1e-08) "Specify the virtual time increment that is added for each nonlinear output state" // EVERY_ITERATION_RESERVE_DIGITS (default: 4) "Specify the maximum digits in the number of iterations that shall be written. This only affects the number of leading zeros in the output file names." // OUTPUT_WRITER (default: vtu_per_rank) "Specify which output writer shall be used to write the visualization data to disk" .. _SECio_runtimevtkoutput_beams: IO/RUNTIME VTK OUTPUT/BEAMS --------------------------- IO/RUNTIME VTK OUTPUT/BEAMS :: // IO/RUNTIME VTK OUTPUT/BEAMS: // OUTPUT_BEAMS (default: false) "write special output for beam elements" // DISPLACEMENT (default: false) "write displacement output" // USE_ABSOLUTE_POSITIONS (default: true) "use absolute positions or initial positions for vtu geometry (i.e. point coordinates)" // INTERNAL_ENERGY_ELEMENT (default: false) "write internal (elastic) energy for each element" // KINETIC_ENERGY_ELEMENT (default: false) "write kinetic energy for each element" // TRIAD_VISUALIZATIONPOINT (default: false) "write triads at every visualization point" // STRAINS_GAUSSPOINT (default: false) "write material cross-section strains at the Gauss points" // STRAINS_CONTINUOUS (default: false) "write material cross-section strains at the visualization points" // MATERIAL_FORCES_GAUSSPOINT (default: false) "write material cross-section stresses at the Gauss points" // MATERIAL_FORCES_CONTINUOUS (default: false) "write material cross-section stresses at the visualization points" // SPATIAL_FORCES_GAUSSPOINT (default: false) "write material cross-section stresses at the Gauss points" // BEAMFILAMENTCONDITION (default: false) "write element filament numbers" // ORIENTATION_PARAMETER (default: false) "write element filament numbers" // RVE_CROSSSECTION_FORCES (default: false) "get sum of all internal forces of" // REF_LENGTH (default: false) "write reference length of all beams" // ELEMENT_GID (default: false) "write the 4C internal element GIDs" // ELEMENT_GHOSTING (default: false) "write which processors ghost the elements" // NUMBER_SUBSEGMENTS (default: 5) "Number of subsegments along a single beam element for visualization" .. _SECio_runtimevtkoutput_fluid: IO/RUNTIME VTK OUTPUT/FLUID --------------------------- IO/RUNTIME VTK OUTPUT/FLUID :: // IO/RUNTIME VTK OUTPUT/FLUID: // OUTPUT_FLUID (default: false) "write fluid output" // VELOCITY (default: false) "write velocity output" // PRESSURE (default: false) "write pressure output" // ACCELERATION (default: false) "write acceleration output" // DISPLACEMENT (default: false) "write displacement output" // GRIDVELOCITY (default: false) "write grid velocity output" // ELEMENT_OWNER (default: false) "write element owner" // ELEMENT_GID (default: false) "write 4C internal element GIDs" // NODE_GID (default: false) "write 4C internal node GIDs" .. _SECio_runtimevtkoutput_structure: IO/RUNTIME VTK OUTPUT/STRUCTURE ------------------------------- IO/RUNTIME VTK OUTPUT/STRUCTURE :: // IO/RUNTIME VTK OUTPUT/STRUCTURE: // OUTPUT_STRUCTURE (default: false) "write structure output" // DISPLACEMENT (default: false) "write displacement output" // VELOCITY (default: false) "write velocity output" // ACCELERATION (default: false) "write acceleration output" // ELEMENT_OWNER (default: false) "write element owner" // ELEMENT_GID (default: false) "write 4C internal element GIDs" // ELEMENT_GHOSTING (default: false) "write which processors ghost the elements" // NODE_GID (default: false) "write 4C internal node GIDs" // ELEMENT_MAT_ID (default: false) "Output of the material id of each element" // STRESS_STRAIN (default: false) "Write element stress and / or strain data. The type of stress / strain has to be selected in the --IO input section" // GAUSS_POINT_DATA_OUTPUT_TYPE (default: none) "Where to write gauss point data. (none, projected to nodes, projected to element center, raw at gauss points)" // OPTIONAL_QUANTITY (default: No) (choices: No|membranethickness) "Output of an optional quantity" .. _SECio_runtimevtpoutputstructure: IO/RUNTIME VTP OUTPUT STRUCTURE ------------------------------- IO/RUNTIME VTP OUTPUT STRUCTURE :: // IO/RUNTIME VTP OUTPUT STRUCTURE: // INTERVAL_STEPS (default: -1) "write VTP output at runtime every INTERVAL_STEPS steps" // STEP_OFFSET (default: 0) "An offset added to the current step to shift the steps to be written." // EVERY_ITERATION (default: false) "write output in every iteration of the nonlinear solver" // OWNER (default: false) "write owner of every point" // ORIENTATIONANDLENGTH (default: false) "write orientation at every point" // NUMBEROFBONDS (default: false) "write number of bonds of every point" // LINKINGFORCE (default: false) "write force acting in linker" .. _SEClevel-setcontrol: LEVEL-SET CONTROL ----------------- LEVEL-SET CONTROL :: // LEVEL-SET CONTROL: // NUMSTEP (default: 24) "Total number of time steps" // TIMESTEP (default: 0.1) "Time increment dt" // MAXTIME (default: 1000) "Total simulation time" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // CALCERROR (default: No) (choices: InitialField|No) "compute error compared to analytical solution" // EXTRACT_INTERFACE_VEL (default: false) "replace computed velocity at nodes of given distance of interface by approximated interface velocity" // NUM_CONVEL_LAYERS (default: -1) "number of layers around the interface which keep their computed convective velocity" .. _SEClevel-setcontrol_reinitialization: LEVEL-SET CONTROL/REINITIALIZATION ---------------------------------- LEVEL-SET CONTROL/REINITIALIZATION :: // LEVEL-SET CONTROL/REINITIALIZATION: // REINITIALIZATION (default: None) (choices: EllipticEq|None|Signed_Distance_Function|Sussman) "Type of reinitialization strategy for level set function" // REINIT_INITIAL (default: false) "Has level set field to be reinitialized before first time step?" // REINITINTERVAL (default: 1) "reinitialization interval" // REINITBAND (default: false) "reinitialization only within a band around the interface, or entire domain?" // REINITBANDWIDTH (default: 1) "level-set value defining band width for reinitialization" // NUMSTEPSREINIT (default: 1) "(maximal) number of pseudo-time steps" // TIMESTEPREINIT (default: 1) "pseudo-time step length (usually a * characteristic element length of discretization with a>0)" // THETAREINIT (default: 1) "theta for time discretization of reinitialization equation" // STABTYPEREINIT (default: SUPG) (choices: GLS|SUPG|USFEM|no_stabilization) "Type of stabilization (if any). No stabilization is only reasonable for low-Peclet-number." // DEFINITION_TAU_REINIT (default: Taylor_Hughes_Zarins) (choices: Codina|Codina_wo_dt|Exact_1D|Franca_Valentin|Franca_Valentin_wo_dt|Shakib_Hughes_Codina|Shakib_Hughes_Codina_wo_dt|Taylor_Hughes_Zarins|Taylor_Hughes_Zarins_wo_dt|Zero) "Definition of tau" // ARTDIFFREINIT (default: no) (choices: crosswind|isotropic|no) "potential incorporation of all-scale subgrid diffusivity (a.k.a. discontinuity-capturing) term" // DEFINITION_ARTDIFFREINIT (default: artificial_linear) (choices: Almeida_Silva_97_nonlinear|Codina_nonlinear|Hughes_etal_86_nonlinear|Tezduyar_Park_86_nonlinear|Tezduyar_Park_86_nonlinear_wo_phizero|YZbeta_nonlinear|artificial_linear|artificial_linear_reinit|doCarmo_Galeao_91_nonlinear) "Definition of (all-scale) subgrid diffusivity" // SMOOTHED_SIGN_TYPE (default: SussmanSmerekaOsher1994) (choices: NonSmoothed|PengEtAl1999|SussmanFatemi1999|SussmanSmerekaOsher1994) "sign function for reinitialization equation" // CHARELELENGTHREINIT (default: root_of_volume) (choices: root_of_volume|streamlength) "characteristic element length for sign function" // INTERFACE_THICKNESS (default: 1) "factor for interface thickness (multiplied by element length)" // VELREINIT (default: integration_point_based) (choices: integration_point_based|node_based) "evaluate velocity at integration point or compute node-based velocity" // LINEARIZATIONREINIT (default: newton) (choices: fixed_point|newton) "linearization scheme for nonlinear convective term of reinitialization equation" // CORRECTOR_STEP (default: true) "correction of interface position via volume constraint according to Sussman & Fatemi" // CONVTOL_REINIT (default: -1) "tolerance for convergence check according to Sussman et al. 1994 (turned off negative)" // REINITVOLCORRECTION (default: false) "volume correction after reinitialization" // PENALTY_PARA (default: -1) "penalty parameter for elliptic reinitialization" // DIMENSION (default: 3D) (choices: 2Dx|2Dy|2Dz|3D) "number of space dimensions for handling of quasi-2D problems with 3D elements" // PROJECTION (default: true) "use L2-projection for grad phi and related quantities" // PROJECTION_DIFF (default: 0) "use diffusive term for L2-projection" // LUMPING (default: false) "use lumped mass matrix for L2-projection" // DIFF_FUNC (default: hyperbolic) (choices: hyperbolic|hyperbolic_clipped_05|hyperbolic_clipped_1|hyperbolic_smoothed_positive) "function for diffusivity" .. _SEClomacontrol: LOMA CONTROL ------------ LOMA CONTROL :: // LOMA CONTROL: // MONOLITHIC (default: false) "monolithic solver" // NUMSTEP (default: 24) "Total number of time steps" // TIMESTEP (default: 0.1) "Time increment dt" // MAXTIME (default: 1000) "Total simulation time" // ITEMAX (default: 10) "Maximum number of outer iterations" // ITEMAX_BEFORE_SAMPLING (default: 1) "Maximum number of outer iterations before sampling (for turbulent flows only)" // CONVTOL (default: 1e-06) "Tolerance for convergence check" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // CONSTHERMPRESS (default: Yes) (choices: No_energy|No_mass|Yes) "treatment of thermodynamic pressure in time" // SGS_MATERIAL_UPDATE (default: false) "update material by adding subgrid-scale scalar field" // LINEAR_SOLVER (default: -1) "number of linear solver used for LOMA problem" .. _SEClubricationdynamic: LUBRICATION DYNAMIC ------------------- LUBRICATION DYNAMIC :: // LUBRICATION DYNAMIC: // MAXTIME (default: 1000) "Total simulation time" // NUMSTEP (default: 20) "Total number of time steps" // TIMESTEP (default: 0.1) "Time increment dt" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // CALCERROR (default: No) (choices: No|error_by_function) "compute error compared to analytical solution" // CALCERRORNO (default: -1) "function number for lubrication error computation" // VELOCITYFIELD (default: zero) (choices: EHL|function|zero) "type of velocity field used for lubrication problems" // VELFUNCNO (default: -1) "function number for lubrication velocity field" // HEIGHTFEILD (default: zero) (choices: EHL|function|zero) "type of height field used for lubrication problems" // HFUNCNO (default: -1) "function number for lubrication height field" // OUTMEAN (default: false) "Output of mean values for scalars and density" // OUTPUT_GMSH (default: false) "Do you want to write Gmsh postprocessing files?" // MATLAB_STATE_OUTPUT (default: false) "Do you want to write the state solution to Matlab file?" // LINEAR_SOLVER (default: -1) "number of linear solver used for the Lubrication problem" // ITEMAX (default: 10) "max. number of nonlin. iterations" // ABSTOLRES (default: 1e-14) "Absolute tolerance for deciding if residual of nonlinear problem is already zero" // CONVTOL (default: 1e-13) "Tolerance for convergence check" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // NORM_PRE (default: Abs) (choices: Abs|Mix|Rel) "type of norm for temperature convergence check" // NORM_RESF (default: Abs) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // ITERNORM (default: L2) (choices: Inf|L1|L2|Rms) "type of norm to be applied to residuals" // TOLPRE (default: 1e-06) "tolerance in the temperature norm of the Newton iteration" // TOLRES (default: 1e-06) "tolerance in the residual norm for the Newton iteration" // PENALTY_CAVITATION (default: 0) "penalty parameter for regularized cavitation" // GAP_OFFSET (default: 0) "Additional offset to the fluid gap" // ROUGHNESS_STD_DEVIATION (default: 0) "standard deviation of surface roughness" // MODIFIED_REYNOLDS_EQU (default: false) "the lubrication problem will use the modified reynolds equ. in order to consider surface roughness" // ADD_SQUEEZE_TERM (default: false) "the squeeze term will also be considered in the Reynolds Equation" // PURE_LUB (default: false) "the problem is pure lubrication" .. _SEClubricationgeometry: LUBRICATION GEOMETRY -------------------- LUBRICATION GEOMETRY :: // LUBRICATION GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECmeshpartitioning: MESH PARTITIONING ----------------- MESH PARTITIONING :: // MESH PARTITIONING: // METHOD (default: hypergraph) "Type of rebalance/partition algorithm to be used for decomposing the entire mesh into subdomains for parallel computing." // IMBALANCE_TOL (default: 1.1) "Tolerance for relative imbalance of subdomain sizes for graph partitioning of unstructured meshes read from input files." // MIN_ELE_PER_PROC (default: 0) "This parameter defines the minimum number of elements to be assigned to any MPI rank during redistribution. Use 0 to not interfere with the minimal size of a subdomain." .. _SECmor: MOR --- MOR :: // MOR: // POD_MATRIX (default: none) "filename of file containing projection matrix" .. _SECmortarcoupling: MORTAR COUPLING --------------- MORTAR COUPLING :: // MORTAR COUPLING: // LM_SHAPEFCN (default: Dual) (choices: Dual|PetrovGalerkin|Standard|dual|petrovgalerkin|pg|standard|std) "Type of employed set of shape functions" // SEARCH_ALGORITHM (default: BinaryTree) (choices: BinaryTree|Binarytree|BruteForce|BruteForceEleBased|binarytree|bruteforce|bruteforceelebased) "Type of contact search" // BINARYTREE_UPDATETYPE (default: BottomUp) (choices: BottomUp|TopDown) "Type of binary tree update, which is either a bottom up or a top down approach." // SEARCH_PARAM (default: 0.3) "Radius / Bounding volume inflation for contact search" // SEARCH_USE_AUX_POS (default: true) "If chosen auxiliary position is used for computing dops" // LM_QUAD (default: undefined) (choices: const|lin|linear|piecewiselinear|pwlin|quad|quadratic|undefined) "Type of LM interpolation for quadratic FE" // CROSSPOINTS (default: false) "If chosen, multipliers are removed from crosspoints / edge nodes" // LM_DUAL_CONSISTENT (default: boundary) (choices: all|boundary|none) "For which elements should the dual basis be calculated on EXACTLY the same GPs as the contact terms" // MESH_RELOCATION (default: Initial) (choices: Every_Timestep|Initial|None) "Type of mesh relocation" // ALGORITHM (default: Mortar) (choices: GPTS|LTL|LTS|Mortar|NTS|STL|gpts|ltl|lts|mortar|nts|stl) "Type of meshtying/contact algorithm" // INTTYPE (default: Segments) (choices: Elements|Elements_BS|Segments|elements|elements_BS|segments) "Type of numerical integration scheme" // NUMGP_PER_DIM (default: 0) "Number of employed integration points per dimension" // TRIANGULATION (default: Delaunay) (choices: Center|Delaunay|center|delaunay) "Type of triangulation for segment-based integration" // RESTART_WITH_MESHTYING (default: false) "Must be chosen if a non-meshtying simulation is to be restarted with meshtying" // OUTPUT_INTERFACES (default: false) "Write output for each mortar interface separately. This is an additional feature, purely to enhance visualization. Currently, this is limited to solid meshtying and contact w/o friction." .. _SECmortarcoupling_parallelredistribution: MORTAR COUPLING/PARALLEL REDISTRIBUTION --------------------------------------- MORTAR COUPLING/PARALLEL REDISTRIBUTION :: // MORTAR COUPLING/PARALLEL REDISTRIBUTION: // EXPLOIT_PROXIMITY (default: true) "Exploit information on geometric proximity to split slave interface into close and non-close parts and redistribute them independently. [Contact only]" // GHOSTING_STRATEGY (default: redundant_master) (choices: binning|redundant_all|redundant_master|round_robin) "Type of interface ghosting and ghosting extension algorithm" // IMBALANCE_TOL (default: 1.1) "Max. relative imbalance of subdomain size after redistribution" // MAX_BALANCE_EVAL_TIME (default: 2) "Max-to-min ratio of contact evaluation time per processor to trigger parallel redistribution" // MAX_BALANCE_SLAVE_ELES (default: 0.5) "Max-to-min ratio of mortar slave elements per processor to trigger parallel redistribution" // MIN_ELEPROC (default: 0) "Minimum no. of elements per processor for parallel redistribution" // PARALLEL_REDIST (default: Static) (choices: Dynamic|None|Static) "Type of redistribution algorithm" // PRINT_DISTRIBUTION (default: true) "Print details of the parallel distribution, i.e. number of nodes/elements for each rank." .. _SECmultipointconstraints: MULTI POINT CONSTRAINTS ----------------------- MULTI POINT CONSTRAINTS :: // MULTI POINT CONSTRAINTS: // RVE_REFERENCE_POINTS (default: automatic) "Method of definition of the reference points of an RVE" // ENFORCEMENT (default: penalty_method) (choices: lagrange_multiplier_method|penalty_method) "Method to enforce the multi point constraint" // PENALTY_PARAM (default: 100000) "Value of the penalty parameter" .. _SECnurbs: NURBS ----- NURBS :: // NURBS: // DO_LS_DBC_PROJECTION (default: false) "Determines if a projection is needed for least square Dirichlet boundary conditions." // SOLVER_LS_DBC_PROJECTION (default: -1) "Number of linear solver for the projection of least squares Dirichlet boundary conditions for NURBS discretizations" .. _SECparticledynamic: PARTICLE DYNAMIC ---------------- PARTICLE DYNAMIC :: // PARTICLE DYNAMIC: // DYNAMICTYPE (default: VelocityVerlet) (choices: SemiImplicitEuler|VelocityVerlet) "type of particle time integration" // INTERACTION (default: None) (choices: DEM|None|SPH) "type of particle interaction" // RESULTSEVERY (default: 1) "write particle runtime output every RESULTSEVERY steps" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // WRITE_GHOSTED_PARTICLES (default: false) "write ghosted particles (debug feature)" // TIMESTEP (default: 0.01) "time step size" // NUMSTEP (default: 100) "maximum number of steps" // MAXTIME (default: 1) "maximum time" // GRAVITY_ACCELERATION (default: 0.0 0.0 0.0) "acceleration due to gravity" // GRAVITY_RAMP_FUNCT (default: -1) "number of function governing gravity ramp" // VISCOUS_DAMPING (default: -1) "apply viscous damping force to determine static equilibrium solutions" // TRANSFER_EVERY (default: false) "transfer particles to new bins every time step" // PHASE_TO_DYNLOADBALFAC (default: none) "considered particle phases with dynamic load balance weighting factor" // PHASE_TO_MATERIAL_ID (default: none) "relate particle phase to material id" // INITIAL_POSITION_AMPLITUDE (default: 0.0 0.0 0.0) "amplitude of noise added to initial position for each spatial direction" // PARTICLE_WALL_SOURCE (default: NoParticleWall) (choices: BoundingBox|DiscretCondition|NoParticleWall) "type of particle wall source" // PARTICLE_WALL_MAT (default: -1) "material id for particle wall from bounding box source" // PARTICLE_WALL_MOVING (default: false) "consider a moving particle wall" // PARTICLE_WALL_LOADED (default: false) "consider loading on particle wall" // RIGID_BODY_MOTION (default: false) "consider rigid body motion" // RIGID_BODY_PHASECHANGE_RADIUS (default: -1) "search radius for neighboring rigid bodies in case of phase change" .. _SECparticledynamic_dem: PARTICLE DYNAMIC/DEM -------------------- PARTICLE DYNAMIC/DEM :: // PARTICLE DYNAMIC/DEM: // WRITE_PARTICLE_ENERGY (default: false) "write particle energy output" // WRITE_PARTICLE_WALL_INTERACTION (default: false) "write particle-wall interaction output" // NORMALCONTACTLAW (default: NormalLinearSpring) (choices: NormalHertz|NormalKuwabaraKono|NormalLeeHerrmann|NormalLinearSpring|NormalLinearSpringDamp|NormalTsuji) "normal contact law for particles" // TANGENTIALCONTACTLAW (default: NoTangentialContact) (choices: NoTangentialContact|TangentialLinSpringDamp) "tangential contact law for particles" // ROLLINGCONTACTLAW (default: NoRollingContact) (choices: NoRollingContact|RollingCoulomb|RollingViscous) "rolling contact law for particles" // ADHESIONLAW (default: NoAdhesion) (choices: AdhesionRegDMT|AdhesionVdWDMT|NoAdhesion) "type of adhesion law for particles" // ADHESION_SURFACE_ENERGY_DISTRIBUTION (default: ConstantSurfaceEnergy) (choices: ConstantSurfaceEnergy|LogNormalSurfaceEnergyDistribution|NormalSurfaceEnergyDistribution) "type of (random) surface energy distribution" // MIN_RADIUS (default: 0) "minimum allowed particle radius" // MAX_RADIUS (default: 0) "maximum allowed particle radius" // MAX_VELOCITY (default: -1) "maximum expected particle velocity" // INITIAL_RADIUS (default: RadiusFromParticleMaterial) (choices: LogNormalRadiusDistribution|NormalRadiusDistribution|RadiusFromParticleInput|RadiusFromParticleMaterial) "type of initial particle radius assignment" // RADIUSDISTRIBUTION_SIGMA (default: -1) "sigma of random particle radius distribution" // REL_PENETRATION (default: -1) "maximum allowed relative penetration" // NORMAL_STIFF (default: -1) "normal contact stiffness" // NORMAL_DAMP (default: -1) "normal contact damping parameter" // COEFF_RESTITUTION (default: -1) "coefficient of restitution" // DAMP_REG_FAC (default: -1) "linearly regularized damping normal force in the interval $|g| < (\\text{DAMP_REG_FAC} \\cdot r_{\\min})$" // TENSION_CUTOFF (default: true) "evaluate tension cutoff of normal contact force" // POISSON_RATIO (default: -1) "poisson ratio" // YOUNG_MODULUS (default: -1) "young's modulus" // FRICT_COEFF_TANG (default: -1) "friction coefficient for tangential contact" // FRICT_COEFF_ROLL (default: -1) "friction coefficient for rolling contact" // ADHESION_DISTANCE (default: -1) "adhesion distance between interacting surfaces" // ADHESION_MAX_CONTACT_PRESSURE (default: 0) "adhesion maximum contact pressure" // ADHESION_MAX_CONTACT_FORCE (default: 0) "adhesion maximum contact force" // ADHESION_USE_MAX_CONTACT_FORCE (default: false) "use maximum contact force instead of maximum contact pressure" // ADHESION_VDW_CURVE_SHIFT (default: false) "shifts van-der-Waals-curve to g = 0" // ADHESION_HAMAKER (default: -1) "hamaker constant of van-der-Waals interaction" // ADHESION_SURFACE_ENERGY (default: -1) "adhesion surface energy for the calculation of the pull-out force" // ADHESION_SURFACE_ENERGY_DISTRIBUTION_VAR (default: -1) "variance of adhesion surface energy distribution" // ADHESION_SURFACE_ENERGY_DISTRIBUTION_CUTOFF_FACTOR (default: -1) "adhesion surface energy distribution limited by multiple of variance" // ADHESION_SURFACE_ENERGY_FACTOR (default: 1) "factor to calculate minimum adhesion surface energy" .. _SECparticledynamic_initialandboundaryconditions: PARTICLE DYNAMIC/INITIAL AND BOUNDARY CONDITIONS ------------------------------------------------ PARTICLE DYNAMIC/INITIAL AND BOUNDARY CONDITIONS :: // PARTICLE DYNAMIC/INITIAL AND BOUNDARY CONDITIONS: // INITIAL_TEMP_FIELD (default: none) "Refer to the function ID describing the initial temperature field of particle phase" // INITIAL_VELOCITY_FIELD (default: none) "Refer to the function ID describing the initial velocity field of particle phase" // INITIAL_ANGULAR_VELOCITY_FIELD (default: none) "Refer to the function ID describing the initial angular velocity field of rigid body phase/DEM particle" // INITIAL_ACCELERATION_FIELD (default: none) "Refer to the function ID describing the initial acceleration field of particle phase" // INITIAL_ANGULAR_ACCELERATION_FIELD (default: none) "Refer to the function ID describing the initial angular acceleration field of rigid body phase/DEM particle" // DIRICHLET_BOUNDARY_CONDITION (default: none) "Refer to the function ID describing the dirichlet boundary condition of particle phase" // TEMPERATURE_BOUNDARY_CONDITION (default: none) "Refer to the function ID describing the temperature boundary condition of particle phase" .. _SECparticledynamic_sph: PARTICLE DYNAMIC/SPH -------------------- PARTICLE DYNAMIC/SPH :: // PARTICLE DYNAMIC/SPH: // WRITE_PARTICLE_WALL_INTERACTION (default: false) "write particle-wall interaction output" // KERNEL (default: CubicSpline) (choices: CubicSpline|QuinticSpline) "type of smoothed particle hydrodynamics kernel" // KERNEL_SPACE_DIM (default: Kernel3D) (choices: Kernel1D|Kernel2D|Kernel3D) "kernel space dimension number" // INITIALPARTICLESPACING (default: 0) "initial spacing of particles" // EQUATIONOFSTATE (default: GenTait) (choices: GenTait|IdealGas) "type of smoothed particle hydrodynamics equation of state" // MOMENTUMFORMULATION (default: AdamiMomentumFormulation) (choices: AdamiMomentumFormulation|MonaghanMomentumFormulation) "type of smoothed particle hydrodynamics momentum formulation" // DENSITYEVALUATION (default: DensitySummation) (choices: DensityIntegration|DensityPredictCorrect|DensitySummation) "type of density evaluation scheme" // DENSITYCORRECTION (default: NoCorrection) (choices: InteriorCorrection|NoCorrection|NormalizedCorrection|RandlesCorrection) "type of density correction scheme" // BOUNDARYPARTICLEFORMULATION (default: NoBoundaryFormulation) (choices: AdamiBoundaryFormulation|NoBoundaryFormulation) "type of boundary particle formulation" // BOUNDARYPARTICLEINTERACTION (default: NoSlipBoundaryParticle) (choices: FreeSlipBoundaryParticle|NoSlipBoundaryParticle) "type of boundary particle interaction" // WALLFORMULATION (default: NoWallFormulation) (choices: NoWallFormulation|VirtualParticleWallFormulation) "type of wall formulation" // TRANSPORTVELOCITYFORMULATION (default: NoTransportVelocity) (choices: GeneralizedTransportVelocity|NoTransportVelocity|StandardTransportVelocity) "type of transport velocity formulation" // TEMPERATUREEVALUATION (default: NoTemperatureEvaluation) (choices: NoTemperatureEvaluation|TemperatureIntegration) "type of temperature evaluation scheme" // TEMPERATUREGRADIENT (default: false) "evaluate temperature gradient" // HEATSOURCETYPE (default: NoHeatSource) (choices: NoHeatSource|SurfaceHeatSource|VolumeHeatSource) "type of heat source" // HEATSOURCE_FUNCT (default: -1) "number of function governing heat source" // HEATSOURCE_DIRECTION (default: 0.0 0.0 0.0) "direction of surface heat source" // VAPOR_HEATLOSS (default: false) "evaluate evaporation induced heat loss" // VAPOR_HEATLOSS_LATENTHEAT (default: 0) "latent heat in heat loss formula" // VAPOR_HEATLOSS_ENTHALPY_REFTEMP (default: 0) "enthalpy reference temperature in heat loss formula" // VAPOR_HEATLOSS_PFAC (default: 0) "pressure factor in heat loss formula" // VAPOR_HEATLOSS_TFAC (default: 0) "temperature factor in heat loss formula" // VAPOR_RECOIL (default: false) "evaluate evaporation induced recoil pressure" // VAPOR_RECOIL_BOILINGTEMPERATURE (default: 0) "boiling temperature in recoil pressure formula" // VAPOR_RECOIL_PFAC (default: 0) "pressure factor in recoil pressure formula" // VAPOR_RECOIL_TFAC (default: 0) "temperature factor in recoil pressure formula" // SURFACETENSIONFORMULATION (default: NoSurfaceTension) (choices: ContinuumSurfaceForce|NoSurfaceTension) "type of surface tension formulation" // SURFACETENSION_RAMP_FUNCT (default: -1) "number of function governing surface tension ramp" // SURFACETENSIONCOEFFICIENT (default: -1) "constant part of surface tension coefficient" // SURFACETENSIONMINIMUM (default: 0) "minimum surface tension coefficient in case of temperature dependence" // SURFACETENSIONTEMPFAC (default: 0) "factor of dependence of surface tension coefficient on temperature" // SURFACETENSIONREFTEMP (default: 0) "reference temperature for surface tension coefficient" // STATICCONTACTANGLE (default: 0) "static contact angle in degree with wetting effects" // TRIPLEPOINTNORMAL_CORR_CF_LOW (default: 0) "triple point normal correction wall color field low" // TRIPLEPOINTNORMAL_CORR_CF_UP (default: 0) "triple point normal correction wall color field up" // INTERFACE_VISCOSITY (default: false) "evaluate interface viscosity" // INTERFACE_VISCOSITY_LIQUIDGAS (default: 0) "artificial viscosity on liquid-gas interface" // INTERFACE_VISCOSITY_SOLIDLIQUID (default: 0) "artificial viscosity on solid-liquid interface" // BARRIER_FORCE (default: false) "evaluate barrier force" // BARRIER_FORCE_DISTANCE (default: 0) "barrier force distance" // BARRIER_FORCE_TEMPSCALE (default: 0) "barrier force temperature scaling" // BARRIER_FORCE_STIFF_HEAVY (default: -1) "barrier force stiffness of heavy phase" // BARRIER_FORCE_DAMP_HEAVY (default: 0) "barrier force damping parameter of heavy phase" // BARRIER_FORCE_STIFF_GAS (default: -1) "barrier force stiffness of gas phase" // BARRIER_FORCE_DAMP_GAS (default: 0) "barrier force damping parameter of gas phase" // TRANS_REF_TEMPERATURE (default: 0) "transition reference temperature" // TRANS_DT_SURFACETENSION (default: 0) "transition temperature difference for surface tension evaluation" // TRANS_DT_MARANGONI (default: 0) "transition temperature difference for marangoni evaluation" // TRANS_DT_CURVATURE (default: 0) "transition temperature difference for curvature evaluation" // TRANS_DT_WETTING (default: 0) "transition temperature difference for wetting evaluation" // TRANS_DT_INTVISC (default: 0) "transition temperature difference for interface viscosity evaluation" // TRANS_DT_BARRIER (default: 0) "transition temperature difference for barrier force evaluation" // DIRICHLETBOUNDARYTYPE (default: NoDirichletOpenBoundary) (choices: DirichletNormalToPlane|NoDirichletOpenBoundary) "type of dirichlet open boundary" // DIRICHLET_FUNCT (default: -1) "number of function governing velocity condition on dirichlet open boundary" // DIRICHLET_OUTWARD_NORMAL (default: 0.0 0.0 0.0) "direction of outward normal on dirichlet open boundary" // DIRICHLET_PLANE_POINT (default: 0.0 0.0 0.0) "point on dirichlet open boundary plane" // NEUMANNBOUNDARYTYPE (default: NoNeumannOpenBoundary) (choices: NeumannNormalToPlane|NoNeumannOpenBoundary) "type of neumann open boundary" // NEUMANN_FUNCT (default: -1) "number of function governing pressure condition on neumann open boundary" // NEUMANN_OUTWARD_NORMAL (default: 0.0 0.0 0.0) "direction of outward normal on neumann open boundary" // NEUMANN_PLANE_POINT (default: 0.0 0.0 0.0) "point on neumann open boundary plane" // PHASECHANGETYPE (default: NoPhaseChange) (choices: NoPhaseChange|OneWayScalarAboveToBelowPhaseChange|OneWayScalarBelowToAbovePhaseChange|TwoWayScalarPhaseChange) "type of phase change" // PHASECHANGEDEFINITION (default: none) "phase change definition" // RIGIDPARTICLECONTACTTYPE (default: NoRigidParticleContact) (choices: ElasticRigidParticleContact|NoRigidParticleContact) "type of rigid particle contact" // RIGIDPARTICLECONTACTSTIFF (default: -1) "rigid particle contact stiffness" // RIGIDPARTICLECONTACTDAMP (default: 0) "rigid particle contact damping parameter" .. _SECpasidynamic: PASI DYNAMIC ------------ PASI DYNAMIC :: // PASI DYNAMIC: // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // TIMESTEP (default: 0.01) "Time increment dt" // NUMSTEP (default: 100) "Total number of Timesteps" // MAXTIME (default: 1) "Total simulation time" // COUPLING (default: partitioned_onewaycoup) "partitioned coupling strategies for particle structure interaction" // ITEMAX (default: 10) "maximum number of partitioned iterations over fields" // CONVTOLSCALEDDISP (default: -1) "tolerance of dof and dt scaled interface displacement increments in partitioned iterations" // CONVTOLRELATIVEDISP (default: -1) "tolerance of relative interface displacement increments in partitioned iterations" // CONVTOLSCALEDFORCE (default: -1) "tolerance of dof and dt scaled interface force increments in partitioned iterations" // CONVTOLRELATIVEFORCE (default: -1) "tolerance of relative interface force increments in partitioned iterations" // IGNORE_CONV_CHECK (default: false) "ignore convergence check and proceed simulation" // STARTOMEGA (default: 1) "fixed relaxation parameter" // MAXOMEGA (default: 10) "largest omega allowed for Aitken relaxation" // MINOMEGA (default: 0.1) "smallest omega allowed for Aitken relaxation" .. _SECperiodicboundingboxgeometry: PERIODIC BOUNDINGBOX GEOMETRY ----------------------------- PERIODIC BOUNDINGBOX GEOMETRY :: // PERIODIC BOUNDINGBOX GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECporoelasticitydynamic: POROELASTICITY DYNAMIC ---------------------- POROELASTICITY DYNAMIC :: // POROELASTICITY DYNAMIC: // COUPALGO (default: poro_monolithic) (choices: poro_monolithic|poro_monolithicfluidsplit|poro_monolithicmeshtying|poro_monolithicnopenetrationsplit|poro_monolithicstructuresplit|poro_partitioned) "Coupling strategies for poroelasticity solvers" // PHYSICAL_TYPE (default: Poro) (choices: Poro|Poro_P1) "Physical Type of Porofluid" // TRANSIENT_TERMS (default: all) (choices: all|continuity|momentum|none) "which equation includes transient terms" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: 0.05) "time step size dt" // ITEMAX (default: 10) "maximum number of iterations over fields" // ITEMIN (default: 1) "minimal number of iterations over fields" // RESULTSEVERY (default: 1) "increment for writing solution" // TOLRES_GLOBAL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_GLOBAL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_DISP (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_DISP (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_PORO (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_PORO (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_VEL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_VEL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_PRES (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_PRES (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_NCOUP (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // NORM_INC (default: AbsSingleFields) (choices: AbsGlobal|AbsSingleFields) "type of norm for primary variables convergence check" // NORM_RESF (default: AbsSingleFields) (choices: AbsGlobal|AbsSingleFields) "type of norm for residual convergence check" // NORMCOMBI_RESFINC (default: And) (choices: And|Or) "binary operator to combine primary variables and residual force values" // VECTORNORM_RESF (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // VECTORNORM_INC (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // SECONDORDER (default: true) "Second order coupling at the interface." // CONTIPARTINT (default: false) "Partial integration of porosity gradient in continuity equation" // CONTACT_NO_PENETRATION (default: false) "No-Penetration Condition on active contact surface in case of poro contact problem!" // MATCHINGGRID (default: true) "is matching grid" // CONVECTIVE_TERM (default: false) "convective term" // LINEAR_SOLVER (default: -1) "number of linear solver used for poroelasticity problems" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" .. _SECporofluidmultiphasedynamic: POROFLUIDMULTIPHASE DYNAMIC --------------------------- POROFLUIDMULTIPHASE DYNAMIC :: // POROFLUIDMULTIPHASE DYNAMIC: // MAXTIME (default: 1000) "Total simulation time" // NUMSTEP (default: 20) "Total number of time steps" // TIMESTEP (default: 0.1) "Time increment dt" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // THETA (default: 0.5) "One-step-theta time integration factor" // TIMEINTEGR (default: One_Step_Theta) (choices: One_Step_Theta) "Time Integration Scheme" // CALCERROR (default: No) (choices: No|error_by_function) "compute error compared to analytical solution" // CALCERRORNO (default: -1) "function number for porofluidmultiphase error computation" // LINEAR_SOLVER (default: -1) "number of linear solver used for the porofluidmultiphase problem" // ITEMAX (default: 10) "max. number of nonlin. iterations" // ABSTOLRES (default: 1e-14) "Absolute tolerance for deciding if residual of nonlinear problem is already zero" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // FDCHECK (default: none) (choices: global|none) "flag for finite difference check: none, local, or global" // FDCHECKEPS (default: 1e-06) "dof perturbation magnitude for finite difference check (1.e-6 seems to work very well, whereas smaller values don't)" // FDCHECKTOL (default: 1e-06) "relative tolerance for finite difference check" // SKIPINITDER (default: true) "Flag to skip computation of initial time derivative" // OUTPUT_SATANDPRESS (default: true) "Flag if output of saturations and pressures should be calculated" // OUTPUT_SOLIDPRESS (default: true) "Flag if output of solid pressure should be calculated" // OUTPUT_POROSITY (default: true) "Flag if output of porosity should be calculated" // OUTPUT_PHASE_VELOCITIES (default: true) "Flag if output of phase velocities should be calculated" // STAB_BIOT (default: false) "Flag to (de)activate BIOT stabilization." // STAB_BIOT_SCALING (default: 1) "Scaling factor for stabilization parameter for biot stabilization of porous flow." // VECTORNORM_RESF (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // VECTORNORM_INC (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // TOLRES (default: 1e-06) "tolerance in the residual norm for the Newton iteration" // TOLINC (default: 1e-06) "tolerance in the increment norm for the Newton iteration" // INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "Initial Field for the porofluid problem" // INITFUNCNO (default: -1) "function number for scalar transport initial field" // DIVERCONT (default: stop) (choices: continue|stop) "What to do with time integration when Newton-Raphson iteration failed" // FLUX_PROJ_SOLVER (default: -1) "Number of linear solver used for L2 projection" // FLUX_PROJ_METHOD (default: none) (choices: L2_projection|none) "Flag to (de)activate flux reconstruction." // DOMAININT_FUNCT (default: -1.0) "functions used for domain integrals" // ARTERY_COUPLING (default: false) "Coupling with 1D blood vessels." // STARTING_DBC_TIME_END (default: -1) "End time for the starting Dirichlet BC." // STARTING_DBC_ONOFF (default: 0) "Switching the starting Dirichlet BC on or off." // STARTING_DBC_FUNCT (default: 0) "Function prescribing the starting Dirichlet BC." .. _SECporofluidmultiphasedynamic_arterycoupling: POROFLUIDMULTIPHASE DYNAMIC/ARTERY COUPLING ------------------------------------------- POROFLUIDMULTIPHASE DYNAMIC/ARTERY COUPLING :: // POROFLUIDMULTIPHASE DYNAMIC/ARTERY COUPLING: // MAXNUMSEGPERARTELE (default: 5) "maximum number of segments per artery element for 1D-3D artery coupling" // PENALTY (default: 1000) "Penalty parameter for line-based coupling" // ARTERY_COUPLING_METHOD (default: None) (choices: GPTS|MP|NTP|Nodal|None) "Coupling method for artery coupling." // COUPLEDDOFS_ART (default: -1.0) "coupled artery dofs for mesh tying" // COUPLEDDOFS_PORO (default: -1.0) "coupled porofluid dofs for mesh tying" // REACFUNCT_ART (default: -1) "functions for coupling (artery part)" // SCALEREAC_ART (default: 0) "scale for coupling (artery part)" // REACFUNCT_CONT (default: -1) "functions for coupling (porofluid part)" // SCALEREAC_CONT (default: 0) "scale for coupling (porofluid part)" // EVALUATE_IN_REF_CONFIG (default: true) "Flag if artery elements are evaluated in reference or current configuration" // LATERAL_SURFACE_COUPLING (default: false) "Flag if 1D-3D coupling should be evaluated on lateral (cylinder) surface of embedded artery elements" // NUMPATCH_AXI (default: 1) "Number of integration patches per 1D element in axial direction for lateral surface coupling" // NUMPATCH_RAD (default: 1) "Number of integration patches per 1D element in radial direction for lateral surface coupling" // OUTPUT_BLOODVESSELVOLFRAC (default: false) "Flag if output of blood vessel volume fraction should be calculated" // PRINT_OUT_SUMMARY_PAIRS (default: false) "Flag if summary of coupling-pairs should be printed" // DELETE_FREE_HANGING_ELES (default: false) "Flag if free-hanging elements (after blood vessel collapse) should be deleted" // DELETE_SMALL_FREE_HANGING_COMPS (default: -1) "Small connected components whose size is smaller than this fraction of the overall network size are additionally deleted (a valid choice of this parameter should lie between 0 and 1)" .. _SECporomultiphasedynamic: POROMULTIPHASE DYNAMIC ---------------------- POROMULTIPHASE DYNAMIC :: // POROMULTIPHASE DYNAMIC: // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: -1) "time step size dt" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // SOLVE_STRUCTURE (default: true) "Flag to skip computation of structural field" // COUPALGO (default: twoway_partitioned) (choices: twoway_monolithic|twoway_partitioned) "Coupling strategies for porofluid-elasticity solvers" // ARTERY_COUPLING (default: false) "Coupling with 1D blood vessels." .. _SECporomultiphasedynamic_monolithic: POROMULTIPHASE DYNAMIC/MONOLITHIC --------------------------------- POROMULTIPHASE DYNAMIC/MONOLITHIC :: // POROMULTIPHASE DYNAMIC/MONOLITHIC: // TOLRES_GLOBAL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_GLOBAL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // LINEAR_SOLVER (default: -1) "number of linear solver used for poroelasticity problems" // FDCHECK (default: none) (choices: global|none) "flag for finite difference check: none or global" // VECTORNORM_RESF (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // VECTORNORM_INC (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.001) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" .. _SECporomultiphasedynamic_partitioned: POROMULTIPHASE DYNAMIC/PARTITIONED ---------------------------------- POROMULTIPHASE DYNAMIC/PARTITIONED :: // POROMULTIPHASE DYNAMIC/PARTITIONED: // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteration" // RELAXATION (default: none) (choices: Aitken|Constant|none) "flag for relaxation of partitioned scheme" // STARTOMEGA (default: 1) "fixed relaxation parameter" // MINOMEGA (default: 0.1) "smallest omega allowed for Aitken relaxation" // MAXOMEGA (default: 10) "largest omega allowed for Aitken relaxation" .. _SECporomultiphasescatradynamic: POROMULTIPHASESCATRA DYNAMIC ---------------------------- POROMULTIPHASESCATRA DYNAMIC :: // POROMULTIPHASESCATRA DYNAMIC: // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: 0.05) "time step size dt" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // ITEMIN (default: 1) "minimal number of iterations over fields" // COUPALGO (default: twoway_partitioned_nested) (choices: twoway_monolithic|twoway_partitioned_nested|twoway_partitioned_sequential) "Coupling strategies for poroscatra solvers" // ARTERY_COUPLING (default: false) "Coupling with 1D blood vessels." // DIVERCONT (default: stop) (choices: continue|stop) "What to do with time integration when Poromultiphase-Scatra iteration failed" .. _SECporomultiphasescatradynamic_monolithic: POROMULTIPHASESCATRA DYNAMIC/MONOLITHIC --------------------------------------- POROMULTIPHASESCATRA DYNAMIC/MONOLITHIC :: // POROMULTIPHASESCATRA DYNAMIC/MONOLITHIC: // VECTORNORM_RESF (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // VECTORNORM_INC (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.001) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // TOLRES_GLOBAL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_GLOBAL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // LINEAR_SOLVER (default: -1) "number of linear solver used for monolithic poroscatra problems" // FDCHECK (default: none) (choices: global|none) "flag for finite difference check: none or global" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" .. _SECporomultiphasescatradynamic_partitioned: POROMULTIPHASESCATRA DYNAMIC/PARTITIONED ---------------------------------------- POROMULTIPHASESCATRA DYNAMIC/PARTITIONED :: // POROMULTIPHASESCATRA DYNAMIC/PARTITIONED: // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteration" .. _SECporoscatracontrol: POROSCATRA CONTROL ------------------ POROSCATRA CONTROL :: // POROSCATRA CONTROL: // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: 0.05) "time step size dt" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // ITEMIN (default: 1) "minimal number of iterations over fields" // TOLRES_GLOBAL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_GLOBAL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_DISP (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_DISP (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_VEL (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_VEL (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_PRES (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_PRES (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // TOLRES_SCALAR (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // TOLINC_SCALAR (default: 1e-08) "tolerance in the increment norm for the Newton iteration" // NORM_INC (default: AbsSingleFields) (choices: AbsGlobal|AbsSingleFields) "type of norm for primary variables convergence check" // NORM_RESF (default: AbsSingleFields) (choices: AbsGlobal|AbsSingleFields) "type of norm for residual convergence check" // NORMCOMBI_RESFINC (default: And) (choices: And|Or) "binary operator to combine primary variables and residual force values" // VECTORNORM_RESF (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // VECTORNORM_INC (default: L2) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // LINEAR_SOLVER (default: -1) "number of linear solver used for monolithic poroscatra problems" // COUPALGO (default: solid_to_scatra) (choices: monolithic|scatra_to_solid|solid_to_scatra|two_way) "Coupling strategies for poroscatra solvers" // MATCHINGGRID (default: true) "is matching grid" .. _SECproblemsize: PROBLEM SIZE ------------ PROBLEM SIZE :: // PROBLEM SIZE: // DIM (default: 3) "2d or 3d problem" // ELEMENTS (default: 0) "Total number of elements" // NODES (default: 0) "Total number of nodes" // NPATCHES (default: 0) "number of nurbs patches" // MATERIALS (default: 0) "number of materials" // NUMDF (default: 3) "maximum number of degrees of freedom" .. _SECproblemtype: PROBLEM TYPE ------------ PROBLEM TYPE :: // PROBLEM TYPE: // PROBLEMTYPE (choices: Ale|ArterialNetwork|Biofilm_Fluid_Structure_Interaction|Cardiac_Monodomain|Elastohydrodynamic_Lubrication|Electrochemistry|Fluid|Fluid_Ale|Fluid_Beam_Interaction|Fluid_Poro_Structure_Interaction_XFEM|Fluid_Porous_Structure_Interaction|Fluid_Porous_Structure_Scalar_Scalar_Interaction|Fluid_RedModels|Fluid_Structure_Interaction|Fluid_Structure_Interaction_RedModels|Fluid_Structure_Interaction_XFEM|Fluid_XFEM|Fluid_XFEM_LevelSet|Gas_Fluid_Structure_Interaction|Level_Set|Low_Mach_Number_Flow|Lubrication|Multiphase_Poroelasticity|Multiphase_Poroelasticity_ScaTra|Multiphase_Porous_Flow|NP_Supporting_Procs|Particle|Particle_Structure_Interaction|Polymer_Network|Poroelastic_scalar_transport|Poroelasticity|RedAirways_Tissue|ReducedDimensionalAirWays|Reduced_Lung|Scalar_Thermo_Interaction|Scalar_Transport|Structure|Structure_Scalar_Interaction|Structure_Scalar_Thermo_Interaction|Thermo|Thermo_Fluid_Structure_Interaction|Thermo_Structure_Interaction) "Type of the problem" // SHAPEFCT (default: Polynomial) (choices: HDG|Nurbs|Polynomial) "Defines the function spaces for the spatial approximation" // RESTART (default: 0) // RANDSEED (default: -1) "Set the random seed. If < 0 use current time." .. _SECreduceddairwaysgeometry: REDUCED D AIRWAYS GEOMETRY -------------------------- REDUCED D AIRWAYS GEOMETRY :: // REDUCED D AIRWAYS GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECreduceddimensionalairwaysdynamic: REDUCED DIMENSIONAL AIRWAYS DYNAMIC ----------------------------------- REDUCED DIMENSIONAL AIRWAYS DYNAMIC :: // REDUCED DIMENSIONAL AIRWAYS DYNAMIC: // DYNAMICTYPE (default: OneStepTheta) (choices: OneStepTheta) "OneStepTheta Scheme" // SOLVERTYPE (default: Linear) (choices: Linear|Nonlinear) "Solver type" // TIMESTEP (default: 0.01) "Time increment dt" // NUMSTEP (default: 0) "Number of Time Steps" // RESTARTEVERY (default: 1) "Increment for writing restart" // RESULTSEVERY (default: 1) "Increment for writing solution" // THETA (default: 1) "One-step-theta time integration factor" // MAXITERATIONS (default: 1) "maximum iteration steps" // TOLERANCE (default: 1e-06) "tolerance" // LINEAR_SOLVER (default: -1) "number of linear solver used for reduced dim arterial dynamics" // SOLVESCATRA (default: false) "Flag to (de)activate solving scalar transport in blood" // COMPAWACINTER (default: false) "Flag to (de)activate computation of airway-acinus interdependency" // CALCV0PRESTRESS (default: false) "Flag to (de)activate initial acini volume adjustment with pre-stress condition" // TRANSPULMPRESS (default: 800) "Transpulmonary pressure needed for recalculation of acini volumes" .. _SECscalartransportdynamic: SCALAR TRANSPORT DYNAMIC ------------------------ SCALAR TRANSPORT DYNAMIC :: // SCALAR TRANSPORT DYNAMIC: // SOLVERTYPE (default: linear_full) (choices: linear_full|linear_incremental|nonlinear|nonlinear_multiscale_macrotomicro|nonlinear_multiscale_macrotomicro_aitken|nonlinear_multiscale_macrotomicro_aitken_dofsplit|nonlinear_multiscale_microtomacro) "type of scalar transport solver" // TIMEINTEGR (default: One_Step_Theta) (choices: BDF2|Gen_Alpha|One_Step_Theta|Stationary) "Time Integration Scheme" // MAXTIME (default: 1000) "Total simulation time" // NUMSTEP (default: 20) "Total number of time steps" // TIMESTEP (default: 0.1) "Time increment dt" // THETA (default: 0.5) "One-step-theta time integration factor" // ALPHA_M (default: 0.5) "Generalized-alpha time integration factor" // ALPHA_F (default: 0.5) "Generalized-alpha time integration factor" // GAMMA (default: 0.5) "Generalized-alpha time integration factor" // RESULTSEVERY (default: 1) "Increment for writing solution" // RESTARTEVERY (default: 1) "Increment for writing restart" // MATID (default: -1) "Material ID for automatic mesh generation" // VELOCITYFIELD (default: zero) (choices: Navier_Stokes|function|zero) "type of velocity field used for scalar transport problems" // VELFUNCNO (default: -1) "function number for scalar transport velocity field" // INITIALFIELD (default: zero_field) (choices: 1D_DISCONTPV|FLAME_VORTEX_INTERACTION|L_shaped_domain|RAYTAYMIXFRAC|algebraic_field_dependence|disturbed_field_by_function|facing_flame_fronts|field_by_condition|field_by_function|high_forced_hit|low_forced_hit|oracles_flame|zero_field) "Initial Field for transport problem" // INITFUNCNO (default: -1) "function number for scalar transport initial field" // SPHERICALCOORDS (default: false) "use of spherical coordinates" // CALCERROR (default: No) (choices: AnalyticSeries|ConcentricCylinders|Electroneutrality|Kwok_Wu|No|SphereDiffusion|error_by_condition|error_by_function) "compute error compared to analytical solution" // CALCERRORNO (default: -1) "function number for scalar transport error computation" // CALCFLUX_DOMAIN (default: No) (choices: No|diffusive|total) "output of diffusive/total flux vectors inside domain" // CALCFLUX_DOMAIN_LUMPED (default: true) "perform approximate domain flux calculation involving matrix lumping" // CALCFLUX_BOUNDARY (default: No) (choices: No|convective|diffusive|total) "output of convective/diffusive/total flux vectors on boundary" // CALCFLUX_BOUNDARY_LUMPED (default: true) "perform approximate boundary flux calculation involving matrix lumping" // WRITEFLUX_IDS (default: -1) "Write diffusive/total flux vector fields for these scalar fields only (starting with 1)" // OUTPUTSCALARS (default: none) (choices: by_condition|entire_domain|entire_domain_and_by_condition|none) "Output of total and mean values for transported scalars" // OUTPUTSCALARSMEANGRAD (default: false) "Output of mean gradient of scalars" // OUTINTEGRREAC (default: false) "Output of integral reaction values" // OUTPUT_GMSH (default: false) "Do you want to write Gmsh postprocessing files?" // MATLAB_STATE_OUTPUT (default: false) "Do you want to write the state solution to Matlab file?" // CONVFORM (default: convective) (choices: conservative|convective) "form of convective term" // NEUMANNINFLOW (default: false) "Flag to (de)activate potential Neumann inflow term(s)" // CONV_HEAT_TRANS (default: false) "Flag to (de)activate potential convective heat transfer boundary conditions" // SKIPINITDER (default: false) "Flag to skip computation of initial time derivative" // FSSUGRDIFF (default: No) (choices: No|Smagorinsky_all|Smagorinsky_small|artificial) "fine-scale subgrid diffusivity" // MESHTYING (default: no) (choices: Condensed_Bmat|Condensed_Bmat_merged|Condensed_Smat|no) "Flag to (de)activate mesh tying algorithm" // FIELDCOUPLING (default: matching) (choices: matching|volmortar) "Type of coupling strategy between fields" // LINEAR_SOLVER (default: -1) "number of linear solver used for scalar transport/elch..." // L2_PROJ_LINEAR_SOLVER (default: -1) "number of linear solver used for l2-projection sub-problems" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" // MATRIXTYPE (default: sparse) "type of global system matrix in global system of equations" // NATURAL_CONVECTION (default: false) "Include natural convection effects" // FDCHECK (default: none) (choices: global|global_extended|local|none) "flag for finite difference check: none, local, or global" // FDCHECKEPS (default: 1e-06) "dof perturbation magnitude for finite difference check (1.e-6 seems to work very well, whereas smaller values don't)" // FDCHECKTOL (default: 1e-06) "relative tolerance for finite difference check" // COMPUTEINTEGRALS (default: none) (choices: initial|none|repeated) "flag for optional computation of domain integrals" // PADAPTIVITY (default: false) "Flag to (de)activate p-adativity" // PADAPTERRORTOL (default: 1e-06) "The error tolerance to calculate the variation of the elemental degree" // PADAPTERRORBASE (default: 1.66) "The error tolerance base to calculate the variation of the elemental degree" // PADAPTDEGREEMAX (default: 4) "The max. degree of the shape functions" // SEMIIMPLICIT (default: false) "Flag to (de)activate semi-implicit calculation of the reaction term" // OUTPUTLINSOLVERSTATS (default: false) "flag for output of performance statistics associated with linear solver into csv file" // OUTPUTNONLINSOLVERSTATS (default: false) "flag for output of performance statistics associated with nonlinear solver into csv file" // NULLSPACE_POINTBASED (default: false) "flag for point-based null space calculation" // ADAPTIVE_TIMESTEPPING (default: false) "flag for adaptive time stepping" .. _SECscalartransportdynamic_arterycoupling: SCALAR TRANSPORT DYNAMIC/ARTERY COUPLING ---------------------------------------- SCALAR TRANSPORT DYNAMIC/ARTERY COUPLING :: // SCALAR TRANSPORT DYNAMIC/ARTERY COUPLING: // ARTERY_COUPLING_METHOD (default: None) (choices: GPTS|MP|NTP|Nodal|None) "Coupling method for artery coupling." // PENALTY (default: 1000) "Penalty parameter for line-based coupling" // COUPLEDDOFS_ARTSCATRA (default: -1.0) "coupled artery dofs for mesh tying" // COUPLEDDOFS_SCATRA (default: -1.0) "coupled porofluid dofs for mesh tying" // REACFUNCT_ART (default: -1) "functions for coupling (arteryscatra part)" // SCALEREAC_ART (default: 0) "scale for coupling (arteryscatra part)" // REACFUNCT_CONT (default: -1) "functions for coupling (scatra part)" // SCALEREAC_CONT (default: 0) "scale for coupling (scatra part)" .. _SECscalartransportdynamic_externalforce: SCALAR TRANSPORT DYNAMIC/EXTERNAL FORCE --------------------------------------- SCALAR TRANSPORT DYNAMIC/EXTERNAL FORCE :: // SCALAR TRANSPORT DYNAMIC/EXTERNAL FORCE: // EXTERNAL_FORCE (default: false) "Flag to activate external force" // FORCE_FUNCTION_ID (default: -1) "Function ID for external force" // INTRINSIC_MOBILITY_FUNCTION_ID (default: -1) "Function ID for intrinsic mobility" .. _SECscalartransportdynamic_nonlinear: SCALAR TRANSPORT DYNAMIC/NONLINEAR ---------------------------------- SCALAR TRANSPORT DYNAMIC/NONLINEAR :: // SCALAR TRANSPORT DYNAMIC/NONLINEAR: // ITEMAX (default: 10) "max. number of nonlin. iterations" // CONVTOL (default: 1e-06) "Tolerance for convergence check" // ITEMAX_OUTER (default: 10) "Maximum number of outer iterations in partitioned coupling schemes (natural convection, multi-scale simulations etc.)" // CONVTOL_OUTER (default: 1e-06) "Convergence check tolerance for outer loop in partitioned coupling schemes (natural convection, multi-scale simulations etc.)" // EXPLPREDICT (default: false) "do an explicit predictor step before starting nonlinear iteration" // ABSTOLRES (default: 1e-14) "Absolute tolerance for deciding if residual of nonlinear problem is already zero" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" .. _SECscalartransportdynamic_s2icoupling: SCALAR TRANSPORT DYNAMIC/S2I COUPLING ------------------------------------- SCALAR TRANSPORT DYNAMIC/S2I COUPLING :: // SCALAR TRANSPORT DYNAMIC/S2I COUPLING: // COUPLINGTYPE (default: Undefined) (choices: CondensedMortar_Bubnov|CondensedMortar_Petrov|MatchingNodes|SaddlePointMortar_Bubnov|SaddlePointMortar_Petrov|StandardMortar|StandardNodeToSegment|Undefined) "type of mortar meshtying" // LMSIDE (default: slave) (choices: master|slave) "flag for interface side underlying Lagrange multiplier definition" // SLAVEONLY (default: false) "flag for evaluation of interface linearizations and residuals on slave side only" // NTSPROJTOL (default: 0) "node-to-segment projection tolerance" // INTLAYERGROWTH_EVALUATION (default: none) (choices: monolithic|none|semi-implicit) "flag for evaluation of scatra-scatra interface coupling involving interface layer growth" // INTLAYERGROWTH_CONVTOL (default: 1e-12) "local Newton-Raphson convergence tolerance for scatra-scatra interface coupling involving interface layer growth" // INTLAYERGROWTH_ITEMAX (default: 5) "maximum number of local Newton-Raphson iterations for scatra-scatra interface coupling involving interface layer growth" // INTLAYERGROWTH_LINEAR_SOLVER (default: -1) "ID of linear solver for monolithic scatra-scatra interface coupling involving interface layer growth" // INTLAYERGROWTH_TIMESTEP (default: -1) "modified time step size for scatra-scatra interface coupling involving interface layer growth" // MESHTYING_CONDITIONS_INDEPENDENT_SETUP (default: false) "mesh tying for different conditions should be setup independently" // OUTPUT_INTERFACE_FLUX (default: false) "evaluate integral of coupling flux on slave side for each s2i condition and write it to csv file" .. _SECscalartransportdynamic_stabilization: SCALAR TRANSPORT DYNAMIC/STABILIZATION -------------------------------------- SCALAR TRANSPORT DYNAMIC/STABILIZATION :: // SCALAR TRANSPORT DYNAMIC/STABILIZATION: // STABTYPE (default: SUPG) (choices: GLS|SUPG|USFEM|centered|no_stabilization|upwind) "Type of stabilization (if any). No stabilization is only reasonable for low-Peclet-number." // SUGRVEL (default: false) "potential incorporation of subgrid-scale velocity" // ASSUGRDIFF (default: false) "potential incorporation of all-scale subgrid diffusivity (a.k.a. discontinuity-capturing) term" // DEFINITION_TAU (default: Franca_Valentin) (choices: Codina|Codina_wo_dt|Exact_1D|Franca_Madureira_Valentin|Franca_Madureira_Valentin_wo_dt|Franca_Valentin|Franca_Valentin_wo_dt|Numerical_Value|Shakib_Hughes_Codina|Shakib_Hughes_Codina_wo_dt|Taylor_Hughes_Zarins|Taylor_Hughes_Zarins_wo_dt|Zero) "Definition of tau" // CHARELELENGTH (default: streamlength) "Characteristic element length for tau" // DEFINITION_ASSGD (default: artificial_linear) (choices: Almeida_Silva_97_nonlinear|Codina_nonlinear|Hughes_etal_86_nonlinear|Tezduyar_Park_86_nonlinear|Tezduyar_Park_86_nonlinear_wo_phizero|YZbeta_nonlinear|artificial_linear|artificial_linear_reinit|doCarmo_Galeao_91_nonlinear) "Definition of (all-scale) subgrid diffusivity" // EVALUATION_TAU (default: element_center) (choices: element_center|integration_point) "Location where tau is evaluated" // EVALUATION_MAT (default: element_center) (choices: element_center|integration_point) "Location where material law is evaluated" // CONSISTENCY (default: no) (choices: L2_projection_lumped|no) "improvement of consistency for stabilization" // TAU_VALUE (default: 0) "Numerical value for tau for stabilization" .. _SECsearchtree: SEARCH TREE ----------- SEARCH TREE :: // SEARCH TREE: // TREE_TYPE (default: notree) (choices: notree|octree3d|quadtree2d|quadtree3d) "set tree type" .. _SECsemi-smoothplasticity: SEMI-SMOOTH PLASTICITY ---------------------- SEMI-SMOOTH PLASTICITY :: // SEMI-SMOOTH PLASTICITY: // SEMI_SMOOTH_CPL (default: 1) "Weighting factor cpl for semi-smooth PDASS" // STABILIZATION_S (default: 1) "Stabilization factor s for semi-smooth PDASS" // NORMCOMBI_RESFPLASTCONSTR (default: And) (choices: And|Or) "binary operator to combine plasticity constraints and residual force values" // NORMCOMBI_DISPPLASTINCR (default: And) (choices: And|Or) "binary operator to combine displacement increments and plastic flow (Delta Lp) increment values" // TOLPLASTCONSTR (default: 1e-08) "tolerance in the plastic constraint norm for the newton iteration" // TOLDELTALP (default: 1e-08) "tolerance in the plastic flow (Delta Lp) norm for the Newton iteration" // NORMCOMBI_EASRES (default: And) (choices: And|Or) "binary operator to combine EAS-residual and residual force values" // NORMCOMBI_EASINCR (default: And) (choices: And|Or) "binary operator to combine displacement increments and EAS increment values" // TOLEASRES (default: 1e-08) "tolerance in the EAS residual norm for the newton iteration" // TOLEASINCR (default: 1e-08) "tolerance in the EAS increment norm for the Newton iteration" // DISSIPATION_MODE (default: pl_multiplier) (choices: Taylor_Quinney|pl_flow|pl_multiplier) "method to calculate the plastic dissipation" .. _SECsolver1: SOLVER 1 -------- SOLVER 1 :: // SOLVER 1: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver2: SOLVER 2 -------- SOLVER 2 :: // SOLVER 2: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver3: SOLVER 3 -------- SOLVER 3 :: // SOLVER 3: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver4: SOLVER 4 -------- SOLVER 4 :: // SOLVER 4: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver5: SOLVER 5 -------- SOLVER 5 :: // SOLVER 5: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver6: SOLVER 6 -------- SOLVER 6 :: // SOLVER 6: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver7: SOLVER 7 -------- SOLVER 7 :: // SOLVER 7: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver8: SOLVER 8 -------- SOLVER 8 :: // SOLVER 8: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECsolver9: SOLVER 9 -------- SOLVER 9 :: // SOLVER 9: // SOLVER (default: undefined) (choices: Belos|Superlu|UMFPACK|undefined) "The solver to attack the system of linear equations arising of FE approach with." // AZSOLVE (default: GMRES) (choices: BiCGSTAB|CG|GMRES) "Type of linear solver algorithm to use." // AZPREC (default: ILU) (choices: AMGnxn|ILU|MueLu|Teko) "Type of internal preconditioner to use. Note! this preconditioner will only be used if the input operator supports the Epetra_RowMatrix interface and the client does not pass in an external preconditioner!" // IFPACKOVERLAP (default: 0) "The amount of overlap used for the ifpack \"ilu\" preconditioner." // IFPACKGFILL (default: 0) "The amount of fill allowed for an internal \"ilu\" preconditioner." // IFPACKCOMBINE (default: Add) (choices: Add|Insert|Zero) "Combine mode for Ifpack Additive Schwarz" // AZITER (default: 1000) "The maximum number of iterations the underlying iterative solver is allowed to perform" // AZTOL (default: 1e-08) "The level the residual norms must reach to decide about successful convergence" // AZCONV (default: AZ_r0) (choices: AZ_noscaled|AZ_r0) "The implicit residual norm scaling type to use for terminating the iterative solver." // AZOUTPUT (default: 0) "The number of iterations between each output of the solver's progress is written to screen" // AZREUSE (default: 0) "The number specifying how often to recompute some preconditioners" // AZSUB (default: 50) "The maximum size of the Krylov subspace used with \"GMRES\" before a restart is performed." // THROW_IF_UNCONVERGED (default: true) "If set to true, the iterative linear solver will throw an exception if it does not converge. To only issue a warning without stopping the simulation, set this parameter to false." // SOLVER_XML_FILE > (default: none) "xml file defining any iterative solver" // MUELU_XML_FILE > (default: none) "xml file defining any MueLu preconditioner" // TEKO_XML_FILE > (default: none) "xml file defining any Teko preconditioner" // NAME (default: No_name) "User specified name for solver block" // AMGNXN_TYPE (default: AMG(BGS)) "Name of the pre-built preconditioner to be used. If set to\"XML\" the preconditioner is defined using a xml file" // AMGNXN_XML_FILE > (default: none) "xml file defining the AMGnxn preconditioner" .. _SECssicontrol: SSI CONTROL ----------- SSI CONTROL :: // SSI CONTROL: // RESTARTEVERYTIME (default: 0) "write restart possibility every RESTARTEVERY steps" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: -1) "time step size dt" // DIFFTIMESTEPSIZE (default: false) "use different step size for scatra and solid" // RESULTSEVERYTIME (default: 0) "increment for writing solution" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // SCATRA_FROM_RESTART_FILE (default: false) "read scatra result from restart files (use option 'restartfromfile' during execution of 4C)" // SCATRA_FILENAME (default: nil) "Control-file name for reading scatra results in SSI" // FIELDCOUPLING (default: volume_matching) (choices: boundary_nonmatching|volume_matching|volume_nonmatching|volumeboundary_matching) "Type of coupling strategy between fields" // COUPALGO (default: ssi_IterStagg) "Coupling strategies for SSI solvers" // SCATRATIMINTTYPE (default: Standard) (choices: Cardiac_Monodomain|Elch|Standard) "scalar transport time integration type is needed to instantiate correct scalar transport time integration scheme for ssi problems" // RESTART_FROM_STRUCTURE (default: false) "restart from structure problem (e.g. from prestress calculations) instead of ssi" // ADAPTIVE_TIMESTEPPING (default: false) "flag for adaptive time stepping" // REDISTRIBUTE_SOLID (default: false) "redistribution by binning of solid mechanics discretization" .. _SECssicontrol_elch: SSI CONTROL/ELCH ---------------- SSI CONTROL/ELCH :: // SSI CONTROL/ELCH: // INITPOTCALC (default: false) "Automatically calculate initial field for electric potential" .. _SECssicontrol_manifold: SSI CONTROL/MANIFOLD -------------------- SSI CONTROL/MANIFOLD :: // SSI CONTROL/MANIFOLD: // ADD_MANIFOLD (default: false) "activate additional manifold?" // MESHTYING_MANIFOLD (default: false) "activate meshtying between all manifold fields in case they intersect?" // INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "Initial field for scalar transport on manifold" // INITFUNCNO (default: -1) "function number for scalar transport on manifold initial field" // LINEAR_SOLVER (default: -1) "linear solver for scalar transport on manifold" // OUTPUT_INFLOW (default: false) "write output of inflow of scatra manifold - scatra coupling into scatra manifold to csv file" .. _SECssicontrol_monolithic: SSI CONTROL/MONOLITHIC ---------------------- SSI CONTROL/MONOLITHIC :: // SSI CONTROL/MONOLITHIC: // ABSTOLRES (default: 1e-14) "absolute tolerance for deciding if global residual of nonlinear problem is already zero" // CONVTOL (default: 1e-06) "tolerance for convergence check of Newton-Raphson iteration within monolithic SSI" // LINEAR_SOLVER (default: -1) "ID of linear solver for global system of equations" // MATRIXTYPE (default: undefined) (choices: block|sparse|undefined) "type of global system matrix in global system of equations" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" // EQUILIBRATION_STRUCTURE (default: none) "flag for equilibration of structural equations" // EQUILIBRATION_SCATRA (default: none) "flag for equilibration of scatra equations" // PRINT_MAT_RHS_MAP_MATLAB (default: false) "print system matrix, rhs vector, and full map to matlab readable file after solution of time step" // RELAX_LIN_SOLVER_TOLERANCE (default: 1) "relax the tolerance of the linear solver in case it is an iterative solver by scaling the convergence tolerance with factor RELAX_LIN_SOLVER_TOLERANCE" // RELAX_LIN_SOLVER_STEP (default: -1) "relax the tolerance of the linear solver within the first RELAX_LIN_SOLVER_STEP steps" .. _SECssicontrol_partitioned: SSI CONTROL/PARTITIONED ----------------------- SSI CONTROL/PARTITIONED :: // SSI CONTROL/PARTITIONED: // MAXOMEGA (default: 10) "largest omega allowed for Aitken relaxation" // MINOMEGA (default: 0.1) "smallest omega allowed for Aitken relaxation" // STARTOMEGA (default: 1) "fixed relaxation parameter" // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteration within partitioned SSI" .. _SECssticontrol: SSTI CONTROL ------------ SSTI CONTROL :: // SSTI CONTROL: // RESTARTEVERYTIME (default: 0) "write restart possibility every RESTARTEVERY steps" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: -1) "time step size dt" // RESULTSEVERYTIME (default: 0) "increment for writing solution" // RESULTSEVERY (default: 1) "increment for writing solution" // ITEMAX (default: 10) "maximum number of iterations over fields" // SCATRA_FROM_RESTART_FILE (default: false) "read scatra result from restart files (use option 'restartfromfile' during execution of 4C)" // SCATRA_FILENAME (default: nil) "Control-file name for reading scatra results in SSTI" // COUPALGO (default: ssti_Monolithic) (choices: ssti_Monolithic) "Coupling strategies for SSTI solvers" // SCATRATIMINTTYPE (default: Elch) (choices: Elch) "scalar transport time integration type is needed to instantiate correct scalar transport time integration scheme for ssi problems" // ADAPTIVE_TIMESTEPPING (default: false) "flag for adaptive time stepping" .. _SECssticontrol_monolithic: SSTI CONTROL/MONOLITHIC ----------------------- SSTI CONTROL/MONOLITHIC :: // SSTI CONTROL/MONOLITHIC: // ABSTOLRES (default: 1e-14) "absolute tolerance for deciding if global residual of nonlinear problem is already zero" // CONVTOL (default: 1e-06) "tolerance for convergence check of Newton-Raphson iteration within monolithic SSI" // LINEAR_SOLVER (default: -1) "ID of linear solver for global system of equations" // MATRIXTYPE (default: undefined) (choices: block|sparse|undefined) "type of global system matrix in global system of equations" // EQUILIBRATION (default: none) "flag for equilibration of global system of equations" // EQUILIBRATION_STRUCTURE (default: none) "flag for equilibration of structural equations" // EQUILIBRATION_SCATRA (default: none) "flag for equilibration of scatra equations" // EQUILIBRATION_THERMO (default: none) "flag for equilibration of scatra equations" // EQUILIBRATION_INIT_SCATRA (default: false) "use equilibration method of ScaTra to equilibrate initial calculation of potential" .. _SECssticontrol_thermo: SSTI CONTROL/THERMO ------------------- SSTI CONTROL/THERMO :: // SSTI CONTROL/THERMO: // INITTHERMOFUNCT (default: -1) "initial function for thermo field" // LINEAR_SOLVER (default: -1) "linear solver for thermo field" // INITIALFIELD (default: field_by_function) (choices: field_by_condition|field_by_function) "defines, how to set the initial field" .. _SECstidynamic: STI DYNAMIC ----------- STI DYNAMIC :: // STI DYNAMIC: // SCATRATIMINTTYPE (default: Standard) (choices: Elch|Standard) "scalar transport time integration type is needed to instantiate correct scalar transport time integration scheme for scatra-thermo interaction problems" // COUPLINGTYPE (default: Undefined) (choices: Monolithic|OneWay_ScatraToThermo|OneWay_ThermoToScatra|TwoWay_ScatraToThermo|TwoWay_ScatraToThermo_Aitken|TwoWay_ScatraToThermo_Aitken_Dofsplit|TwoWay_ThermoToScatra|TwoWay_ThermoToScatra_Aitken|Undefined) "type of coupling between scatra and thermo fields" // THERMO_INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "initial temperature field for scatra-thermo interaction problems" // THERMO_INITFUNCNO (default: -1) "function number for initial temperature field for scatra-thermo interaction problems" // THERMO_LINEAR_SOLVER (default: -1) "ID of linear solver for temperature field" // THERMO_CONDENSATION (default: false) "flag for double condensation of linear equations associated with temperature field" .. _SECstidynamic_monolithic: STI DYNAMIC/MONOLITHIC ---------------------- STI DYNAMIC/MONOLITHIC :: // STI DYNAMIC/MONOLITHIC: // LINEAR_SOLVER (default: -1) "ID of linear solver for global system of equations" // MATRIXTYPE (default: block) (choices: block|sparse) "type of global system matrix in global system of equations" .. _SECstidynamic_partitioned: STI DYNAMIC/PARTITIONED ----------------------- STI DYNAMIC/PARTITIONED :: // STI DYNAMIC/PARTITIONED: // OMEGA (default: 1) "relaxation parameter" // OMEGAMAX (default: 0) "maximum value of Aitken relaxation parameter (0.0 = no constraint)" .. _SECstructnox: STRUCT NOX ---------- STRUCT NOX :: // STRUCT NOX: // Nonlinear Solver (default: Line Search Based) (choices: Inexact Trust Region Based|Line Search Based|Pseudo Transient|Single Step|Tensor Based|Trust Region Based) "Choose a nonlinear solver method." .. _SECstructnox_direction: STRUCT NOX/Direction -------------------- STRUCT NOX/Direction :: // STRUCT NOX/Direction: // Method (default: Newton) (choices: Broyden|Newton|NonlinearCG|Steepest Descent|User Defined) "Choose a direction method for the nonlinear solver." // User Defined Method (default: Modified Newton) (choices: Modified Newton|Newton) "Choose a user-defined direction method." .. _SECstructnox_direction_newton: STRUCT NOX/Direction/Newton --------------------------- STRUCT NOX/Direction/Newton :: // STRUCT NOX/Direction/Newton: // Forcing Term Method (default: Constant) (choices: Constant|Type 1|Type 2) // Forcing Term Initial Tolerance (default: 0.1) "initial linear solver tolerance" // Forcing Term Minimum Tolerance (default: 1e-06) // Forcing Term Maximum Tolerance (default: 0.01) // Forcing Term Alpha (default: 1.5) "used only by \"Type 2\"" // Forcing Term Gamma (default: 0.9) "used only by \"Type 2\"" // Rescue Bad Newton Solve (default: true) "If set to true, we will use the computed direction even if the linear solve does not achieve the tolerance specified by the forcing term" .. _SECstructnox_direction_newton_linearsolver: STRUCT NOX/Direction/Newton/Linear Solver ----------------------------------------- STRUCT NOX/Direction/Newton/Linear Solver :: // STRUCT NOX/Direction/Newton/Linear Solver: // Adaptive Control (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // Adaptive Control Objective (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // Zero Initial Guess (default: true) "Zero out the delta X vector if requested." // Computing Scaling Manually (default: false) "Allows the manually scaling of your linear system (not supported at the moment)." // Output Solver Details (default: true) "Switch the linear solver output on and off." .. _SECstructnox_direction_steepestdescent: STRUCT NOX/Direction/Steepest Descent ------------------------------------- STRUCT NOX/Direction/Steepest Descent :: // STRUCT NOX/Direction/Steepest Descent: // Scaling Type (default: None) (choices: 2-Norm|F 2-Norm|None|Quadratic Model Min) .. _SECstructnox_linesearch: STRUCT NOX/Line Search ---------------------- STRUCT NOX/Line Search :: // STRUCT NOX/Line Search: // Method (default: Full Step) (choices: Backtrack|Full Step|More'-Thuente|Polynomial|User Defined) // Inner Status Test Check Type (default: Minimal) (choices: Complete|Minimal|None) "Specify the check type for the inner status tests." .. _SECstructnox_linesearch_backtrack: STRUCT NOX/Line Search/Backtrack -------------------------------- STRUCT NOX/Line Search/Backtrack :: // STRUCT NOX/Line Search/Backtrack: // Default Step (default: 1) "starting step length" // Minimum Step (default: 1e-12) "minimum acceptable step length" // Recovery Step (default: 1) "step to take when the line search fails (defaults to value for \"Default Step\")" // Max Iters (default: 50) "maximum number of iterations" // Reduction Factor (default: 0.5) "A multiplier between zero and one that reduces the step size between line search iterations" // Allow Exceptions (default: false) "Set to true, if exceptions during the force evaluation and backtracking routine should be allowed." .. _SECstructnox_linesearch_fullstep: STRUCT NOX/Line Search/Full Step -------------------------------- STRUCT NOX/Line Search/Full Step :: // STRUCT NOX/Line Search/Full Step: // Full Step (default: 1) "length of a full step" .. _SECstructnox_linesearch_more'-thuente: STRUCT NOX/Line Search/More'-Thuente ------------------------------------ STRUCT NOX/Line Search/More'-Thuente :: // STRUCT NOX/Line Search/More'-Thuente: // Sufficient Decrease (default: 0.0001) "The ftol in the sufficient decrease condition" // Curvature Condition (default: 0.9999) "The gtol in the curvature condition" // Interval Width (default: 1e-15) "The maximum width of the interval containing the minimum of the modified function" // Maximum Step (default: 1e+06) "maximum allowable step length" // Minimum Step (default: 1e-12) "minimum allowable step length" // Max Iters (default: 20) "maximum number of right-hand-side and corresponding Jacobian evaluations" // Default Step (default: 1) "starting step length" // Recovery Step Type (default: Constant) (choices: Constant|Last Computed Step) "Determines the step size to take when the line search fails" // Recovery Step (default: 1) "The value of the step to take when the line search fails. Only used if the \"Recovery Step Type\" is set to \"Constant\"" // Sufficient Decrease Condition (default: Armijo-Goldstein) (choices: Ared/Pred|Armijo-Goldstein|None) "Choice to use for the sufficient decrease condition" // Optimize Slope Calculation (default: false) "Boolean value. If set to true the value of $s^T J^T F$ is estimated using a directional derivative in a call to ::NOX::LineSearch::Common::computeSlopeWithOutJac. If false the slope computation is computed with the ::NOX::LineSearch::Common::computeSlope method. Setting this to true eliminates having to compute the Jacobian at each inner iteration of the More'-Thuente line search" .. _SECstructnox_linesearch_polynomial: STRUCT NOX/Line Search/Polynomial --------------------------------- STRUCT NOX/Line Search/Polynomial :: // STRUCT NOX/Line Search/Polynomial: // Default Step (default: 1) "Starting step length" // Max Iters (default: 100) "Maximum number of line search iterations. The search fails if the number of iterations exceeds this value" // Minimum Step (default: 1e-12) "Minimum acceptable step length. The search fails if the computed $\\lambda_k$ is less than this value" // Recovery Step Type (default: Constant) (choices: Constant|Last Computed Step) "Determines the step size to take when the line search fails" // Recovery Step (default: 1) "The value of the step to take when the line search fails. Only used if the \"Recovery Step Type\" is set to \"Constant\"" // Interpolation Type (default: Cubic) (choices: Cubic|Quadratic|Quadratic3) "Type of interpolation that should be used" // Min Bounds Factor (default: 0.1) "Choice for $\\gamma_{\\min}$, i.e., the factor that limits the minimum size of the new step based on the previous step" // Max Bounds Factor (default: 0.5) "Choice for $\\gamma_{\\max}$, i.e., the factor that limits the maximum size of the new step based on the previous step" // Sufficient Decrease Condition (default: Armijo-Goldstein) (choices: Ared/Pred|Armijo-Goldstein|None) "Choice to use for the sufficient decrease condition" // Alpha Factor (default: 0.0001) "Parameter choice for sufficient decrease condition" // Force Interpolation (default: false) "Set to true if at least one interpolation step should be used. The default is false which means that the line search will stop if the default step length satisfies the convergence criteria" // Use Counters (default: true) "Set to true if we should use counters and then output the result to the parameter list as described in Output Parameters" // Maximum Iteration for Increase (default: 0) "Maximum index of the nonlinear iteration for which we allow a relative increase" // Allowed Relative Increase (default: 100) .. _SECstructnox_printing: STRUCT NOX/Printing ------------------- STRUCT NOX/Printing :: // STRUCT NOX/Printing: // Error (default: false) // Warning (default: true) // Outer Iteration (default: true) // Inner Iteration (default: true) // Parameters (default: false) // Details (default: false) // Outer Iteration StatusTest (default: true) // Linear Solver Details (default: false) // Test Details (default: false) // Debug (default: false) .. _SECstructnox_pseudotransient: STRUCT NOX/Pseudo Transient --------------------------- STRUCT NOX/Pseudo Transient :: // STRUCT NOX/Pseudo Transient: // deltaInit (default: -1) "Initial time step size. If its negative, the initial time step is calculated automatically." // deltaMax (default: 1.79769e+308) "Maximum time step size. If the new step size is greater than this value, the transient terms will be eliminated from the Newton iteration resulting in a full Newton solve." // deltaMin (default: 1e-05) "Minimum step size." // Max Number of PTC Iterations (default: 2147483647) // SER_alpha (default: 1) "Exponent of SET." // ScalingFactor (default: 1) "Scaling Factor for ptc matrix." // Time Step Control (default: SER) (choices: MRR|Model Reduction Ratio|SER|Switched Evolution Relaxation|TTE|Temporal Truncation Error) // Norm Type for TSC (default: Max Norm) (choices: Max Norm|One Norm|Two Norm) "Norm Type for the time step control" // Scaling Type (default: Identity) (choices: CFL Diagonal|Element based|Identity|Lumped Mass) "Type of the scaling matrix for the PTC method." // Build scaling operator (default: every timestep) (choices: every iter|every timestep) "Build scaling operator in every iteration or timestep" .. _SECstructnox_solveroptions: STRUCT NOX/Solver Options ------------------------- STRUCT NOX/Solver Options :: // STRUCT NOX/Solver Options: // Merit Function (default: Sum of Squares) (choices: Sum of Squares) // Status Test Check Type (default: Complete) (choices: Complete|Minimal|None) .. _SECstructnox_statustest: STRUCT NOX/Status Test ---------------------- STRUCT NOX/Status Test :: // STRUCT NOX/Status Test: // XML File > (default: none) "Filename of XML file with configuration of nox status test" .. _SECstructnox_trustregion: STRUCT NOX/Trust Region ----------------------- STRUCT NOX/Trust Region :: // STRUCT NOX/Trust Region: // Minimum Trust Region Radius (default: 1e-06) "Minimum allowable trust region radius" // Maximum Trust Region Radius (default: 1e+09) "Maximum allowable trust region radius" // Minimum Improvement Ratio (default: 0.0001) "Minimum improvement ratio to accept the step" // Contraction Trigger Ratio (default: 0.1) "If the improvement ratio is less than this value, then the trust region is contracted by the amount specified by the \"Contraction Factor\". Must be larger than \"Minimum Improvement Ratio\"" // Contraction Factor (default: 0.25) // Expansion Trigger Ratio (default: 0.75) "If the improvement ratio is greater than this value, then the trust region is contracted by the amount specified by the \"Expansion Factor\"" // Expansion Factor (default: 4) // Recovery Step (default: 1) .. _SECstructuraldynamic: STRUCTURAL DYNAMIC ------------------ STRUCTURAL DYNAMIC :: // STRUCTURAL DYNAMIC: // INT_STRATEGY (default: Old) (choices: Old|Standard) "global type of the used integration strategy" // TIME_ADAPTIVITY (default: false) "Enable adaptive time integration" // DYNAMICTYPE (default: GenAlpha) (choices: AdamsBashforth2|AdamsBashforth4|CentrDiff|ExplicitEuler|GenAlpha|GenAlphaLieGroup|OneStepTheta|Statics) "type of the specific dynamic time integration scheme" // PRESTRESS (default: NONE) (choices: MATERIAL_ITERATIVE|MULF|Material_Iterative|Mulf|NONE|None|material_iterative|mulf|none) "prestressing takes values none mulf material_iterative" // PRESTRESSTIME (default: 0) "time to switch from pre to post stressing" // PRESTRESSTOLDISP (default: 1e-09) "tolerance in the displacement norm during prestressing" // PRESTRESSMINLOADSTEPS (default: 0) "Minimum number of load steps during prestressing" // RESULTSEVERY (default: 1) "save displacements and contact forces every RESULTSEVERY steps" // RESEVERYERGY (default: 0) "write system energies every requested step" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // CALC_ACC_ON_RESTART (default: false) "Compute the initial state for a restart dynamics analysis" // OUTPUT_STEP_OFFSET (default: 0) "An offset added to the current step to shift the steps to be written." // TIMESTEP (default: 0.05) "time step size" // NUMSTEP (default: 200) "maximum number of steps" // TIMEINIT (default: 0) "initial time" // MAXTIME (default: 5) "maximum time" // DAMPING (default: None) (choices: Material|None|Rayleigh) "type of damping: (1) Rayleigh damping matrix and use it from M_DAMP x M + K_DAMP x K, (2) Material based and calculated in elements" // M_DAMP (default: -1) // K_DAMP (default: -1) // TOLDISP (default: 1e-10) "tolerance in the displacement norm for the newton iteration" // NORM_DISP (default: Abs) (choices: Abs|Mix|Rel) "type of norm for displacement convergence check" // TOLRES (default: 1e-08) "tolerance in the residual norm for the newton iteration" // NORM_RESF (default: Abs) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // TOLPRE (default: 1e-08) "tolerance in pressure norm for the newton iteration" // NORM_PRES (default: Abs) (choices: Abs) "type of norm for pressure convergence check" // TOLINCO (default: 1e-08) "tolerance in the incompressible residual norm for the newton iteration" // NORM_INCO (default: Abs) (choices: Abs) "type of norm for incompressible residual convergence check" // NORMCOMBI_DISPPRES (default: And) (choices: And|Or) "binary operator to combine pressure and displacement values" // NORMCOMBI_RESFINCO (default: And) (choices: And|Or) "binary operator to combine force and incompressible residual" // NORMCOMBI_RESFDISP (default: And) (choices: And|Or) "binary operator to combine displacement and residual force values" // STC_SCALING (default: Inactive) (choices: Inactive|Right|Symmetric) "Scaled director conditioning for thin shell structures" // STC_LAYER (default: 1) "number of STC layers for multilayer case" // PTCDT (default: 0.1) "pseudo time step for pseudo transient continuation (PTC) stabilized Newton procedure" // TOLCONSTR (default: 1e-08) "tolerance in the constr error norm for the newton iteration" // TOLCONSTRINCR (default: 1e-08) "tolerance in the constr lm incr norm for the newton iteration" // MAXITER (default: 50) "maximum number of iterations allowed for Newton-Raphson iteration before failure" // MINITER (default: 0) "minimum number of iterations to be done within Newton-Raphson loop" // ITERNORM (default: L2) (choices: Inf|L1|L2|Rms) "type of norm to be applied to residuals" // DIVERCONT (default: stop) (choices: adapt_3D0Dptc_ele_err|adapt_penaltycontact|adapt_step|continue|halve_step|rand_adapt_step|rand_adapt_step_ele_err|repeat_simulation|repeat_step|stop) "What to do with time integration when Newton-Raphson iteration failed" // MAXDIVCONREFINEMENTLEVEL (default: 10) "number of times timestep is halved in case nonlinear solver diverges" // NLNSOL (default: fullnewton) (choices: NoxNewtonLineSearch|augmentedlagrange|fullnewton|lsnewton|modnewton|newtonlinuzawa|noxgeneral|noxnln|ptc|singlestep|vague) "Nonlinear solution technique" // LSMAXITER (default: 30) "maximum number of line search steps" // ALPHA_LS (default: 0.5) "step reduction factor alpha in (Newton) line search scheme" // SIGMA_LS (default: 0.0001) "sufficient descent factor in (Newton) line search scheme" // MATERIALTANGENT (default: analytical) (choices: analytical|finitedifferences) "way of evaluating the constitutive matrix" // LOADLIN (default: false) "Use linearization of external follower load in Newton" // MASSLIN (default: none) (choices: none|rotations) "Application of nonlinear inertia terms" // NEGLECTINERTIA (default: false) "Neglect inertia" // PREDICT (default: ConstDis) (choices: ConstAcc|ConstDis|ConstDisPres|ConstDisVelAcc|ConstDisVelAccPres|ConstVel|TangDis|TangDisConstFext|Vague) "Type of predictor" // UZAWAPARAM (default: 1) "Parameter for Uzawa algorithm dealing with lagrange multipliers" // UZAWATOL (default: 1e-08) "Tolerance for iterative solve with Uzawa algorithm" // UZAWAMAXITER (default: 50) "maximum number of iterations allowed for uzawa algorithm before failure going to next newton step" // UZAWAALGO (default: direct) (choices: direct|simple|uzawa) // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // LUMPMASS (default: false) "Lump the mass matrix for explicit time integration" // MODIFIEDEXPLEULER (default: true) "Use the modified explicit Euler time integration scheme" // LINEAR_SOLVER (default: -1) "number of linear solver used for structural problems" // MIDTIME_ENERGY_TYPE (default: vague) (choices: imrLike|trLike|vague) "Specify the mid-averaging type for the structural energy contributions" // INITIALDISP (default: zero_displacement) (choices: displacement_by_function|zero_displacement) "Initial displacement for structure problem" // STARTFUNCNO (default: -1) "Function for Initial displacement" .. _SECstructuraldynamic_errorevaluation: STRUCTURAL DYNAMIC/ERROR EVALUATION ----------------------------------- STRUCTURAL DYNAMIC/ERROR EVALUATION :: // STRUCTURAL DYNAMIC/ERROR EVALUATION: // EVALUATE_ERROR_ANALYTICAL_REFERENCE (default: false) "Calculate error with respect to analytical solution defined by a function" // ANALYTICAL_DISPLACEMENT_FUNCTION (default: -1) "function ID of the analytical solution" .. _SECstructuraldynamic_genalpha: STRUCTURAL DYNAMIC/GENALPHA --------------------------- STRUCTURAL DYNAMIC/GENALPHA :: // STRUCTURAL DYNAMIC/GENALPHA: // GENAVG (default: TrLike) (choices: ImrLike|TrLike|Vague) "mid-average type of internal forces" // BETA (default: -1) "Generalised-alpha factor in (0,1/2]" // GAMMA (default: -1) "Generalised-alpha factor in (0,1]" // ALPHA_M (default: -1) "Generalised-alpha factor in [0,1)" // ALPHA_F (default: -1) "Generalised-alpha factor in [0,1)" // RHO_INF (default: 1) "Spectral radius for generalised-alpha time integration, valid range is [0,1]" .. _SECstructuraldynamic_onesteptheta: STRUCTURAL DYNAMIC/ONESTEPTHETA ------------------------------- STRUCTURAL DYNAMIC/ONESTEPTHETA :: // STRUCTURAL DYNAMIC/ONESTEPTHETA: // THETA (default: 0.5) "One-step-theta factor in (0,1]" .. _SECstructuraldynamic_timeadaptivity: STRUCTURAL DYNAMIC/TIMEADAPTIVITY --------------------------------- STRUCTURAL DYNAMIC/TIMEADAPTIVITY :: // STRUCTURAL DYNAMIC/TIMEADAPTIVITY: // KIND (default: None) (choices: AdamsBashforth2|CentralDifference|ExplicitEuler|JointExplicit|None|ZienkiewiczXie) "Method for time step size adaptivity" // OUTSYSPERIOD (default: 0) "Write system vectors (displacements, velocities, etc) every given period of time" // OUTSTRPERIOD (default: 0) "Write stress/strain every given period of time" // OUTENEPERIOD (default: 0) "Write energy every given period of time" // OUTRESTPERIOD (default: 0) "Write restart data every given period of time" // OUTSIZEEVERY (default: 0) "Write step size every given time step" // STEPSIZEMAX (default: 0) "Limit maximally permitted time step size (>0)" // STEPSIZEMIN (default: 0) "Limit minimally allowed time step size (>0)" // SIZERATIOMAX (default: 0) "Limit maximally permitted change of time step size compared to previous size, important for multi-step schemes (>0)" // SIZERATIOMIN (default: 0) "Limit minimally permitted change of time step size compared to previous size, important for multi-step schemes (>0)" // SIZERATIOSCALE (default: 0.9) "This is a safety factor to scale theoretical optimal step size, should be lower than 1 and must be larger than 0" // LOCERRNORM (default: Vague) (choices: Inf|L1|L2|Rms|Vague) "Vector norm to treat error vector with" // LOCERRTOL (default: 0) "Target local error tolerance (>0)" // ADAPTSTEPMAX (default: 0) "Limit maximally allowed step size reduction attempts (>0)" .. _SECstructuraldynamic_timeadaptivity_jointexplicit: STRUCTURAL DYNAMIC/TIMEADAPTIVITY/JOINT EXPLICIT ------------------------------------------------ STRUCTURAL DYNAMIC/TIMEADAPTIVITY/JOINT EXPLICIT :: // STRUCTURAL DYNAMIC/TIMEADAPTIVITY/JOINT EXPLICIT: // LINEAR_SOLVER (default: -1) "number of linear solver used for auxiliary integrator" // INT_STRATEGY (default: Standard) (choices: Standard) "global type of the used integration strategy" // DYNAMICTYPE (default: CentrDiff) (choices: AdamsBashforth2|AdamsBashforth4|CentrDiff|ExplicitEuler) "type of the specific auxiliary dynamic time integration scheme" // LUMPMASS (default: false) "Lump the mass matrix for explicit time integration" // DAMPING (default: None) (choices: Material|None|Rayleigh) "type of damping: (1) Rayleigh damping matrix and use it from M_DAMP x M + K_DAMP x K, (2) Material based and calculated in elements" // M_DAMP (default: -1) // K_DAMP (default: -1) .. _SECstructuregeometry: STRUCTURE GEOMETRY ------------------ STRUCTURE GEOMETRY :: // STRUCTURE GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECthermaldynamic: THERMAL DYNAMIC --------------- THERMAL DYNAMIC :: // THERMAL DYNAMIC: // DYNAMICTYPE (default: OneStepTheta) "type of time integration control" // RESULTSEVERY (default: 1) "save temperature and other global quantities every RESULTSEVERY steps" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // INITIALFIELD (default: zero_field) (choices: field_by_condition|field_by_function|zero_field) "Initial Field for thermal problem" // INITFUNCNO (default: -1) "function number for thermal initial field" // TIMESTEP (default: 0.05) "time step size" // NUMSTEP (default: 200) "maximum number of steps" // MAXTIME (default: 5) "maximum time" // TOLTEMP (default: 1e-10) "tolerance in the temperature norm of the Newton iteration" // NORM_TEMP (default: Abs) (choices: Abs|Mix|Rel) "type of norm for temperature convergence check" // TOLRES (default: 1e-08) "tolerance in the residual norm for the Newton iteration" // NORM_RESF (default: Abs) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // NORMCOMBI_RESFTEMP (default: And) (choices: And|Or) "binary operator to combine temperature and residual force values" // MAXITER (default: 50) "maximum number of iterations allowed for Newton-Raphson iteration before failure" // MINITER (default: 0) "minimum number of iterations to be done within Newton-Raphson loop" // ITERNORM (default: L2) (choices: Inf|L1|L2|Rms) "type of norm to be applied to residuals" // DIVERCONT (default: stop) (choices: continue|halve_step|repeat_simulation|repeat_step|stop) "What to do with time integration when Newton-Raphson iteration failed" // MAXDIVCONREFINEMENTLEVEL (default: 10) "number of times timestep is halved in case nonlinear solver diverges" // NLNSOL (default: fullnewton) (choices: fullnewton|vague) "Nonlinear solution technique" // PREDICT (default: ConstTemp) (choices: ConstTemp|ConstTempRate|TangTemp|Vague) "Predictor of iterative solution techniques" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // LUMPCAPA (default: false) "Lump the capacity matrix for explicit time integration" // LINEAR_SOLVER (default: -1) "number of linear solver used for thermal problems" // CALCERROR (default: No) (choices: No|byfunct) "compute error compared to analytical solution" // CALCERRORFUNCNO (default: -1) "Function for Error Calculation" .. _SECthermaldynamic_genalpha: THERMAL DYNAMIC/GENALPHA ------------------------ THERMAL DYNAMIC/GENALPHA :: // THERMAL DYNAMIC/GENALPHA: // GENAVG (default: TrLike) (choices: ImrLike|TrLike|Vague) "mid-average type of internal forces" // GAMMA (default: 0.5) "Generalised-alpha factor in (0,1]" // ALPHA_M (default: 0.5) "Generalised-alpha factor in [0.5,1)" // ALPHA_F (default: 0.5) "Generalised-alpha factor in [0.5,1)" // RHO_INF (default: -1) "Generalised-alpha factor in [0,1]" .. _SECthermaldynamic_onesteptheta: THERMAL DYNAMIC/ONESTEPTHETA ---------------------------- THERMAL DYNAMIC/ONESTEPTHETA :: // THERMAL DYNAMIC/ONESTEPTHETA: // THETA (default: 0.5) "One-step-theta factor in (0,1]" .. _SECthermaldynamic_runtimecsvoutput: THERMAL DYNAMIC/RUNTIME CSV OUTPUT ---------------------------------- THERMAL DYNAMIC/RUNTIME CSV OUTPUT :: // THERMAL DYNAMIC/RUNTIME CSV OUTPUT: // OUTPUT_THERMO (default: false) "write thermo output" // ENERGY (default: false) "write energy output" .. _SECthermaldynamic_runtimevtkoutput: THERMAL DYNAMIC/RUNTIME VTK OUTPUT ---------------------------------- THERMAL DYNAMIC/RUNTIME VTK OUTPUT :: // THERMAL DYNAMIC/RUNTIME VTK OUTPUT: // OUTPUT_THERMO (default: false) "write thermo output" // TEMPERATURE (default: false) "write temperature output" // HEATFLUX (default: false) "write heatflux output" // TEMPGRAD (default: false) "write temperature gradient output" // ELEMENT_OWNER (default: false) "write element owner" // ELEMENT_GID (default: false) "write 4C internal element GIDs" // NODE_GID (default: false) "write 4C internal node GIDs" .. _SECthermogeometry: THERMO GEOMETRY --------------- THERMO GEOMETRY :: // THERMO GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECtransportgeometry: TRANSPORT GEOMETRY ------------------ TRANSPORT GEOMETRY :: // TRANSPORT GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECtransport2geometry: TRANSPORT2 GEOMETRY ------------------- TRANSPORT2 GEOMETRY :: // TRANSPORT2 GEOMETRY: // FILE "Path to the exodus geometry file. Either absolute or relative to the input file." // list 'ELEMENT_BLOCKS' with entries: // : // ID "ID of the element block in the exodus file." // ELEMENT_NAME "The name of the element that should be assigned to the block." // ELEMENT_DATA "A dat-style string of parameters for the element." .. _SECtsicontact: TSI CONTACT ----------- TSI CONTACT :: // TSI CONTACT: // HEATTRANSSLAVE (default: 0) "Heat transfer parameter for slave side in thermal contact" // HEATTRANSMASTER (default: 0) "Heat transfer parameter for master side in thermal contact" // TEMP_DAMAGE (default: 1e+12) "damage temperature at contact interface: friction coefficient zero there" // TEMP_REF (default: 0) "reference temperature at contact interface: friction coefficient equals the given value" // NITSCHE_THETA_TSI (default: 0) "+1: symmetric, 0: non-symmetric, -1: skew-symmetric" // NITSCHE_WEIGHTING_TSI (default: harmonic) "how to weight consistency terms in Nitsche contact formulation" // NITSCHE_PENALTY_ADAPTIVE_TSI (default: true) "adapt penalty parameter after each converged time step" // PENALTYPARAM_THERMO (default: 0) "Penalty parameter for Nitsche solution strategy" .. _SECtsidynamic: TSI DYNAMIC ----------- TSI DYNAMIC :: // TSI DYNAMIC: // COUPALGO (default: tsi_monolithic) (choices: tsi_iterstagg|tsi_iterstagg_aitken|tsi_iterstagg_aitkenirons|tsi_iterstagg_fixedrelax|tsi_monolithic|tsi_oneway|tsi_sequstagg) "Coupling strategies for TSI solvers" // MATCHINGGRID (default: true) "is matching grid" // RESTARTEVERY (default: 1) "write restart possibility every RESTARTEVERY steps" // NUMSTEP (default: 200) "maximum number of Timesteps" // MAXTIME (default: 1000) "total simulation time" // TIMESTEP (default: 0.05) "time step size dt" // ITEMAX (default: 10) "maximum number of iterations over fields" // ITEMIN (default: 1) "minimal number of iterations over fields" // RESULTSEVERY (default: 1) "increment for writing solution" // NORM_INC (default: Abs) (choices: Abs|Mix|Rel) "type of norm for convergence check of primary variables in TSI" .. _SECtsidynamic_monolithic: TSI DYNAMIC/MONOLITHIC ---------------------- TSI DYNAMIC/MONOLITHIC :: // TSI DYNAMIC/MONOLITHIC: // CONVTOL (default: 1e-06) "tolerance for convergence check of TSI" // TOLINC (default: 1e-06) "tolerance for convergence check of TSI-increment in monolithic TSI" // NORM_RESF (default: Abs) (choices: Abs|Mix|Rel) "type of norm for residual convergence check" // NORMCOMBI_RESFINC (default: Coupl_And_Single) (choices: And|And_Single|Coupl_And_Single|Coupl_Or_Single|Or|Or_Single) "binary operator to combine primary variables and residual force values" // ITERNORM (default: Rms) (choices: Inf|L1|L1_Scaled|L2|Rms) "type of norm to be applied to residuals" // NLNSOL (default: fullnewton) (choices: fullnewton|ptc) "Nonlinear solution technique" // PTCDT (default: 0.1) "pseudo time step for pseudo-transient continuation (PTC) stabilised Newton procedure" // LINEAR_SOLVER (default: -1) "number of linear solver used for monolithic TSI problems" // ADAPTCONV (default: false) "Switch on adaptive control of linear solver tolerance for nonlinear solution" // ADAPTCONV_BETTER (default: 0.1) "The linear solver shall be this much better than the current nonlinear residual in the nonlinear convergence limit" // INFNORMSCALING (default: true) "Scale blocks of matrix with row infnorm?" // MERGE_TSI_BLOCK_MATRIX (default: false) "Merge TSI block matrix" // CALC_NECKING_TSI_VALUES (default: false) "Calculate nodal values for evaluation and validation of necking" // TSI_LINE_SEARCH (default: none) (choices: and|none|or|structure|thermo) "line-search strategy" .. _SECtsidynamic_partitioned: TSI DYNAMIC/PARTITIONED ----------------------- TSI DYNAMIC/PARTITIONED :: // TSI DYNAMIC/PARTITIONED: // COUPVARIABLE (default: Displacement) (choices: Displacement|Temperature) "Coupling variable" // MAXOMEGA (default: 0) "largest omega allowed for Aitken relaxation (0.0 means no constraint)" // FIXEDOMEGA (default: 1) "fixed relaxation parameter" // CONVTOL (default: 1e-06) "tolerance for convergence check of outer iteraiton within partitioned TSI" .. _SECvolmortarcoupling: VOLMORTAR COUPLING ------------------ VOLMORTAR COUPLING :: // VOLMORTAR COUPLING: // INTTYPE (default: Elements) (choices: Elements|Segments|elements|segments) "Type of numerical integration scheme" // COUPLINGTYPE (default: Volmortar) (choices: Volmortar|consint|consistentinterpolation|volmortar) "Type of coupling" // SHAPEFCN (default: Dual) (choices: Dual|Standard|dual|standard|std) "Type of employed set of shape functions" // CUTTYPE (default: DirectDivergence) (choices: DirectDivergence|Tessellation|dd|directdivergence|t|tessellation) "Type of cut procedure/ integration point calculation" // DUALQUAD (default: nm) (choices: lin_mod|lm|nm|nomod|qm|quad_mod) "Type of dual shape function for weighting function for quadr. problems" // MESH_INIT (default: false) "If chosen, mesh initialization procedure is performed" // KEEP_EXTENDEDGHOSTING (default: true) "If chosen, extended ghosting is kept for simulation" .. _SECwear: WEAR ---- WEAR :: // WEAR: // WEARLAW (default: None) (choices: Archard|None|archard|none) "Type of wear law" // MATCHINGGRID (default: true) "is matching grid" // WEAR_SHAPEFCN (default: Standard) (choices: Dual|Standard|dual|standard|std) "Type of employed set of shape functions for wear" // WEARCOEFF (default: 0) "Wear coefficient for slave surface" // WEARCOEFF_MASTER (default: 0) "Wear coefficient for master surface" // WEAR_TIMERATIO (default: 1) "Time step ratio between wear and spatial time scale" // SSSLIP (default: 1) "Fixed slip for steady state wear" // SSWEAR (default: false) "flag for steady state wear" // WEAR_SIDE (default: Slave) (choices: Slave|both|s|slave|slave_master|sm) "Definition of wear side" // WEARTYPE (default: internal_state) (choices: internal_state|intstate|is|primary_variable|primvar|pv) "Definition of wear algorithm" // WEARTIMINT (default: e) (choices: e|expl|explicit|i|impl|implicit) "Definition of wear time integration" // WEAR_TIMESCALE (default: e) (choices: d|different|e|equal) "Definition wear time scale compares to std. time scale" .. _SECxfemgeneral: XFEM GENERAL ------------ XFEM GENERAL :: // XFEM GENERAL: // GMSH_DEBUG_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // GMSH_DEBUG_OUT_SCREEN (default: false) "Do you want to be informed, if Gmsh output is written?" // GMSH_SOL_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // GMSH_TIMINT_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // GMSH_EOS_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // GMSH_DISCRET_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // GMSH_CUT_OUT (default: false) "Do you want to write extended Gmsh output for each timestep?" // PRINT_OUTPUT (default: false) "Is the output of the cut process desired?" // MAX_NUM_DOFSETS (default: 3) "Maximum number of volumecells in the XFEM element" // NODAL_DOFSET_STRATEGY (default: full) (choices: ConnectGhostDofsets_PerNodeAndPosition|OneDofset_PerNodeAndPosition|full) "Strategy used for the nodal dofset management per node" // VOLUME_GAUSS_POINTS_BY (default: Tessellation) (choices: DirectDivergence|MomentFitting|Tessellation) "Method for finding Gauss Points for the cut volumes" // BOUNDARY_GAUSS_POINTS_BY (default: Tessellation) (choices: DirectDivergence|MomentFitting|Tessellation) "Method for finding Gauss Points for the boundary cells" .. _SECxfluiddynamic_general: XFLUID DYNAMIC/GENERAL ---------------------- XFLUID DYNAMIC/GENERAL :: // XFLUID DYNAMIC/GENERAL: // XFLUIDFLUID (default: false) "Use an embedded fluid patch." // RELAXING_ALE_EVERY (default: 1) "Relaxing Ale after how many monolithic steps" // RELAXING_ALE (default: true) "switch on/off for relaxing Ale in monolithic fluid-fluid-fsi" // XFLUIDFLUID_SEARCHRADIUS (default: 1) "Radius of the search tree" // MONOLITHIC_XFFSI_APPROACH (default: xffsi_fixedALE_partitioned) (choices: xffsi_fixedALE_interpolation|xffsi_fixedALE_partitioned|xffsi_full_newton) "The monolithic approach for xfluidfluid-fsi" // XFLUIDFLUID_TIMEINT (default: Xff_TimeInt_FullProj) (choices: Xff_TimeInt_FullProj|Xff_TimeInt_IncompProj|Xff_TimeInt_KeepGhostValues|Xff_TimeInt_ProjIfMoved) "The xfluidfluid-timeintegration approach" // XFLUID_TIMEINT (default: STD=COPY/SL_and_GHOST=COPY/GP) (choices: STD=COPY/PROJ_and_GHOST=COPY/PROJ/GP|STD=COPY/SL_and_GHOST=COPY/GP|STD=COPY_and_GHOST=COPY/GP|STD=SL(boundary-zone)_and_GHOST=GP) "The xfluid time integration approach" // ALE_XFluid (default: false) "XFluid is Ale Fluid?" // INTERFACE_TERMS_PREVIOUS_STATE (default: PreviousState_only_consistency) (choices: PreviousState_full|PreviousState_only_consistency) "how to treat interface terms from previous time step (new OST)" // XFLUID_TIMEINT_CHECK_INTERFACETIPS (default: true) "Xfluid TimeIntegration Special Check if node has changed the side!" // XFLUID_TIMEINT_CHECK_SLIDINGONSURFACE (default: true) "Xfluid TimeIntegration Special Check if node is sliding on surface!" .. _SECxfluiddynamic_stabilization: XFLUID DYNAMIC/STABILIZATION ---------------------------- XFLUID DYNAMIC/STABILIZATION :: // XFLUID DYNAMIC/STABILIZATION: // COUPLING_METHOD (default: Nitsche) (choices: Hybrid_LM_Cauchy_stress|Hybrid_LM_viscous_stress|Nitsche) "method how to enforce embedded boundary/coupling conditions at the interface" // HYBRID_LM_L2_PROJ (default: part_ele_proj) (choices: full_ele_proj|part_ele_proj) "perform the L2 projection between stress fields on whole element or on fluid part?" // VISC_ADJOINT_SYMMETRY (default: sym) (choices: no|none|skew|sym|yes) "viscous and adjoint viscous interface terms with matching sign?" // NIT_STAB_FAC (default: 35) "( stabilization parameter for Nitsche's penalty term" // NIT_STAB_FAC_TANG (default: 35) "( stabilization parameter for Nitsche's penalty tangential term" // VISC_STAB_TRACE_ESTIMATE (default: CT_div_by_hk) (choices: CT_div_by_hk|eigenvalue) "how to estimate the scaling from the trace inequality in Nitsche's method" // UPDATE_EIGENVALUE_TRACE_ESTIMATE (default: every_iter) (choices: every_iter|every_timestep|once) "how often should the local eigenvalue problem be updated" // VISC_STAB_HK (default: ele_vol_div_by_max_ele_surf) (choices: cut_vol_div_by_cut_surf|ele_vol_div_by_cut_surf|ele_vol_div_by_ele_surf|ele_vol_div_by_max_ele_surf|vol_equivalent) "how to define the characteristic element length in cut elements" // CONV_STAB_SCALING (default: none) (choices: abs_inflow|inflow|none) "scaling factor for viscous interface stabilization (Nitsche, MSH)" // XFF_CONV_STAB_SCALING (default: none) (choices: averaged|inflow|none) "scaling factor for convective interface stabilization of fluid-fluid Coupling" // MASS_CONSERVATION_COMBO (default: max) (choices: max|sum) "choose the maximum from viscous and convective contributions or just sum both up" // MASS_CONSERVATION_SCALING (default: only_visc) (choices: full|only_visc) "apply additional scaling of penalty term to enforce mass conservation for convection-dominated flow" // GHOST_PENALTY_STAB (default: false) "switch on/off ghost penalty interface stabilization" // GHOST_PENALTY_TRANSIENT_STAB (default: false) "switch on/off ghost penalty transient interface stabilization" // GHOST_PENALTY_2nd_STAB (default: false) "switch on/off ghost penalty interface stabilization for 2nd order derivatives" // GHOST_PENALTY_2nd_STAB_NORMAL (default: false) "switch between ghost penalty interface stabilization for 2nd order derivatives in normal or all spatial directions" // GHOST_PENALTY_FAC (default: 0.1) "define stabilization parameter ghost penalty interface stabilization" // GHOST_PENALTY_TRANSIENT_FAC (default: 0.001) "define stabilization parameter ghost penalty transient interface stabilization" // GHOST_PENALTY_2nd_FAC (default: 0.05) "define stabilization parameter ghost penalty 2nd order viscous interface stabilization" // GHOST_PENALTY_PRESSURE_2nd_FAC (default: 0.05) "define stabilization parameter ghost penalty 2nd order pressure interface stabilization" // XFF_EOS_PRES_EMB_LAYER (default: false) "switch on/off edge-based pressure stabilization on interface-contributing elements of the embedded fluid" // IS_PSEUDO_2D (default: false) "modify viscous interface stabilization due to the vanishing polynomial in third dimension when using strong Dirichlet conditions to block polynomials in one spatial dimension" // GHOST_PENALTY_ADD_INNER_FACES (default: false) "Apply ghost penalty stabilization also for inner faces if this is possible due to the dofsets" .. _SECxfluiddynamic_xfpsimonolithic: XFLUID DYNAMIC/XFPSI MONOLITHIC ------------------------------- XFLUID DYNAMIC/XFPSI MONOLITHIC :: // XFLUID DYNAMIC/XFPSI MONOLITHIC: // ITEMIN (default: 1) "How many iterations are performed minimal" // ITEMAX_OUTER (default: 5) "How many outer iterations are performed maximal" // ND_NEWTON_DAMPING (default: false) "Activate Newton damping based on residual and increment" // ND_MAX_DISP_ITERINC (default: -1) "Maximal displacement increment to apply full newton --> otherwise damp newton" // ND_MAX_VEL_ITERINC (default: -1) "Maximal fluid velocity increment to apply full newton --> otherwise damp newton" // ND_MAX_PRES_ITERINC (default: -1) "Maximal fluid pressure increment to apply full newton --> otherwise damp newton" // ND_MAX_PVEL_ITERINC (default: -1) "Maximal porofluid velocity increment to apply full newton --> otherwise damp newton" // ND_MAX_PPRES_ITERINC (default: -1) "Maximal porofluid pressure increment to apply full newton --> otherwise damp newton" // CUT_EVALUATE_MINTOL (default: 0) "Minimal value of the maximal structural displacement for which the CUT is evaluate in this iteration!" // EXTRAPOLATE_TO_ZERO (default: false) "the extrapolation of the fluid stress in the contact zone is relaxed to zero after a certain distance" // CUT_EVALUATE_MINITER (default: 0) "Minimal number of nonlinear iterations, before the CUT is potentially not evaluated" // POROCONTACTFPSI_HFRACTION (default: 1) "factor of element size, when transition between FPSI and PSCI is started!" // POROCONTACTFPSI_FULLPCFRACTION (default: 0) "ration of gap/(POROCONTACTFPSI_HFRACTION*h) when full PSCI is started!" // USE_PORO_PRESSURE (default: true) "the extrapolation of the fluid stress in the contact zone is relaxed to zero after a certtain distance"