mp-heat-transfer: Parameter descriptionο
Contentsο
π· baseο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Sets the base name for the application that will be fed to the problem type. |
|
|
|
Defines the dimension of the problem |
|
|
|
Floating point number format. Currently, only βdoubleβ is explicitely instantiated. |
|
|
|
Defines the number of initial global refinements |
|
|
|
Set this parameter to true to list parameters in output |
|
|
|
Sets the verbosity level of the console output: 0: silent: for non-robust tests and benchmark runs; 1: minimal: for robust tests; 2: detailed; 3: full |
|
base: feο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Finite Element.FE_Q: hexahedral continuous finite element with polynomial degree p; FE_SimplexP: tetrahedral continuous finite element with polynomial degree p; FE_Q_iso_Q1: hexahedral continuous finite element with p subdivisions containing linear elements; FE_DGQ: hexahedral discontinuous finite element with polynomial degree p |
|
|
|
Defines the degree p of the finite element type. If βtypeβ is βFE_Q_iso_Q1β this parameter defines the number of subdivisions. |
π· time steppingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Defines the start time for the solution of the levelset problem |
|
|
|
Sets the end time for the solution of the levelset problem |
|
|
|
Sets the step size for time stepping. For non-uniform time stepping, this parameter determines the size of the first time step. |
|
|
|
Sets the maximum number of melt_pool steps |
|
|
|
Set an analytical function to determine the time step size. For the prediction of the new time increment, the old time is used. |
π· adaptive meshingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to activate adaptive meshing |
|
|
|
Set this parameter to true to not refine/coarsen along boundaries. |
|
|
|
Defines the (upper) percentage of elements that should be refined |
|
|
|
Defines the (lower) percentage of elements that should be coarsened |
|
|
|
Defines the number of maximum refinement steps one grid cell will be undergone. |
|
|
|
Defines the number of minimum refinement steps one grid cell will be undergone. |
|
|
|
Defines the number of initial refinements. |
|
|
|
Defines at every nth step the amr should be performed. |
|
|
|
Minimum number of cells that must be marked for refinement/coarsening before the mesh is updated. |
|
|
|
Minimum indicator value required for a cell to be considered for refinement. |
|
|
|
Set this parameter to true to average the contribututions to the same DoF coming from different cells during solution transfer. |
π· heatο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|||
|
|
|
Choose the heat operator implementation. Options: diffuse, cut |
|
|||
|
|
|
Set this parameter to true to enable time-dependent bc. |
|
|||
|
|||
|
|||
|
|
|
Sets the maximum verbosity level of the console output. |
|
|||
|
|||
|
heat: feο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Finite Element.FE_Q: hexahedral continuous finite element with polynomial degree p; FE_SimplexP: tetrahedral continuous finite element with polynomial degree p; FE_Q_iso_Q1: hexahedral continuous finite element with p subdivisions containing linear elements; FE_DGQ: hexahedral discontinuous finite element with polynomial degree p |
|
|
|
Defines the degree p of the finite element type. If βtypeβ is βFE_Q_iso_Q1β this parameter defines the number of subdivisions. |
heat: cutο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to βfalseβ to ignore the gas phase. |
|
|
|
Parameter for one step theta time integration. |
|
|
|
Set this parameter to true to consider the explicit symmetry term. Note: this parameter only applies if the setup is two-phase. |
|
heat: cut: stabilizationο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Nitsche stabilization parameter. |
|
heat: cut: stabilization: ghost-penaltyο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Mass matrix ghost-penalty parameter for degree 0. |
|
|
|
Mass matrix ghost-penalty parameter for degree 1. |
|
|
|
Mass matrix ghost-penalty parameter for degree 2. |
|
|
|
Stiffness matrix ghost-penalty parameter for degree 0. |
|
|
|
Stiffness matrix ghost-penalty parameter for degree 1. |
|
|
|
Stiffness matrix ghost-penalty parameter for degree 2. |
heat: diffuseο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Perform phase interpolation via the volumetric thermal capacity (product of density and capacity) instead of interpolating density and thermal capacity individually. |
heat: radiative boundary conditionο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Emissivity. |
|
|
|
Infinity temperature. |
heat: convective boundary conditionο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Convection coefficient. |
|
|
|
Infinity temperature. |
heat: nlsolveο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set the number of maximum nonlinear iterations with standard tolerances. |
|
|
|
Set the tolerance for the maximum allowed correction of the unknown field. |
|
|
|
Set the tolerance for the maximum allowed residual of the nonlinear system. |
|
|
|
Set the number of maximum nonlinear iterations with alternative tolerances. |
|
|
|
Set the alternative tolerance for the maximum allowed correction of the unknown field. |
|
|
|
Set the alternative tolerance for the maximum allowed residual of the nonlinear system. |
|
|
|
Set to one for detailed solver output. |
heat: linear solverο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter for choosing an iterative linear solver type. |
|
|
|
Set this parameter for choosing a preconditioner type. |
|
|
|
Set the maximum number of iterations for solving the linear system of equations. |
|
|
|
Set the relative tolerance for a successful solution of the linear system of equations. |
|
|
|
Set the absolute tolerance for a successful solution of the linear system of equations. |
|
|
|
Set this parameter if a matrix free solution procedure should be performed. |
|
|
|
Set the monitor type of the linear solver. |
heat: predictorο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Choose a predictor type: none: use old value as initial guess; zero: se zeros as initial guess; linear_extrapolation: calculate the predictor by a linear combination from the two old solution vectors; least_squares_projection: least squares projection (WIP) |
|
|
|
Choose the number of old solution vectors considered.This parameter is only relevant for least squares projection.For all other predictors, this parameter will be set appropriately. |
π· materialο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
If this parameter is initialized, the material parameters of the specified material will be used as template. Individual properties can be modified. However, be aware to put |
|
|||
|
|||
|
|||
|
|
|
Solidus temperature (K). |
|
|
|
Liquidus temperature (K). |
|
|
|
Function type for the apparent capacity method to model latent heat during solidification. constant: apparent capacity is constant between the solidus and liquidus temperature; qlq: apparent capacity is given by a quadratic/quadratic function of temperature between the solidus and liquidus temperature (default); poly4_bell: apparent capacity is given by a bell-shaped quartic polynomial function of temperature between the solidus and liquidus temperature. |
|
|
|
Latent heat of fusion (J/kg) |
|
|
|
Boiling temperature (K). |
|
|
|
Latent heat of evaporation (J/kg). |
|
|
|
Molar mass (mol/kg). |
|
|
|
Reference temperature of the specific enthalpy |
|
|
|
Choose how to interpolate the properties over the interface. sharp: properties jump at heaviside = 0.5; smooth: properties are smeared between the phases proportional to the heaviside (default); consistent_with_evaporation: same as βsmoothβ, but the density is interpolated proportional by the harmonic mean. |
|
|
|
Choose how to interpolate the properties over between the liquid and the solid phase. mushy_zone: solid and liquid properties are interpolated between the solidus and liquidus temperature (default); sharp: the solid and liquid properties jump at the melting point, which is set via the solidus temperature. |
material: gasο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
thermal conductivity of the gas phase |
|
|
|
specific heat capacity of the gas phase |
|
|
|
density of the gas phase |
|
|
|
dynamic viscosity of the gas phase |
material: liquidο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
thermal conductivity of the liquid phase |
|
|
|
specific heat capacity of the liquid phase |
|
|
|
density of the liquid phase |
|
|
|
dynamic viscosity of the liquid phase |
material: solidο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
thermal conductivity of the solid phase |
|
|
|
specific heat capacity of the solid phase |
|
|
|
density of the solid phase |
|
|
|
dynamic viscosity of the solid phase |
π· laserο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Laser model. analytical_temperature: see Mirkoohi et al. (2019); volumetric: volumetric heat source, the intensity is defined by ββintensity profileββ; interface_projection: projection-based regularized continuum surface flux in βdirectionβ, the intensity is defined by ββintensity profileββ; interface_projection_sharp: projection-based sharp surface flux in βdirectionβ, the intensity is defined by ββintensity profileββ; interface_projection_sharp_conforming: projection-based sharp surface flux in βdirectionβ on a conforming mesh, the intensity is defined by ββintensity profileββ; RTE: continuum surface flux projected using the radiative transport equation in βdirectionβ, supporting shadowing of undercuts, the intensity is defined by ββintensity profileββ; |
|
|
|
Laser intensity profile. uniform: note that the βpowerβ input is treated as the uniform power density in the whole domain; Gauss: Gaussian laser intensity shape with βradiusβ that retains the βpowerβ; Gusarov: see Gusarov et al. (2009); |
|
|
|
Laser power |
|
|
|
Temporal distribution of the laser power |
|
|
|
In case of time-dependent laser power: activation time of |
|
|
|
In case of time-dependent laser power: end time of |
|
|
|
Laser energy absorptivity of the gaseous part of the domain. |
|
|
|
Laser energy absorptivity of the liquid part of the domain. |
|
|
`` |
Center coordinates of the laser beam starting position on the interface melt/gas. |
|
|
|
Scan speed of the laser |
|
|
`` |
Direction of laser motion as a vector |
|
|
`` |
Laser beam direction. |
|
|
`` |
Axis around which the initial laser beam direction will be rotated. Relevant only in 3D. |
|
|
|
Rotation angle applied to the laser beam direction (in 3D about βbeam rotation axisβ following the right-hand rule; in 2D: as defined by the 2D rotation matrix |
|
|
|
Laser beam radius. |
|
|||
|
|||
|
laser: gusarovο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Reflectivity of the material. |
|
|
|
Extinction coefficient in [1/m]. |
|
|
|
Layer thickness |
laser: analyticalο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Ambient temperature in the inert gas. |
|
|
|
Maximum temperature arising in the melt pool. If this temperature is lower than the boiling temperature, this value is corrected to correspond to the boiling temperature + 500 K. |
|
|
|
This factor scales the analytical temperature field to be anisotropic. |
laser: dirac delta function approximationο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Choose how to smear a parameter over the interface. |
|
|
|
Choose if weights should be computed automatically. |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the heavy phase (level set = 1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the heavy liquid/solid phase (level set = 1). |
π· rteο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|||
|
|
|
Sets the maximum verbosity level of the console output. The maximum level with respect to the base value is decisive. |
|
|
|
Choose a predictor type. |
|
|
|
Chooses the formulation of the absorptivity coefficient |
|
|
|
Minimum value for absorptivity to ensure a non-singular matrix for RTE. |
|
|||
|
|||
|
rte: feο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Finite Element.FE_Q: hexahedral continuous finite element with polynomial degree p; FE_SimplexP: tetrahedral continuous finite element with polynomial degree p; FE_Q_iso_Q1: hexahedral continuous finite element with p subdivisions containing linear elements; FE_DGQ: hexahedral discontinuous finite element with polynomial degree p |
|
|
|
Defines the degree p of the finite element type. If βtypeβ is βFE_Q_iso_Q1β this parameter defines the number of subdivisions. |
rte: linear solverο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter for choosing an iterative linear solver type. |
|
|
|
Set this parameter for choosing a preconditioner type. |
|
|
|
Set the maximum number of iterations for solving the linear system of equations. |
|
|
|
Set the relative tolerance for a successful solution of the linear system of equations. |
|
|
|
Set the absolute tolerance for a successful solution of the linear system of equations. |
|
|
|
Set this parameter if a matrix free solution procedure should be performed. |
|
|
|
Set the monitor type of the linear solver. |
rte: pseudo time steppingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Scaling parameter of diffusion term. |
|
|
|
Scaling parameter of advection term. |
|
|
|
Determine the pseudo-time step as the product of this scaling and minimum cell size. |
|
|
|
Pseudo-time stepping relative tolerance. |
|
|||
|
rte: pseudo time stepping: time steppingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Defines the start time for the solution of the levelset problem |
|
|
|
Sets the end time for the solution of the levelset problem |
|
|
|
Sets the step size for time stepping. For non-uniform time stepping, this parameter determines the size of the first time step. |
|
|
|
Sets the maximum number of melt_pool steps |
|
|
|
Set an analytical function to determine the time step size. For the prediction of the new time increment, the old time is used. |
rte: pseudo time stepping: linear solverο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter for choosing an iterative linear solver type. |
|
|
|
Set this parameter for choosing a preconditioner type. |
|
|
|
Set the maximum number of iterations for solving the linear system of equations. |
|
|
|
Set the relative tolerance for a successful solution of the linear system of equations. |
|
|
|
Set the absolute tolerance for a successful solution of the linear system of equations. |
|
|
|
Set this parameter if a matrix free solution procedure should be performed. |
|
|
|
Set the monitor type of the linear solver. |
rte: absorptivityο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Sets the absorptivity of the gas phase. |
|
|
|
Sets the absorptivity of the liquid phase. |
|
|
|
Sets the absorptivity of the gas phase. |
π· evaporationο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Choose the formulation how the evaporative mass flux mDot (kg/(m2s)) will be calculated. |
|
|
|
Choose the formulation how the (local) evaporative mass flux will be converted to a DoF vector.will be calculated. When the CutFEM heat transfer operator is used, this input parameter is ignored and the temperature is evaluated at the sharp interface which is equivalent to βsharpβ. |
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|
|
Select the type how the evaporative mass flux should be considered in the level set equation. |
|
|
|
Set if the level set gradient for computing the delta function within the evaporative mass flux source terms should be computed based on an interpolation to the pressure space. This is only implemented for evapor_level_set_source_term_type = rhs. |
evaporation: analyticalο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
For evapor evaporation model == analytical, prescribe a spatially constant mass flux due to evaporation (SI unit in kg/mΒ²s), as a function over time t , e.g. min(2.*t,0.01). |
evaporation: hardt wondraο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Evaporation coefficient for the model by Hardt and Wondra. |
evaporation: pressure awareο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
`` |
Fitting parameters for the evaporative mass flux function with pressure-aware boundary conditions. |
|
|
|
Ambient gas pressure for the pressure-aware model. |
evaporation: evaporative dilation rateο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to consider the evaporative dilation rate in the Navier-Stokes equation. This results in an evaporation-induced jump in the normal velocity component. |
|
|
|
Select how the additional source term due to evaporation in the continuity equation (=evaporative dilation rate) is computed. |
evaporation: evaporative coolingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to consider evaporative cooling in the heat equation |
|
|
|
Enable a linear activation ramp for evaporative cooling between the activation temperature and the boiling temperature. If enabled, the mass flux increases smoothly and linearly within this temperature range. Otherwise, the mass flux is computed directly without applying a ramp. |
|
|
|
Set this parameter to true to account for the enthalpy transported by the vapor mass flux in the heat equation. This is only recommended if the vapor mass flux is not considered in the Navier-Stokes equations. |
|
|
|
Activation temperature for the evaporative cooling. It must be smaller than or equal to the boiling temperature. By default, it will be chosen such that the transition from the linear activation ramp is kink-free. |
|
|
|
Select how the additional source term due to evaporation in the heat equation (evaporative cooling) is computed. |
|
evaporation: evaporative cooling: dirac delta function approximationο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Choose how to smear a parameter over the interface. |
|
|
|
Choose if weights should be computed automatically. |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the heavy phase (level set = 1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the heavy liquid/solid phase (level set = 1). |
evaporation: recoil pressureο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to prescribe the evaporation-induced jump in the pressure field (i.e. recoil pressure), considered as an interfacial force in the momentum balance equation.If βevaporative dilation rateβ is enabled, this pressure jump will be added to the one resulting from the discontinuous normal velocity field. |
|
|
|
Enable a linear activation ramp for recoil pressure between the activation temperature and the boiling temperature. If enabled, the recoil pressure increases smoothly and linearly within this temperature range. Otherwise, the recoil pressure is computed directly without applying a ramp. |
|
|
|
Subtract ambient pressure from the recoil pressure. This can be used to ensure that the recoil pressure is zero at the boiling temperature. |
|
|
|
Ambient gas pressure for the recoil pressure model. |
|
|
|
Pressure coefficient for the recoil pressure model. |
|
|
|
Temperature constant for the recoil pressure model. If this parameter is not set, the value is computed by latent_heat_evaporation * molar_mass / universal_gas_constant; |
|
|
|
Sticking constant. |
|
|
|
Type that determines how the recoil pressure force is computed in the interfacial zone. |
|
|
|
Activation temperature for the recoil pressure. It must be smaller than or equal to the boiling temperature. As default value, the boiling temperature is chosen. |
|
|||
|
|
|
Choose the model to compute the recoil pressure coefficient: phenomenological or hybrid, in case there is also an evaporation-induced velocity jump. |
|
evaporation: recoil pressure: dirac delta function approximationο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Choose how to smear a parameter over the interface. |
|
|
|
Choose if weights should be computed automatically. |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to any phase weighted optionthis parameter controls the (first) weight of the heavy phase (level set = 1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the gas phase (level set = -1). |
|
|
|
If >>> dirac delta function approximation type <<< is set to >>> heaviside_times_heaviside_phase_weighted <<< this parameter controls the second weight of the heavy liquid/solid phase (level set = 1). |
evaporation: recoil pressure: pressure awareο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
`` |
Fitting parameters for the recoil pressure calculation with pressure-aware boundary conditions. |
|
|
|
Ambient gas pressure for the pressure-aware model. |
π· outputο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Sets the base directory for all output. |
|
|
|
Every n timestep that should be written |
|
|
|
Write output output every given time step. If this parameter is set, the output write frequency is deactivated. |
|
|
|
Specify variables that you request to output. |
|
|
|
Set this parameter to true to enable user defined postprocessing. |
|
|||
|
output: paraviewο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to activate paraview output. |
|
|
|
Sets the base name for paraview output files. |
|
|
|
Number of digits for the frame number of the vtu-file. |
|
|
|
Set this parameter to true to output a vtu-file with the boundary id. |
|
|
|
Set this parameter to true to output the subdomain ranks. |
|
|
|
Set to true to output the material id. |
|
|
|
Set this parameter to false to write bi- or trilinear data only. Set this parameter to true to write higher order cell data. Note: higher order cell data can only be written for hexahedron meshes and 2 or 3 dimensions. |
|
|
|
Number of parallel written vtu-files. |
|
|
|
Control number of patches to enable high-order. |
output: particlesο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true to activate particle paraview output. |
|
|
|
Sets the base name for particle output files. |
π· profilingο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true if profiling should be enabled. It will be automaticallyenabled for verbosity level >=1. |
|
|
|
Write profiling output every given time step size. If this parameter is set, the specified parameter for write frequency is overwritten. |
|
|
|
Choose the type of time measure to write profiling information. |
π· application specificο
Parameter |
Type |
Default |
Description |
|---|---|---|---|
|
|
|
Set this parameter to true if you want to consider melting/solidification effects. |
|
|
|
Select the AMR strategy. |