Solids model ============ .. _InputsSolids: Enabling the solids solver and specifying options common to both DEM and PIC models. The following inputs are defined with the prefix ``solids``: +----------------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +============================+==============================================================+==========+==========+ | types | Specified name(s) of the solids types or None to disable | Strings | None | | | the solids solver. The user-defined names are used to | | | | | specify DEM and/or PIC model inputs. | | | +----------------------------+--------------------------------------------------------------+----------+----------+ | species | Specify which species can constitute | Strings | None | | | the solids phase. Defined species must | | | | | be a subset of the species defined in ``species.solve`` | | | +----------------------------+--------------------------------------------------------------+----------+----------+ Molecular Weight ---------------- The following inputs are defined with the prefix ``solids``: +----------------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +============================+==============================================================+==========+==========+ | molecular_weight | Value of constant solid molecular weight. | Real | 0 | +----------------------------+--------------------------------------------------------------+----------+----------+ Specific Heat ------------- The following inputs are defined with the prefix ``solids.specific_heat``: +----------------------------+----------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +============================+================================================================+==========+==========+ | model | Specify which specific heat model to use for solids. | String | constant | | | | | | | | A specific heat model is required if advecting enthalpy. | | | | | Furthermore, the model must be ``mixture`` if species | | | | | equations are solved. | | | | | | | | | | Options: | | | | | | | | | | * ``constant`` - constant solid specific heat | | | | | * ``NASA7-poly`` the solid specific heat is defined by | | | | | NASA-7 polynomials. | | | | | | | | | | NASA7 polynomial format: | | | | | | | | | | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | | | | | | | | Additional information on NASA-7 polynomials is provided | | | | | in :ref:`species specific heats`. | | | | | | | | | | * ``mixture`` - a *mixture* specific heat is computed | | | | | from :ref:`species specific heats` | | | | | and local species mass fractions. | | | | | | | | | | :math:`c_{p,\mathrm{mixture}} = \sum_n X_n c_{p,n}` | | | | | | | | +----------------------------+----------------------------------------------------------------+----------+----------+ | constant | Constant solids specific heat value. | Real | 0 | | | A value is required for ``constant`` specific heat model. | | | +----------------------------+----------------------------------------------------------------+----------+----------+ | ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | | | temperatures (``Tsplit``). More information is provided in | | | | | :ref:`species specific heats`. | | | | | | | | +----------------------------+----------------------------------------------------------------+----------+----------+ | NASA7.a[i] | Specific heat polynomial coefficients. Polynomial | Reals | None | | | coefficients are required if the fluid specific heat | | | | | model is ``NASA7-poly``. | | | | | | | | | | Additional information on NASA-7 polynomials is provided | | | | | in :ref:`species specific heats`. | | | +----------------------------+----------------------------------------------------------------+----------+----------+ | NASA7.Tsplit | Defines the transition temperature between NASA-7 | Reals | 1000 | | | polynomials. | | | | | | | | | | If there are ``N`` sets of coefficients for ``N`` temperature | | | | | ranges, then ``N-1`` transition temperatures must be specified.| | | | | | | | | | In the standard case of two temperature ranges, a | | | | | transition temperature of 1000K is assumed if ``Tsplit`` is | | | | | not provided. | | | +----------------------------+----------------------------------------------------------------+----------+----------+ The following inputs are defined with the prefix ``solids``: +----------------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +============================+==============================================================+==========+==========+ | reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | +----------------------------+--------------------------------------------------------------+----------+----------+ Thermal conductivity -------------------- .. |Sutherland_k_Eq| replace:: :math:`\kappa(T) = \kappa_{ref}\left(\frac{T}{T_{ref}}\right)^{3/2}\frac{T_{ref} + S}{T+S}` .. |PowerLaw_k_Eq| replace:: :math:`\kappa(T) = \left( \frac{T}{T_{ref}} \right)^n` The following inputs are defined using the prefix ``[solid_name].thermal_conductivity``. +-----------------------+----------------------------------------------------------------+--------+----------+ | | Description | Type | Default | +=======================+================================================================+========+==========+ | model | Specify the thermal conductivity model. Options: | String | None | | | | | | | | * ``constant`` - constant thermal conductivity | | | | | * ``Sutherland`` :cite:p:`suth1893` - | | | | | | | | | | |Sutherland_k_Eq| | | | | | | | | | | * ``power_law`` : Power law formulation | | | | | | | | | | |PowerLaw_k_Eq| | | | | | | | | | | A thermal conductivity model is required if advecting | | | | | enthalpy. | | | +-----------------------+----------------------------------------------------------------+--------+----------+ | constant | Value of constant thermal conductivity. | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ | Sutherland.T_ref | Sutherland model reference temperature. | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ | Sutherland.k_ref | Sutherland model reference conductivity at T_ref. | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ | Sutherland.S | Sutherland model temperature. | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ | power_law.T_ref | Power law model reference temperature. | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ | power_law.exponent | Power law model exponent :math:`n` | Real | None | +-----------------------+----------------------------------------------------------------+--------+----------+ The following inputs are defined with the prefix ``solids``: +-----------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +=======================+==============================================================+==========+==========+ | min_conduction_dist | The surface roughness separating two touching particles. | Real | 2.e-8 | | | This value is used to remove the singularity that the | | | | | particle-fluid-particle conduction model develops at | | | | | the contact interface. | | | +-----------------------+--------------------------------------------------------------+----------+----------+ | flpc | The fluid lens proportionality constant (flpc) is used to | Real | 0.4 | | | calculate the radius of the fluid lens surrounding the | | | | | particle for particle-fluid-particle conduction. | | | +-----------------------+--------------------------------------------------------------+----------+----------+ Density ------- The following inputs are defined with the prefix ``solids``: +-----------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +=======================+==============================================================+==========+==========+ | density | Specify which density model to use for solids. | String | None | | | | | | | | Options (case-insensitive): | | | | | | | | | | * ``None`` - solids density specified through initial and | | | | | boundary conditions; | | | | | * ``Mixture`` - solids density computed as a weighted | | | | | average of the solid phase species densities. It requires | | | | | constant :ref:`species densities` to | | | | | be defined by the user. | | | | | | | | +-----------------------+--------------------------------------------------------------+----------+----------+ Newton solver ------------- The following inputs are defined with the prefix ``solids``: +------------------------------+--------------------------------------------------------------+----------+----------+ | | Description | Type | Default | +==============================+==============================================================+==========+==========+ | newton_solver.absolute_tol | Define absolute tolerance for damped Newton solver. | Real | 1.e-6 | +------------------------------+--------------------------------------------------------------+----------+----------+ | newton_solver.relative_tol | Define relative tolerance for damped Newton solver. | Real | 1.e-6 | +------------------------------+--------------------------------------------------------------+----------+----------+ | newton_solver.max_iterations | Define max number of iterations for damped Newton solver. | Int | 100 | +------------------------------+--------------------------------------------------------------+----------+----------+ DEM model settings ------------------ .. |RollingFricModelA| replace:: :math:`\boldsymbol{\tau}_{ij}^{(r)} = -\mu_{r}\hat{r}_{ij} \left| f_{ij}^{(n)} \right| \frac{\omega_{ij}}{\left|\omega_{ij}\right|}` .. |RollingFricModelB| replace:: :math:`\boldsymbol{\tau}_{ij}^{(r)} = -\mu_{r}\hat{r}_{ij} \left| f_{ij}^{(n)} \right| \boldsymbol{u}_{ij}^{\omega}` Enabling the DEM solver and specifying model options. The following keys must be defined with the prefix ``dem``: +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | | Description | Type | Default | +===============================+=========================================================================+==========+===========+ | solve | Specified name(s) of the DEM types or None to disable the DEM solver. | Strings | None | | | The user-defined names are used to specify DEM model inputs. | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | friction_coeff.pp | Friction coefficient :: particle to particle collisions [required] | Real | 0 | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | friction_coeff.pw | Friction coefficient :: particle to wall collisions [required] | Real | 0 | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | spring_const.pp | Normal spring constant :: particle to particle collisions [required] | Real | 0 | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | spring_const.pw | Normal spring constant :: particle to wall collisions [required] | Real | 0 | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | spring_tang_fac.pp | Tangential-to-normal spring constant factor :: particle to particle | Real | 0.2857 | | | collisions | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | spring_tang_fac.pw | Tangential-to-normal spring constant factor :: particle to wall | Real | 0.2857 | | | collisions | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | damping_tang_fac.pp | Factor relating the tangential damping coefficient to the normal | Real | 0.5 | | | damping coefficient :: particle to particle collisions | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | damping_tang_fac.pw | Factor relating the tangential damping coefficient to the normal | Real | 0.5 | | | damping coefficient :: particle to wall collisions | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | implicit_drag | Apply fluid-particle drag force by | Int | 0 | | | :ref:`implicit velocity update`. | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | tan_history | Include tangential history force in the collision model. | Bool | false | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | tan_history.max_contacts | Maximum number of contacts per particle tracked at any point when | Int | 10 | | | tangential history is enabled. | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | rolling_friction | Rolling friction model to use when tangential history is enabled. | String | None | | | :cite:p:`zhou1999,ai2011,wensrich2012,blais2019` | | | | | | | | | | Options: | | | | | | | | | | * ``None`` | | | | | | | | | | * ``ModelA`` | | | | | | | | | | |RollingFricModelA| | | | | | | | | | | * ``ModelB`` | | | | | | | | | | |RollingFricModelB| | | | | | | | | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ | rolling_friction.coefficient | Rolling friction coefficient when using a rolling friction model | Real | 0 | +-------------------------------+-------------------------------------------------------------------------+----------+-----------+ The following inputs use the DEM type names specified using the `dem.solve` input to define restitution coefficients and are proceeded with the prefix ``dem.restitution_coeff``. These must be defined for all DEM solid-solid and solid-wall combinations. +--------------------------+-------------------------------------------------------------------------+----------+-----------+ | | Description | Type | Default | +==========================+=========================================================================+==========+===========+ | [dem_solid].[dem_solid2] | Specifies the restitution coefficient between DEM solids. | Real | 0 | | | Order is not important and this could be defined as | | | | | ``[dem_solid2].[dem_solid]``. If both are defined they must be equal. | | | +--------------------------+-------------------------------------------------------------------------+----------+-----------+ | [dem_solid].wall | Specifies the restitution coefficient between solid and the wall. | Real | 0 | | | Order is not important and this could be defined as | | | | | ``wall.[dem_solid]``. If both are defined they must be equal. | | | +--------------------------+-------------------------------------------------------------------------+----------+-----------+ Below is an example for specifying the inputs for two DEM solids. .. code-block:: none dem.solve = sand char dem.friction_coeff.pp = 0.25 dem.friction_coeff.pw = 0.15 dem.spring_const.pp = 100.0 dem.spring_const.pw = 100.0 dem.spring_tang_fac.pp = 0.2857 dem.spring_tang_fac.pw = 0.2857 dem.damping_tang_fac.pp = 0.5 dem.damping_tang_fac.pw = 0.5 dem.restitution_coeff.sand.sand = 0.85 dem.restitution_coeff.sand.char = 0.88 dem.restitution_coeff.char.char = 0.90 dem.restitution_coeff.sand.wall = 0.85 dem.restitution_coeff.char.wall = 0.89 .. _InputsDEMImplicitDrag: Implicit drag ~~~~~~~~~~~~~ By default, the fluid-particle drag force is applied explicitly to particles, .. math:: u_p^{n^{\prime}+1} = u_p^{n^\prime} + dt_{\mathrm{DEM}}\left(\beta^n (u_f^{n+1} - u_p^{n}) + \mathcal{F}_p^{n^\prime}\right)/m_p where :math:`u_p` is the `p`-th particle velocity at the indicated time level, :math:`\beta^n` is the drag coefficient, :math:`u_f^{n+1}` is the fluid velocity, :math:`\mathcal{F}_p` are all non-drag forces (e.g., gravitational, buoyancy, collision etc.), and :math:`m_p` is particle mass. Terms exist at different time levels due in part to how the fluid and particle models are coupled and partly because particles typically sub-step in time. Specifically, a time step advances the fluid from :math:`t^n` to :math:`t^{n+1}`, then particles take multiple smaller time steps to traverse the same total time. In the above equation, prime markers differentiate quantities updated at each DEM sub-step. The updated particle velocity, :math:`u_p^{n^{\prime}+1}`, is substituted into the drag expression when `implicit drag` is enabled. .. math:: u_p^{n^{\prime}+1} = \frac{u_p^{n^\prime} + dt_{\mathrm{DEM}}\left( \beta^n u_f^{n+1} + \mathcal{F}_p^{n^\prime} \right)/m_p}{1 + dt_{\mathrm{DEM}} \beta^n / m_p} By using the updated velocity, the drag force computed for the fluid and particles is inconsistent, and interphase momentum is not conserved. However, implicit drag may be needed for numerical stability when the particle density is much less than the fluid (such as bubbles). PIC model settings ------------------ The following inputs define the PIC solids model and its associated frictional stress parameters: .. math:: \tau_p = \frac{P_s \varepsilon_p^{\beta}}{\max \left[ (\varepsilon_{cp} - \varepsilon_p), \epsilon (1-\varepsilon_p)\right]} The following inputs are defined with the prefix ``pic``: +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | | Description | Type | Default | +====================================+=========================================================================+==========+===========+ | solve | Specified name(s) of the PIC types or None to disable the PIC solver. | Strings | None | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | close_pack | Solids volume fraction at close pack, :math:`\varepsilon_{cp}`. | Real | 0.55 | | | Below this threshold, the solid stress tensor has minimal impact. | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | pressure_coefficient | Pressure linear scale factor in frictional stress model, :math:`P_s`. | Real | 100.0 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | beta | Volume fraction exponential scale factor in frictional stress model, | Real | 3.0 | | | :math:`\beta`. | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | small_number | Non-singularity term in frictional stress model, :math:`\epsilon`. | Real | 1.0e-7 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | damping_factor | An empirical dampening factor for the frictional stress model. | Real | 0.4 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | damping_factor_wall_normal | Normal coefficient of restitution for parcel-wall collisions. | Real | 0.3 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | damping_factor_wall_tangent | Tangential coefficient of restitution for parcel-wall collisions. | Real | 0.99 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | parcels_per_cell_at_pack | Specifies the number of particles that can occupy a regular | Real | 32.0 | | | computational cell near close packing. | | | | | | | | | | Larger values reduce sensitivity to small changes in parcel positions, | | | | | which helps stabilize the computed solids volume fraction and, by | | | | | extension, the frictional stress model. However, increasing this value | | | | | also results in simulations tracking a greater number of parcels. | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | verbose | Verbosity level for PIC-related output. | Int | 0 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ The solids volume fraction is required to evaluate the frictional stress model when advancing parcels from :math:`t^n` to :math:`t^{n+1}`. Ideally, this evaluation would use either the volume fraction at the midpoint of the time step, :math:`\varepsilon_p^{n+1/2}`, or at the end of the step, :math:`\varepsilon_p^{n+1}`. However, since neither of these values is available a priori, an iterative approach is used. An initial approximation of the volume fraction is used to advance parcels to a tentative position at :math:`t^{n+1}`. These new positions are then used to compute an updated volume fraction field, which is subsequently used in the next iteration to recompute the tentative parcel positions. The following keywords control the iterative PIC algorithm and are defined with the prefix ``pic``: +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | | Description | Type | Default | +====================================+=========================================================================+==========+===========+ | initial_step_type | Specifies the initial solids volume fraction used in PIC iterations to | String | "nth_eps" | | | approximate the volume fraction at the end of the time step, | | | | | :math:`t^{n+1}` | | | | | | | | | | Options (case-insensitive): | | | | | | | | | | * ``nth_eps`` - use the volume fraction from the previous time step, | | | | | :math:`t^{n}` | | | | | * ``zero_eps`` - set the initial volume fraction to zero | | | | | * ``taylor_approx`` - use a Taylor series expansion about parcel | | | | | positions to approximate the volume fraction at :math:`t^{n+1}` | | | | | | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | max_iter | Maximum number of iterations used to approximate the end state | Int | 3 | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | advance_vel_p | Specifies the weighting factor used to compute the averaged particle | Real | 0.5 | | | velocity field during parcel advancement. The field-averaged velocity | | | | | is used when applying forces to parcels. | | | | | | | | | | This parameter takes values between 0 and 1, blending the particle | | | | | velocity at the beginning and end of the time step. A value of ``0.5`` | | | | | (default) yields a time-centered average. | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | velocity_reference_frame | The velocity reference frame scales the field-averaged velocity | Real | 0.5 | | | used to compute the slip velocity between a parcel and the bulk when | | | | | when the parcels is moving towards an area of higher solids | | | | | concentration. This parameter takes values between 0 and 1. | | | +------------------------------------+-------------------------------------------------------------------------+----------+-----------+