From b0b13702db9565f556c5fec6b4d2325f68991c05 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:24:41 -0600 Subject: [PATCH 01/16] clarify --- docs/source_docs/user_guide/inputs/porous_media.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/porous_media.rst b/docs/source_docs/user_guide/inputs/porous_media.rst index 0de4984..6ab3c1f 100644 --- a/docs/source_docs/user_guide/inputs/porous_media.rst +++ b/docs/source_docs/user_guide/inputs/porous_media.rst @@ -3,7 +3,7 @@ Porous media definitions ======================== -Homogeneous porous media are modeled by the addition of a momentum +Homogeneous porous media are modeled by the introduction of a momentum source term, :math:`\boldsymbol{S}_{pm}`, to the fluid momentum equation. .. math:: -- GitLab From 852e0cec603e21f7b502d82f97d9bb44d8aae890 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:24:58 -0600 Subject: [PATCH 02/16] whitespace --- .../references/fluids/FluidEquations.rst | 61 +++++++++---------- 1 file changed, 29 insertions(+), 32 deletions(-) diff --git a/docs/source_docs/references/fluids/FluidEquations.rst b/docs/source_docs/references/fluids/FluidEquations.rst index e74da21..456cbb8 100644 --- a/docs/source_docs/references/fluids/FluidEquations.rst +++ b/docs/source_docs/references/fluids/FluidEquations.rst @@ -22,18 +22,18 @@ Conservation of fluid mass: .. math:: \frac{\partial (\varepsilon_g \rho_g)}{\partial t} + \nabla \cdot (\varepsilon_g \rho_g U_g) = 0 -Unlike two-fluid modeling, :math:`\varepsilon_g = 1 - \varepsilon_s`, is not a solution -variable. Rather, :math:`\varepsilon_s` is evaluated from the particle field through -volume filtering, +Unlike two-fluid modeling, :math:`\varepsilon_g = 1 - \varepsilon_s`, is not a solution +variable. Rather, :math:`\varepsilon_s` is evaluated from the particle field through +volume filtering, -.. math:: - (1 - \varepsilon_g) A (\boldsymbol{x},t) \approx - \sum_{i=1}^{N_p} A_i(\boldsymbol{X}_i,t) {\mathcal G} - (\left| \boldsymbol{x} - \boldsymbol{X}_i \right|) {\mathcal V}_i +.. math:: + (1 - \varepsilon_g) A (\boldsymbol{x},t) \approx + \sum_{i=1}^{N_p} A_i(\boldsymbol{X}_i,t) {\mathcal G} + (\left| \boldsymbol{x} - \boldsymbol{X}_i \right|) {\mathcal V}_i -where :math:`A_i` is a general particle property, :math:`{\mathcal V}_i` is the particle -volume and :math:`\mathcal G` is a transfer kernel with compact support--here linear hat. -Setting :math:`A_i = 1` gives the local void fraction. +where :math:`A_i` is a general particle property, :math:`{\mathcal V}_i` is the particle +volume and :math:`\mathcal G` is a transfer kernel with compact support--here linear hat. +Setting :math:`A_i = 1` gives the local void fraction. Assuming the fluid phase is incompressible, :math:`\frac{D \rho_g}{Dt} = 0`, the conservation of fluid mass is equivalent to the conservation of fluid volume: @@ -42,42 +42,39 @@ Assuming the fluid phase is incompressible, :math:`\frac{D \rho_g}{Dt} = 0`, the The conservation of fluid momentum is: -.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} +.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g = \nabla \cdot \tau + M_{sg} + \varepsilon_g \rho_g g -where :math:`M_{sg} = - M_{gs}` is the generalized interfacial momentum transfer from -the solid particles to the fluid-phase. Like :math:`\varepsilon_s`, -:math:`M_{gs}` is determined from the L-E transfer kernel by setting :math:`A_i = F_{gi}`, -where :math:`F_{gi}` is the force due to the fluid-phase on the ith particle. Following -MFiX classic (and many other CFD-DEM codes designed for high density ratio gas-solids -flows), only buoyancy (pressure gradient) and steady drag are considered: +where :math:`M_{sg} = - M_{gs}` is the generalized interfacial momentum transfer from +the solid particles to the fluid-phase. Like :math:`\varepsilon_s`, +:math:`M_{gs}` is determined from the L-E transfer kernel by setting :math:`A_i = F_{gi}`, +where :math:`F_{gi}` is the force due to the fluid-phase on the ith particle. Following +MFiX classic (and many other CFD-DEM codes designed for high density ratio gas-solids +flows), only buoyancy (pressure gradient) and steady drag are considered: .. math:: - F_{gi} = - \mathcal{V}_i \nabla p_g + F_{gi} = - \mathcal{V}_i \nabla p_g - \frac{1}{2} C_D \rho_g \boldsymbol{V}_{ig} \left|\boldsymbol{V}_{ig}\right| A_i^{(proj)} -where :math:`\boldsymbol{V}_{ig} = \boldsymbol{V}_i - \boldsymbol{U}_g ( \boldsymbol{X}_i )` -is the velocity of ith particle relative to the fluid-phase (at the particle position -:math:`\boldsymbol{X}_i`). :math:`F_{gi}` is closed by the specification of a drag -coefficient, :math:`C_D`. Currently, MFIX-Exa includes Wen-Yu, Gidaspow and BVK2 drag laws. +where :math:`\boldsymbol{V}_{ig} = \boldsymbol{V}_i - \boldsymbol{U}_g ( \boldsymbol{X}_i )` +is the velocity of ith particle relative to the fluid-phase (at the particle position +:math:`\boldsymbol{X}_i`). :math:`F_{gi}` is closed by the specification of a drag +coefficient, :math:`C_D`. Currently, MFIX-Exa includes Wen-Yu, Gidaspow, BVK2, and Syamlal-O'Brien drag laws. -.. wdf todo +.. wdf todo gidaspow form: C_D^{Gidaspow} = \chi C_D^{(Wen-Yu)} + \left(1 - \chi \right) C_D^{(Ergun)} \chi = \frac{\arctan 150 \left( \varepsilon_g - 0.8 \right)}{\pi} + \frac{1}{2} wen-yu: C_D^{(Wen-Yu)} = \max \left[\frac{24}{Re_i}\left(1 + 0.15 Re_i^{0.687}\right),\ 0.44 \right] \left( 1 - \varepsilon_g \right)^{-1.65} - ergun: C_D^{(Ergun)} = \frac{200 \left(1 - \varepsilon_g\right)}{Re_i} + \frac{7}{3} + ergun: C_D^{(Ergun)} = \frac{200 \left(1 - \varepsilon_g\right)}{Re_i} + \frac{7}{3} BVK2: C_D^{BVK2} = ... particle reynolds number: Re_i = \frac{\rho_g \left( 1 - \varepsilon_g \right) d_i \left| \boldsymbol{V}_{ig} \right| }{\mu_g} -In chemical engineering literature, it is common to lump all drag-related terms of -:math:`M_{gs}` into :math:`\beta`. With this simplification and some re-arrangement, -the fluid momentum takes the more convenient form: +In chemical engineering literature, it is common to lump all drag-related terms of +:math:`M_{gs}` into :math:`\beta`. With this simplification and some re-arrangement, +the fluid momentum takes the more convenient form: -.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} - + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g +.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} + + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g = \nabla \cdot \tau + \sum_p \beta_p (V_p - U_g) + \varepsilon_g \rho_g g - - - -- GitLab From 88a5bd5e571b69417d2846f6fb495ea41da44c27 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:25:08 -0600 Subject: [PATCH 03/16] spell out Syamlal --- docs/source_docs/user_guide/inputs/model_options.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 92ae7eb..7983264 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -318,7 +318,7 @@ The following input is defined using the prefix ``mfix.drag``: | | * ``WenYu`` :cite:p:`Wen66` | | | | | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | | | * ``BVK2`` :cite:p:`Bee07,Bee071,tenneti11,tang15` | | | -| | * ``SyamOBrien`` :cite:p:`syam88` | | | +| | * ``SyamOBrien`` -- Syamlal-O'Brien :cite:p:`syam88` | | | | | * ``UserDrag`` | | | | | | | | +--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -- GitLab From f38abdbd45296335f17cff3e008db978eac5fc06 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:07:45 -0600 Subject: [PATCH 04/16] move mfix.solve_species to species section --- docs/source_docs/user_guide/inputs/model_options.rst | 2 -- 1 file changed, 2 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 7983264..658b9c3 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -13,8 +13,6 @@ The following inputs are defined using the prefix ``mfix``: +------------------------+----------------------------------------------------------------------------+----------+---------------------+ | advect_enthalpy | Enable time evolution of fluid temperature and enthalpy. | Int | 0 | +------------------------+----------------------------------------------------------------------------+----------+---------------------+ -| solve_species | Enable the time evolution of fluid species mass fractions. | Int | 0 | -+------------------------+----------------------------------------------------------------------------+----------+---------------------+ | constraint | Select low Mach number constraint. | String | IncompressibleFluid | | | | | | | | Options (case-insensitive): | | | -- GitLab From aae29571634f74f57e14fc3e12f7cdb6e051a1a4 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:11:56 -0600 Subject: [PATCH 05/16] update docs move mfix.solve_species to species section rename ignore_discontinuities to specific_heat.ignore_discontinuities --- .../user_guide/inputs/fluid_model.rst | 420 +++++++++--------- .../user_guide/inputs/initial_conditions.rst | 2 +- .../user_guide/inputs/model_options.rst | 79 ++-- .../user_guide/inputs/solids_model.rst | 172 ++++--- .../source_docs/user_guide/inputs/species.rst | 286 ++++++------ duplicate-keys.txt | 7 - nonexistent-keys.txt | 1 - 7 files changed, 516 insertions(+), 451 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 198894b..2b6d787 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -27,210 +27,226 @@ a single fluid; therefore, it is common to name the fluid ``fluid``. This is ill .. |Sato_Eq| replace:: :math:`\mu_{pit} = C d_s \rho \left|\boldsymbol{u} - \boldsymbol{u_s}\right|` .. |eff_visc| replace:: :math:`\mu_{eff} = \mu_{mol} + \mu_{eddy} + \mu_{susp} + \mu_{pit}` - .. |mix_Eq| replace:: :math:`\mu_{mix} = \sum_{\alpha=1}^{N} \frac{X_{\alpha} \mu_{\alpha}}{\sum_{\beta} X_{\beta} \phi_{\alpha \beta}}` -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| | Description | Type | Default | -+============================================+============================================================+========+==========+ -| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - no viscosity model. | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` - | | | -| | | | | -| | |Sutherland_Eq| | | | -| | | | | -| | * ``Reid`` :cite:p:`reid87` - | | | -| | | | | -| | |Reid_4parm_Eq| | | | -| | | | | -| | * ``mixture`` - a mixture viscosity is computed from | | | -| | :ref:`species viscosities` | | | -| | and local species mass fractions :cite:p:`bird06` | | | -| | | | | -| | |mix_Eq| | | | -| | | | | -| | A viscosity model is required if the fluid solver is | | | -| | enabled. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy | Specify eddy viscosity model. | | | -| | | | | -| | Options: | | | -| | |String | None | -| | * ``None`` - No eddy viscosity model | | | -| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | -| | * ``WALE`` :cite:p:`ducros98` | | | -| | * ``usr`` | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | -| | between 0.1 and 0.2. | | | -| | | | | -| | A value is required when using the ``Smagorinsky-Lilly`` | | | -| | eddy viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension | Specify suspension viscosity model of the form | String | None | -| | | | | -| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - No eddy suspension model | | | -| | | | | -| | :math:`\mu^* = 1` | | | -| | | | | -| | * ``Einstein`` :cite:p:`einstein11` | | | -| | | | | -| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | -| | | | | -| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | -| | | | | -| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | -| | | | | -| | * ``ChengLaw`` :cite:p:`cheng03` | | | -| | | | | -| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | -| | | | | -| | * ``Sato`` :cite:p:`sato81` | | | -| | | | | -| | |Sato_Eq| | | | -| | | | | -| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | -| | | | | -| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Brinkman`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``ChengLaw`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | -| | molecular viscosity | | | -| | | | | -| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | -| | |eff_visc| | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| species | Specify the species that constitute the fluid. | Strings| None | -| | | | | -| | All listed species must be properly defined. See the | | | -| | :ref:`species definition documentation` | | | -| | for additional details. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| molecular_weight | Constant fluid molecular weight. | Real | 0 | -| | | | | -| | A molecular weight should only be defined when using one | | | -| | of the ``IdealGas`` constraints and the species equations | | | -| | are not solved. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat | Specify which specific heat model to use for the fluid. | String | None | -| | Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant specific heat | | | -| | | | | -| | * ``NASA7-poly`` the fluid 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}` | | | -| | | | | -| | A specific heat model is required if advecting enthalpy. | | | -| | Furthermore, the model must be ``mixture`` if species | | | -| | equations are solved. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.constant | Constant fluid specific heat | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.NASA7.a[i] | Specific heat polynomial coefficients. Polynomial | Real | 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`. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | -| | fluid. Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant thermal | | | -| | conductivity model | | | -| | | | | -| | A thermal conductivity model is required if advecting | | | -| | enthalpy. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | -| | A value is required when using the | | | -| | ``IdealGasClosedSystem`` constraint. If using the | | | -| | ``IdealGasOpenSystem`` constraint, then the value must | | | -| | match all specified outflow pressures. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| reference_temperature | Reference temperature used to compute specific | Real | 298. | -| | enthalpy | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| | Description | Type | Default | ++============================================+================================================================+========+==========+ +| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model. | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` - | | | +| | | | | +| | |Sutherland_Eq| | | | +| | | | | +| | * ``Reid`` :cite:p:`reid87` - | | | +| | | | | +| | |Reid_4parm_Eq| | | | +| | | | | +| | * ``mixture`` - a mixture viscosity is computed from | | | +| | :ref:`species viscosities` | | | +| | and local species mass fractions :cite:p:`bird06` | | | +| | | | | +| | |mix_Eq| | | | +| | | | | +| | A viscosity model is required if the fluid solver is | | | +| | enabled. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy | Specify eddy viscosity model. | | | +| | | | | +| | Options: | | | +| | |String | None | +| | * ``None`` - No eddy viscosity model | | | +| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | +| | * ``WALE`` :cite:p:`ducros98` | | | +| | * ``usr`` | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | +| | between 0.1 and 0.2. | | | +| | | | | +| | A value is required when using the ``Smagorinsky-Lilly`` | | | +| | eddy viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension | Specify suspension viscosity model of the form | String | None | +| | | | | +| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - No eddy suspension model | | | +| | | | | +| | :math:`\mu^* = 1` | | | +| | | | | +| | * ``Einstein`` :cite:p:`einstein11` | | | +| | | | | +| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | +| | | | | +| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | +| | | | | +| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | +| | | | | +| | * ``ChengLaw`` :cite:p:`cheng03` | | | +| | | | | +| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | +| | | | | +| | * ``Sato`` :cite:p:`sato81` | | | +| | | | | +| | |Sato_Eq| | | | +| | | | | +| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | +| | | | | +| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Brinkman`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``ChengLaw`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | +| | molecular viscosity | | | +| | | | | +| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | +| | |eff_visc| | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| species | Specify the species that constitute the fluid. | Strings| None | +| | | | | +| | All listed species must be properly defined. See the | | | +| | :ref:`species definition documentation` | | | +| | for additional details. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| molecular_weight | Constant fluid molecular weight. | Real | 0 | +| | | | | +| | A molecular weight should only be defined when using one | | | +| | of the ``IdealGas`` constraints and the species equations | | | +| | are not solved. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat | Specify which specific heat model to use for the fluid. | String | None | +| | Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant specific heat | | | +| | | | | +| | * ``NASA7-poly`` the fluid 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.constant | Constant fluid specific heat | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provided | | | +| | in :ref:`species specific heats`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.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`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single | | | +| | polynomial defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat | | | +| | and enthalpy. | | | +| | | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | +| | fluid. Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant thermal | | | +| | conductivity model | | | +| | | | | +| | A thermal conductivity model is required if advecting | | | +| | enthalpy. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | +| | A value is required when using the | | | +| | ``IdealGasClosedSystem`` constraint. If using the | | | +| | ``IdealGasOpenSystem`` constraint, then the value must | | | +| | match all specified outflow pressures. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| reference_temperature | Reference temperature used to compute specific | Real | 298. | +| | enthalpy | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ Example inputs diff --git a/docs/source_docs/user_guide/inputs/initial_conditions.rst b/docs/source_docs/user_guide/inputs/initial_conditions.rst index 579419c..5c1397c 100644 --- a/docs/source_docs/user_guide/inputs/initial_conditions.rst +++ b/docs/source_docs/user_guide/inputs/initial_conditions.rst @@ -20,7 +20,7 @@ The following inputs are defined using the prefix ``ic``: | | * ``Volume`` - sort by volume of each IC region | | | | | * ``Priority`` - sort by specified priority value | | | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ -| [region_name].priority| Priority value for IC regions if ``ranking_type` == ``Priority`` | Int | Max | +| [region_name].priority| Priority value for IC regions if ``ranking_type`` == ``Priority`` | Int | Max | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 658b9c3..f5c0338 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -308,45 +308,46 @@ The following input is defined using the prefix ``mfix.drag``: .. |VM_Nijssen_Eq| replace:: :math:`\begin{split}C_p,vm =& C_0,vm \\ &+ \sum_{k=1}^3 \left\left[a_{k,0} + a_{k,1} \ln\left(\frac{\rho_p}{\rho_f}\right)\right] \varepsilon_p^k \right\end{split}` -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| | Description | Type | Default | -+==========================+===============================================================================+=============+===========+ -| model | Fluid-particle drag model. Options: | String | None | -| | | | | -| | * ``WenYu`` :cite:p:`Wen66` | | | -| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | -| | * ``BVK2`` :cite:p:`Bee07,Bee071,tenneti11,tang15` | | | -| | * ``SyamOBrien`` -- Syamlal-O'Brien :cite:p:`syam88` | | | -| | * ``UserDrag`` | | | -| | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | -| model.SyamOBrien.d1 | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | -| | fluid-particle drag force. The force is applied to the fluid by | | | -| | multiplying the shear stress by fluid volume fraction. | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | -| | The virtual mass force is not included by default. | | | -| | | | | -| | |VirtualMass_Eq| | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - virtual mass force calculation is skipped | | | -| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | -| | * ``Constant`` - specified constant coefficient | | | -| | * ``Zuber`` - :cite:p:`zuber64` | | | -| | | | | -| | |VM_Zuber_Eq| | | | -| | | | | -| | * ``Nijssen`` - :cite:p:`nijssen20` | | | -| | | | | -| | |VM_Nijssen_Eq| | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| | Description | Type | Default | ++==========================+=========================================================================================+=============+===========+ +| model | Fluid-particle drag model. Options: | String | None | +| | | | | +| | * ``WenYu``- Wen-Yu :cite:p:`Wen66` | | | +| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | +| | * ``BVK2`` - Beetstra, van der Hoef, and Kuipers :cite:p:`Bee07,Bee071,tenneti11,tang15`| | | +| | * ``SyamOBrien`` - Syamlal-O'Brien :cite:p:`syam88` | | | +| | * ``UserDrag`` - User-supplied drag model | | | +| | | | | +| | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | +| model.SyamOBrien.d1 | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | +| | fluid-particle drag force. The force is applied to the fluid by | | | +| | multiplying the shear stress by fluid volume fraction. | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | +| | The virtual mass force is not included by default. | | | +| | | | | +| | |VirtualMass_Eq| | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - virtual mass force calculation is skipped | | | +| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | +| | * ``Constant`` - specified constant coefficient | | | +| | * ``Zuber`` - :cite:p:`zuber64` | | | +| | | | | +| | |VM_Zuber_Eq| | | | +| | | | | +| | * ``Nijssen`` - :cite:p:`nijssen20` | | | +| | | | | +| | |VM_Nijssen_Eq| | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ .. note:: diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index 1436464..9ad6278 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -9,71 +9,113 @@ Solids settings 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| molecular_weight | Value of constant solid molecular weight. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| density | Specify which density model to use for solids. | String | None | -| | | | | -| | Options (case-insensive): | | | -| | | | | -| | * ``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. | | | -| | | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat | Specify which specific heat model to use for solids. | String | constant | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant specific heat | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat.constant | Constant solids specific heat value. | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| species | Specify which species can constitute | Strings | None | -| | the solids phase. Defined species must | | | -| | be a subset of the species defined in ``species.solve`` | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | -| | Options: | | | -| | | | | -| | * ``constant`` - constant thermal conductivity | | | -| | | | | -| | A thermal conductivity model is not required for solids. | | | -| | Defining a thermal conductivity model enables conductive | | | -| | heat transfer for DEM. However, it has no effect for PIC | | | -| | solids heat transfer. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| | 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| molecular_weight | Value of constant solid molecular weight. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| density | Specify which density model to use for solids. | String | None | +| | | | | +| | Options (case-insensive): | | | +| | | | | +| | * ``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. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat | Specify which specific heat model to use for solids. | String | constant | +| | | | | +| | 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.constant | Constant solids specific heat value. | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provied in | | | +| | :ref:`species specific heats`. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.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`. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 | Reals | 1000 | +| | polynomials. | | | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial | | | +| | defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat and | | | +| | enthalpy. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| species | Specify which species can constitute | Strings | None | +| | the solids phase. Defined species must | | | +| | be a subset of the species defined in ``species.solve`` | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | +| | Options: | | | +| | | | | +| | * ``constant`` - constant thermal conductivity | | | +| | | | | +| | A thermal conductivity model is not required for solids. | | | +| | Defining a thermal conductivity model enables conductive | | | +| | heat transfer for DEM. However, it has no effect for PIC | | | +| | solids heat transfer. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ Below is an example for specifying the solids solver model options. diff --git a/docs/source_docs/user_guide/inputs/species.rst b/docs/source_docs/user_guide/inputs/species.rst index 17037ad..f02aa4e 100644 --- a/docs/source_docs/user_guide/inputs/species.rst +++ b/docs/source_docs/user_guide/inputs/species.rst @@ -1,149 +1,163 @@ .. _InputsSpeciesDefs: + +Species solver +============== + +Enabling the species solver +--------------------------- + +The following input is defined using the prefix ``mfix``: + ++-------------------------+------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=========================+========================================================================+==========+===========+ +| solve_species | Enable the MFIX species solver. | Int | 0 | ++-------------------------+------------------------------------------------------------------------+----------+-----------+ + Species definitions -=================== +------------------- The following inputs are defined using the prefix ``species``: -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+==========================+========================================================================+==========+===========+ -| solve | Specified name(s) of the species or None to disable the | Strings | None | -| | species solver. The name assigned to the species | | | -| | solver is used to specify species inputs. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity | Fluid species diffusivity model. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - a constant diffusion coefficient is used for | | | -| | all fluid species. | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | -| | ``constant`` diffusivity model | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat | Species specific heat model. This setting only applies if either | String | | -| | ``fluid.specific_heat = mixture`` or | | | -| | ``solids.specific_heat = mixture`` | | | -| | | | | -| | Options (case-insensitive): | | | -| | | | | -| | * ``constant`` - a constant specific heat is defined for each species | | | -| | and a *mixture* specific heat is computed based on the fluid or | | | -| | particle composition. | | | -| | | | | -| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | -| | or more polynomials that are a function of temperature, then the | | | -| | *mixture* specific heat is computed based on the fluid or particle | | | -| | composition. | | | -| | | | | -| | NASA7 polynomial format: | | | -| | | | | -| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| ignore_discontinuities | MFIX-Exa asserts that the low and high NASA-7 polynomials are | Int | 0 | -| | continuous at the 1000K temperature split. However, some polynomials | | | -| | do not have valid high temperature coefficients, for example, liquid | | | -| | water. By setting this keyword to 1, the simulation will only warn | | | -| | that the polynomials are discontinuous. Simulations that use | | | -| | discontinuous polynomials and cross between the 1000K temperature | | | -| | split may fail when computing temperature from enthalpy. | | | -| | Therefore, this option should only be used when the simulation is not | | | -| | expected to cross over the transition temperature. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | | -| | ``fluid.viscosity.molecular = mixture``. | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` | | | -| | * ``Reid`` :cite:p:`reid87` | | | -| | | | | -| | See :ref:`viscosity model descriptions ` | | | -| | in fluid section. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++======================================+=========================================================================+==========+===========+ +| solve | Specified name(s) of the species or None to disable the | Strings | None | +| | species solver. The name assigned to the species | | | +| | solver is used to specify species inputs. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity | Fluid species diffusivity model. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``constant`` - a constant diffusion coefficient is used for | | | +| | all fluid species. | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | +| | ``constant`` diffusivity model | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat | Species specific heat model. This setting only applies if either | String | | +| | ``fluid.specific_heat = mixture`` or | | | +| | ``solids.specific_heat = mixture`` | | | +| | | | | +| | Options (case-insensitive): | | | +| | | | | +| | * ``constant`` - a constant specific heat is defined for each species | | | +| | and a *mixture* specific heat is computed based on the fluid or | | | +| | particle composition. | | | +| | | | | +| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | +| | or more polynomials that are a function of temperature, then the | | | +| | *mixture* specific heat is computed based on the fluid or particle | | | +| | composition. | | | +| | | | | +| | NASA7 polynomial format: | | | +| | | | | +| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.ignore_discontinuities | MFIX-Exa checks that NASA-7 polynomials are continuous at the specified | Int | 0 | +| | transition temperature(s)(``Tsplit``). However, some polynomials do not | | | +| | have valid high temperature coefficients, for example, liquid water. By | | | +| | setting this keyword to 1, the simulation will only warn that the | | | +| | polynomials are discontinuous. Simulations that use discontinuous | | | +| | polynomials and cross ``Tsplit`` boundaries may fail when computing | | | +| | temperature from enthalpy. Therefore, this option should only be used | | | +| | when the simulation is not expected to cross over the transition | | | +| | temperature. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | None | +| | ``[fluid_name].viscosity.molecular`` = ``mixture``. | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` | | | +| | * ``Reid`` :cite:p:`reid87` | | | +| | | | | +| | See :ref:`viscosity model descriptions ` | | | +| | in fluid section. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ The following inputs are defined for each species using the prefix ``species.[species_name]``: -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+=======================================+========================================================================+==========+===========+ -| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | -| | model and when including chemical reactions. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | -| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | -| | specific heat model is ``constant`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Real | None | -| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | -| | and the specific heat model is ``NASA7-poly`` | | | -| | | | | -| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | -| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | -| | heat while all six coefficients are used to compute enthalpy. | | | -| | The seventh coefficient, ``a6``, may be provided for completeness | | | -| | but is not used as MFIX-Exa presently does not need to compute | | | -| | entropy. | | | -| | | | | -| | * By default, two sets of coefficients are required. The coefficients | | | -| | define two polynomials used to compute the specific heat and | | | -| | enthalpy across low and high temperatures. The split temperature | | | -| | (Tsplit) defines the transition from low-to-high polynomials. | | | -| | | | | -| | * A single polynomial can be used for all temperatures by setting the | | | -| | split temperature to -1. | | | -| | | | | -| | * Three or more polynomials can be used by defining multiple split | | | -| | temperatures. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.Tsplit | Defines the temperature transition between NASA-7 polynomials. | Real | 1000. | -| | | | | -| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | -| | the specific heat and enthalpy for all temperatures. | | | -| | | | | -| | * Specifying ``N`` split temperatures where ``N>1``, signifies ``N+1`` | | | -| | polynomials are used to define the specific heat and enthalpy. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | -| | the specific heat model is ``constant`` and chemical reactions | | | -| | are defined. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| density | Constant species density (only used for solid phase species). | Real | 0. | -| | | | | -| | Needed for all solid species if | | | -| | :ref:`solids density` is set to ``mixture`` or | | | -| | :ref:`density updates` is set to | | | -| | ``Constant_Species_Densities`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=======================================+=========================================================================+==========+===========+ +| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | +| | model and when including chemical reactions. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | +| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | +| | specific heat model is ``constant`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Reals | None | +| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | +| | and the specific heat model is ``NASA7-poly`` | | | +| | | | | +| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | +| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | +| | heat while all six coefficients are used to compute enthalpy. | | | +| | The seventh coefficient, ``a6``, may be provided for completeness | | | +| | but is not used as MFIX-Exa presently does not need to compute | | | +| | entropy. | | | +| | | | | +| | * By default, two sets of coefficients are required. The coefficients | | | +| | define two polynomials used to compute the specific heat and | | | +| | enthalpy across low and high temperatures. The transition temperature | | | +| | (``Tsplit``) defines the transition from low to high polynomials. | | | +| | | | | +| | * A single polynomial can be used for all temperatures by setting the | | | +| | transition temperature to -1. | | | +| | | | | +| | * Three or more polynomials can be used by defining multiple transition | | | +| | temperatures. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | +| | the specific heat and enthalpy for all temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` polynomials are | | | +| | used to define the specific heat and enthalpy. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | +| | the specific heat model is ``constant`` and chemical reactions | | | +| | are defined. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| density | Constant species density (only used for solid phase species). | Real | 0. | +| | | | | +| | Needed for all solid species if | | | +| | :ref:`solids density` is set to ``mixture`` or | | | +| | :ref:`density updates` is set to | | | +| | ``Constant_Species_Densities`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ Example inputs diff --git a/duplicate-keys.txt b/duplicate-keys.txt index efa155e..dd14b43 100644 --- a/duplicate-keys.txt +++ b/duplicate-keys.txt @@ -1,11 +1,4 @@ Duplicate keys: -ascent.actions: ascent.rst, plotting.rst -bc.[region_name].[fluid_name].massflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].velocity: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volfrac: boundary_conditions.rst, boundary_conditions.rst hypre.bamg_num_down_sweeps: multigrid_solvers.rst, multigrid_solvers.rst hypre.bamg_num_sweeps: multigrid_solvers.rst, multigrid_solvers.rst -mfix.ascent_int: ascent.rst, plotting.rst -mfix.ascent_per_approx: ascent.rst, plotting.rst mfix.restart: checkpointing.rst, initialization.rst diff --git a/nonexistent-keys.txt b/nonexistent-keys.txt index 72a209e..de176cc 100644 --- a/nonexistent-keys.txt +++ b/nonexistent-keys.txt @@ -1,4 +1,3 @@ -# This group may be incorrect bc.[region_name].[solid_name].bins bc.[region_name].[solid_name].density.custom bc.[region_name].[solid_name].density.max -- GitLab From fbf3a69e8d1698f332f9728d5366c26c4015981f Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:24:41 -0600 Subject: [PATCH 06/16] clarify --- docs/source_docs/user_guide/inputs/porous_media.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/porous_media.rst b/docs/source_docs/user_guide/inputs/porous_media.rst index 0de4984..6ab3c1f 100644 --- a/docs/source_docs/user_guide/inputs/porous_media.rst +++ b/docs/source_docs/user_guide/inputs/porous_media.rst @@ -3,7 +3,7 @@ Porous media definitions ======================== -Homogeneous porous media are modeled by the addition of a momentum +Homogeneous porous media are modeled by the introduction of a momentum source term, :math:`\boldsymbol{S}_{pm}`, to the fluid momentum equation. .. math:: -- GitLab From 6133692bdc3a214727053c466c2b654661439399 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:24:58 -0600 Subject: [PATCH 07/16] whitespace --- .../references/fluids/FluidEquations.rst | 61 +++++++++---------- 1 file changed, 29 insertions(+), 32 deletions(-) diff --git a/docs/source_docs/references/fluids/FluidEquations.rst b/docs/source_docs/references/fluids/FluidEquations.rst index e74da21..456cbb8 100644 --- a/docs/source_docs/references/fluids/FluidEquations.rst +++ b/docs/source_docs/references/fluids/FluidEquations.rst @@ -22,18 +22,18 @@ Conservation of fluid mass: .. math:: \frac{\partial (\varepsilon_g \rho_g)}{\partial t} + \nabla \cdot (\varepsilon_g \rho_g U_g) = 0 -Unlike two-fluid modeling, :math:`\varepsilon_g = 1 - \varepsilon_s`, is not a solution -variable. Rather, :math:`\varepsilon_s` is evaluated from the particle field through -volume filtering, +Unlike two-fluid modeling, :math:`\varepsilon_g = 1 - \varepsilon_s`, is not a solution +variable. Rather, :math:`\varepsilon_s` is evaluated from the particle field through +volume filtering, -.. math:: - (1 - \varepsilon_g) A (\boldsymbol{x},t) \approx - \sum_{i=1}^{N_p} A_i(\boldsymbol{X}_i,t) {\mathcal G} - (\left| \boldsymbol{x} - \boldsymbol{X}_i \right|) {\mathcal V}_i +.. math:: + (1 - \varepsilon_g) A (\boldsymbol{x},t) \approx + \sum_{i=1}^{N_p} A_i(\boldsymbol{X}_i,t) {\mathcal G} + (\left| \boldsymbol{x} - \boldsymbol{X}_i \right|) {\mathcal V}_i -where :math:`A_i` is a general particle property, :math:`{\mathcal V}_i` is the particle -volume and :math:`\mathcal G` is a transfer kernel with compact support--here linear hat. -Setting :math:`A_i = 1` gives the local void fraction. +where :math:`A_i` is a general particle property, :math:`{\mathcal V}_i` is the particle +volume and :math:`\mathcal G` is a transfer kernel with compact support--here linear hat. +Setting :math:`A_i = 1` gives the local void fraction. Assuming the fluid phase is incompressible, :math:`\frac{D \rho_g}{Dt} = 0`, the conservation of fluid mass is equivalent to the conservation of fluid volume: @@ -42,42 +42,39 @@ Assuming the fluid phase is incompressible, :math:`\frac{D \rho_g}{Dt} = 0`, the The conservation of fluid momentum is: -.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} +.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g = \nabla \cdot \tau + M_{sg} + \varepsilon_g \rho_g g -where :math:`M_{sg} = - M_{gs}` is the generalized interfacial momentum transfer from -the solid particles to the fluid-phase. Like :math:`\varepsilon_s`, -:math:`M_{gs}` is determined from the L-E transfer kernel by setting :math:`A_i = F_{gi}`, -where :math:`F_{gi}` is the force due to the fluid-phase on the ith particle. Following -MFiX classic (and many other CFD-DEM codes designed for high density ratio gas-solids -flows), only buoyancy (pressure gradient) and steady drag are considered: +where :math:`M_{sg} = - M_{gs}` is the generalized interfacial momentum transfer from +the solid particles to the fluid-phase. Like :math:`\varepsilon_s`, +:math:`M_{gs}` is determined from the L-E transfer kernel by setting :math:`A_i = F_{gi}`, +where :math:`F_{gi}` is the force due to the fluid-phase on the ith particle. Following +MFiX classic (and many other CFD-DEM codes designed for high density ratio gas-solids +flows), only buoyancy (pressure gradient) and steady drag are considered: .. math:: - F_{gi} = - \mathcal{V}_i \nabla p_g + F_{gi} = - \mathcal{V}_i \nabla p_g - \frac{1}{2} C_D \rho_g \boldsymbol{V}_{ig} \left|\boldsymbol{V}_{ig}\right| A_i^{(proj)} -where :math:`\boldsymbol{V}_{ig} = \boldsymbol{V}_i - \boldsymbol{U}_g ( \boldsymbol{X}_i )` -is the velocity of ith particle relative to the fluid-phase (at the particle position -:math:`\boldsymbol{X}_i`). :math:`F_{gi}` is closed by the specification of a drag -coefficient, :math:`C_D`. Currently, MFIX-Exa includes Wen-Yu, Gidaspow and BVK2 drag laws. +where :math:`\boldsymbol{V}_{ig} = \boldsymbol{V}_i - \boldsymbol{U}_g ( \boldsymbol{X}_i )` +is the velocity of ith particle relative to the fluid-phase (at the particle position +:math:`\boldsymbol{X}_i`). :math:`F_{gi}` is closed by the specification of a drag +coefficient, :math:`C_D`. Currently, MFIX-Exa includes Wen-Yu, Gidaspow, BVK2, and Syamlal-O'Brien drag laws. -.. wdf todo +.. wdf todo gidaspow form: C_D^{Gidaspow} = \chi C_D^{(Wen-Yu)} + \left(1 - \chi \right) C_D^{(Ergun)} \chi = \frac{\arctan 150 \left( \varepsilon_g - 0.8 \right)}{\pi} + \frac{1}{2} wen-yu: C_D^{(Wen-Yu)} = \max \left[\frac{24}{Re_i}\left(1 + 0.15 Re_i^{0.687}\right),\ 0.44 \right] \left( 1 - \varepsilon_g \right)^{-1.65} - ergun: C_D^{(Ergun)} = \frac{200 \left(1 - \varepsilon_g\right)}{Re_i} + \frac{7}{3} + ergun: C_D^{(Ergun)} = \frac{200 \left(1 - \varepsilon_g\right)}{Re_i} + \frac{7}{3} BVK2: C_D^{BVK2} = ... particle reynolds number: Re_i = \frac{\rho_g \left( 1 - \varepsilon_g \right) d_i \left| \boldsymbol{V}_{ig} \right| }{\mu_g} -In chemical engineering literature, it is common to lump all drag-related terms of -:math:`M_{gs}` into :math:`\beta`. With this simplification and some re-arrangement, -the fluid momentum takes the more convenient form: +In chemical engineering literature, it is common to lump all drag-related terms of +:math:`M_{gs}` into :math:`\beta`. With this simplification and some re-arrangement, +the fluid momentum takes the more convenient form: -.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} - + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g +.. math:: \frac{ \partial (\varepsilon_g \rho_g U_g)}{\partial t} + + \nabla \cdot (\varepsilon_g \rho_g U_g U_g) + \varepsilon_g \nabla p_g = \nabla \cdot \tau + \sum_p \beta_p (V_p - U_g) + \varepsilon_g \rho_g g - - - -- GitLab From 3b9199cb8f2b14d3a8978ee512ee2cf7c6014e5c Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Wed, 10 Dec 2025 13:25:08 -0600 Subject: [PATCH 08/16] spell out Syamlal --- docs/source_docs/user_guide/inputs/model_options.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 92ae7eb..7983264 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -318,7 +318,7 @@ The following input is defined using the prefix ``mfix.drag``: | | * ``WenYu`` :cite:p:`Wen66` | | | | | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | | | * ``BVK2`` :cite:p:`Bee07,Bee071,tenneti11,tang15` | | | -| | * ``SyamOBrien`` :cite:p:`syam88` | | | +| | * ``SyamOBrien`` -- Syamlal-O'Brien :cite:p:`syam88` | | | | | * ``UserDrag`` | | | | | | | | +--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -- GitLab From f7574d14334ba46afa08df4bef99f5d0cae7193f Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:07:45 -0600 Subject: [PATCH 09/16] move mfix.solve_species to species section --- docs/source_docs/user_guide/inputs/model_options.rst | 2 -- 1 file changed, 2 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 7983264..658b9c3 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -13,8 +13,6 @@ The following inputs are defined using the prefix ``mfix``: +------------------------+----------------------------------------------------------------------------+----------+---------------------+ | advect_enthalpy | Enable time evolution of fluid temperature and enthalpy. | Int | 0 | +------------------------+----------------------------------------------------------------------------+----------+---------------------+ -| solve_species | Enable the time evolution of fluid species mass fractions. | Int | 0 | -+------------------------+----------------------------------------------------------------------------+----------+---------------------+ | constraint | Select low Mach number constraint. | String | IncompressibleFluid | | | | | | | | Options (case-insensitive): | | | -- GitLab From 660ed8caaaf925915df5ecb30e838f63821d302a Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:11:56 -0600 Subject: [PATCH 10/16] update docs move mfix.solve_species to species section rename ignore_discontinuities to specific_heat.ignore_discontinuities --- .../user_guide/inputs/fluid_model.rst | 420 +++++++++--------- .../user_guide/inputs/initial_conditions.rst | 2 +- .../user_guide/inputs/model_options.rst | 79 ++-- .../user_guide/inputs/solids_model.rst | 172 ++++--- .../source_docs/user_guide/inputs/species.rst | 286 ++++++------ duplicate-keys.txt | 7 - nonexistent-keys.txt | 1 - 7 files changed, 516 insertions(+), 451 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 198894b..2b6d787 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -27,210 +27,226 @@ a single fluid; therefore, it is common to name the fluid ``fluid``. This is ill .. |Sato_Eq| replace:: :math:`\mu_{pit} = C d_s \rho \left|\boldsymbol{u} - \boldsymbol{u_s}\right|` .. |eff_visc| replace:: :math:`\mu_{eff} = \mu_{mol} + \mu_{eddy} + \mu_{susp} + \mu_{pit}` - .. |mix_Eq| replace:: :math:`\mu_{mix} = \sum_{\alpha=1}^{N} \frac{X_{\alpha} \mu_{\alpha}}{\sum_{\beta} X_{\beta} \phi_{\alpha \beta}}` -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| | Description | Type | Default | -+============================================+============================================================+========+==========+ -| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - no viscosity model. | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` - | | | -| | | | | -| | |Sutherland_Eq| | | | -| | | | | -| | * ``Reid`` :cite:p:`reid87` - | | | -| | | | | -| | |Reid_4parm_Eq| | | | -| | | | | -| | * ``mixture`` - a mixture viscosity is computed from | | | -| | :ref:`species viscosities` | | | -| | and local species mass fractions :cite:p:`bird06` | | | -| | | | | -| | |mix_Eq| | | | -| | | | | -| | A viscosity model is required if the fluid solver is | | | -| | enabled. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy | Specify eddy viscosity model. | | | -| | | | | -| | Options: | | | -| | |String | None | -| | * ``None`` - No eddy viscosity model | | | -| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | -| | * ``WALE`` :cite:p:`ducros98` | | | -| | * ``usr`` | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | -| | between 0.1 and 0.2. | | | -| | | | | -| | A value is required when using the ``Smagorinsky-Lilly`` | | | -| | eddy viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension | Specify suspension viscosity model of the form | String | None | -| | | | | -| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - No eddy suspension model | | | -| | | | | -| | :math:`\mu^* = 1` | | | -| | | | | -| | * ``Einstein`` :cite:p:`einstein11` | | | -| | | | | -| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | -| | | | | -| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | -| | | | | -| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | -| | | | | -| | * ``ChengLaw`` :cite:p:`cheng03` | | | -| | | | | -| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | -| | | | | -| | * ``Sato`` :cite:p:`sato81` | | | -| | | | | -| | |Sato_Eq| | | | -| | | | | -| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | -| | | | | -| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Brinkman`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``ChengLaw`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | -| | molecular viscosity | | | -| | | | | -| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | -| | |eff_visc| | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| species | Specify the species that constitute the fluid. | Strings| None | -| | | | | -| | All listed species must be properly defined. See the | | | -| | :ref:`species definition documentation` | | | -| | for additional details. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| molecular_weight | Constant fluid molecular weight. | Real | 0 | -| | | | | -| | A molecular weight should only be defined when using one | | | -| | of the ``IdealGas`` constraints and the species equations | | | -| | are not solved. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat | Specify which specific heat model to use for the fluid. | String | None | -| | Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant specific heat | | | -| | | | | -| | * ``NASA7-poly`` the fluid 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}` | | | -| | | | | -| | A specific heat model is required if advecting enthalpy. | | | -| | Furthermore, the model must be ``mixture`` if species | | | -| | equations are solved. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.constant | Constant fluid specific heat | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.NASA7.a[i] | Specific heat polynomial coefficients. Polynomial | Real | 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`. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | -| | fluid. Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant thermal | | | -| | conductivity model | | | -| | | | | -| | A thermal conductivity model is required if advecting | | | -| | enthalpy. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | -| | A value is required when using the | | | -| | ``IdealGasClosedSystem`` constraint. If using the | | | -| | ``IdealGasOpenSystem`` constraint, then the value must | | | -| | match all specified outflow pressures. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| reference_temperature | Reference temperature used to compute specific | Real | 298. | -| | enthalpy | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| | Description | Type | Default | ++============================================+================================================================+========+==========+ +| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model. | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` - | | | +| | | | | +| | |Sutherland_Eq| | | | +| | | | | +| | * ``Reid`` :cite:p:`reid87` - | | | +| | | | | +| | |Reid_4parm_Eq| | | | +| | | | | +| | * ``mixture`` - a mixture viscosity is computed from | | | +| | :ref:`species viscosities` | | | +| | and local species mass fractions :cite:p:`bird06` | | | +| | | | | +| | |mix_Eq| | | | +| | | | | +| | A viscosity model is required if the fluid solver is | | | +| | enabled. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy | Specify eddy viscosity model. | | | +| | | | | +| | Options: | | | +| | |String | None | +| | * ``None`` - No eddy viscosity model | | | +| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | +| | * ``WALE`` :cite:p:`ducros98` | | | +| | * ``usr`` | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | +| | between 0.1 and 0.2. | | | +| | | | | +| | A value is required when using the ``Smagorinsky-Lilly`` | | | +| | eddy viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension | Specify suspension viscosity model of the form | String | None | +| | | | | +| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - No eddy suspension model | | | +| | | | | +| | :math:`\mu^* = 1` | | | +| | | | | +| | * ``Einstein`` :cite:p:`einstein11` | | | +| | | | | +| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | +| | | | | +| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | +| | | | | +| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | +| | | | | +| | * ``ChengLaw`` :cite:p:`cheng03` | | | +| | | | | +| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | +| | | | | +| | * ``Sato`` :cite:p:`sato81` | | | +| | | | | +| | |Sato_Eq| | | | +| | | | | +| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | +| | | | | +| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Brinkman`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``ChengLaw`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | +| | molecular viscosity | | | +| | | | | +| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | +| | |eff_visc| | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| species | Specify the species that constitute the fluid. | Strings| None | +| | | | | +| | All listed species must be properly defined. See the | | | +| | :ref:`species definition documentation` | | | +| | for additional details. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| molecular_weight | Constant fluid molecular weight. | Real | 0 | +| | | | | +| | A molecular weight should only be defined when using one | | | +| | of the ``IdealGas`` constraints and the species equations | | | +| | are not solved. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat | Specify which specific heat model to use for the fluid. | String | None | +| | Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant specific heat | | | +| | | | | +| | * ``NASA7-poly`` the fluid 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.constant | Constant fluid specific heat | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provided | | | +| | in :ref:`species specific heats`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.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`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single | | | +| | polynomial defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat | | | +| | and enthalpy. | | | +| | | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | +| | fluid. Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant thermal | | | +| | conductivity model | | | +| | | | | +| | A thermal conductivity model is required if advecting | | | +| | enthalpy. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | +| | A value is required when using the | | | +| | ``IdealGasClosedSystem`` constraint. If using the | | | +| | ``IdealGasOpenSystem`` constraint, then the value must | | | +| | match all specified outflow pressures. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| reference_temperature | Reference temperature used to compute specific | Real | 298. | +| | enthalpy | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ Example inputs diff --git a/docs/source_docs/user_guide/inputs/initial_conditions.rst b/docs/source_docs/user_guide/inputs/initial_conditions.rst index 579419c..5c1397c 100644 --- a/docs/source_docs/user_guide/inputs/initial_conditions.rst +++ b/docs/source_docs/user_guide/inputs/initial_conditions.rst @@ -20,7 +20,7 @@ The following inputs are defined using the prefix ``ic``: | | * ``Volume`` - sort by volume of each IC region | | | | | * ``Priority`` - sort by specified priority value | | | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ -| [region_name].priority| Priority value for IC regions if ``ranking_type` == ``Priority`` | Int | Max | +| [region_name].priority| Priority value for IC regions if ``ranking_type`` == ``Priority`` | Int | Max | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 658b9c3..f5c0338 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -308,45 +308,46 @@ The following input is defined using the prefix ``mfix.drag``: .. |VM_Nijssen_Eq| replace:: :math:`\begin{split}C_p,vm =& C_0,vm \\ &+ \sum_{k=1}^3 \left\left[a_{k,0} + a_{k,1} \ln\left(\frac{\rho_p}{\rho_f}\right)\right] \varepsilon_p^k \right\end{split}` -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| | Description | Type | Default | -+==========================+===============================================================================+=============+===========+ -| model | Fluid-particle drag model. Options: | String | None | -| | | | | -| | * ``WenYu`` :cite:p:`Wen66` | | | -| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | -| | * ``BVK2`` :cite:p:`Bee07,Bee071,tenneti11,tang15` | | | -| | * ``SyamOBrien`` -- Syamlal-O'Brien :cite:p:`syam88` | | | -| | * ``UserDrag`` | | | -| | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | -| model.SyamOBrien.d1 | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | -| | fluid-particle drag force. The force is applied to the fluid by | | | -| | multiplying the shear stress by fluid volume fraction. | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | -| | The virtual mass force is not included by default. | | | -| | | | | -| | |VirtualMass_Eq| | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - virtual mass force calculation is skipped | | | -| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | -| | * ``Constant`` - specified constant coefficient | | | -| | * ``Zuber`` - :cite:p:`zuber64` | | | -| | | | | -| | |VM_Zuber_Eq| | | | -| | | | | -| | * ``Nijssen`` - :cite:p:`nijssen20` | | | -| | | | | -| | |VM_Nijssen_Eq| | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| | Description | Type | Default | ++==========================+=========================================================================================+=============+===========+ +| model | Fluid-particle drag model. Options: | String | None | +| | | | | +| | * ``WenYu``- Wen-Yu :cite:p:`Wen66` | | | +| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | +| | * ``BVK2`` - Beetstra, van der Hoef, and Kuipers :cite:p:`Bee07,Bee071,tenneti11,tang15`| | | +| | * ``SyamOBrien`` - Syamlal-O'Brien :cite:p:`syam88` | | | +| | * ``UserDrag`` - User-supplied drag model | | | +| | | | | +| | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | +| model.SyamOBrien.d1 | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | +| | fluid-particle drag force. The force is applied to the fluid by | | | +| | multiplying the shear stress by fluid volume fraction. | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | +| | The virtual mass force is not included by default. | | | +| | | | | +| | |VirtualMass_Eq| | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - virtual mass force calculation is skipped | | | +| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | +| | * ``Constant`` - specified constant coefficient | | | +| | * ``Zuber`` - :cite:p:`zuber64` | | | +| | | | | +| | |VM_Zuber_Eq| | | | +| | | | | +| | * ``Nijssen`` - :cite:p:`nijssen20` | | | +| | | | | +| | |VM_Nijssen_Eq| | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ .. note:: diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index 1436464..9ad6278 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -9,71 +9,113 @@ Solids settings 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| molecular_weight | Value of constant solid molecular weight. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| density | Specify which density model to use for solids. | String | None | -| | | | | -| | Options (case-insensive): | | | -| | | | | -| | * ``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. | | | -| | | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat | Specify which specific heat model to use for solids. | String | constant | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant specific heat | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat.constant | Constant solids specific heat value. | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| species | Specify which species can constitute | Strings | None | -| | the solids phase. Defined species must | | | -| | be a subset of the species defined in ``species.solve`` | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | -| | Options: | | | -| | | | | -| | * ``constant`` - constant thermal conductivity | | | -| | | | | -| | A thermal conductivity model is not required for solids. | | | -| | Defining a thermal conductivity model enables conductive | | | -| | heat transfer for DEM. However, it has no effect for PIC | | | -| | solids heat transfer. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| | 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| molecular_weight | Value of constant solid molecular weight. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| density | Specify which density model to use for solids. | String | None | +| | | | | +| | Options (case-insensive): | | | +| | | | | +| | * ``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. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat | Specify which specific heat model to use for solids. | String | constant | +| | | | | +| | 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.constant | Constant solids specific heat value. | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provied in | | | +| | :ref:`species specific heats`. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.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`. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 | Reals | 1000 | +| | polynomials. | | | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial | | | +| | defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat and | | | +| | enthalpy. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| species | Specify which species can constitute | Strings | None | +| | the solids phase. Defined species must | | | +| | be a subset of the species defined in ``species.solve`` | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | +| | Options: | | | +| | | | | +| | * ``constant`` - constant thermal conductivity | | | +| | | | | +| | A thermal conductivity model is not required for solids. | | | +| | Defining a thermal conductivity model enables conductive | | | +| | heat transfer for DEM. However, it has no effect for PIC | | | +| | solids heat transfer. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ Below is an example for specifying the solids solver model options. diff --git a/docs/source_docs/user_guide/inputs/species.rst b/docs/source_docs/user_guide/inputs/species.rst index 17037ad..f02aa4e 100644 --- a/docs/source_docs/user_guide/inputs/species.rst +++ b/docs/source_docs/user_guide/inputs/species.rst @@ -1,149 +1,163 @@ .. _InputsSpeciesDefs: + +Species solver +============== + +Enabling the species solver +--------------------------- + +The following input is defined using the prefix ``mfix``: + ++-------------------------+------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=========================+========================================================================+==========+===========+ +| solve_species | Enable the MFIX species solver. | Int | 0 | ++-------------------------+------------------------------------------------------------------------+----------+-----------+ + Species definitions -=================== +------------------- The following inputs are defined using the prefix ``species``: -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+==========================+========================================================================+==========+===========+ -| solve | Specified name(s) of the species or None to disable the | Strings | None | -| | species solver. The name assigned to the species | | | -| | solver is used to specify species inputs. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity | Fluid species diffusivity model. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - a constant diffusion coefficient is used for | | | -| | all fluid species. | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | -| | ``constant`` diffusivity model | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat | Species specific heat model. This setting only applies if either | String | | -| | ``fluid.specific_heat = mixture`` or | | | -| | ``solids.specific_heat = mixture`` | | | -| | | | | -| | Options (case-insensitive): | | | -| | | | | -| | * ``constant`` - a constant specific heat is defined for each species | | | -| | and a *mixture* specific heat is computed based on the fluid or | | | -| | particle composition. | | | -| | | | | -| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | -| | or more polynomials that are a function of temperature, then the | | | -| | *mixture* specific heat is computed based on the fluid or particle | | | -| | composition. | | | -| | | | | -| | NASA7 polynomial format: | | | -| | | | | -| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| ignore_discontinuities | MFIX-Exa asserts that the low and high NASA-7 polynomials are | Int | 0 | -| | continuous at the 1000K temperature split. However, some polynomials | | | -| | do not have valid high temperature coefficients, for example, liquid | | | -| | water. By setting this keyword to 1, the simulation will only warn | | | -| | that the polynomials are discontinuous. Simulations that use | | | -| | discontinuous polynomials and cross between the 1000K temperature | | | -| | split may fail when computing temperature from enthalpy. | | | -| | Therefore, this option should only be used when the simulation is not | | | -| | expected to cross over the transition temperature. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | | -| | ``fluid.viscosity.molecular = mixture``. | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` | | | -| | * ``Reid`` :cite:p:`reid87` | | | -| | | | | -| | See :ref:`viscosity model descriptions ` | | | -| | in fluid section. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++======================================+=========================================================================+==========+===========+ +| solve | Specified name(s) of the species or None to disable the | Strings | None | +| | species solver. The name assigned to the species | | | +| | solver is used to specify species inputs. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity | Fluid species diffusivity model. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``constant`` - a constant diffusion coefficient is used for | | | +| | all fluid species. | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | +| | ``constant`` diffusivity model | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat | Species specific heat model. This setting only applies if either | String | | +| | ``fluid.specific_heat = mixture`` or | | | +| | ``solids.specific_heat = mixture`` | | | +| | | | | +| | Options (case-insensitive): | | | +| | | | | +| | * ``constant`` - a constant specific heat is defined for each species | | | +| | and a *mixture* specific heat is computed based on the fluid or | | | +| | particle composition. | | | +| | | | | +| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | +| | or more polynomials that are a function of temperature, then the | | | +| | *mixture* specific heat is computed based on the fluid or particle | | | +| | composition. | | | +| | | | | +| | NASA7 polynomial format: | | | +| | | | | +| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.ignore_discontinuities | MFIX-Exa checks that NASA-7 polynomials are continuous at the specified | Int | 0 | +| | transition temperature(s)(``Tsplit``). However, some polynomials do not | | | +| | have valid high temperature coefficients, for example, liquid water. By | | | +| | setting this keyword to 1, the simulation will only warn that the | | | +| | polynomials are discontinuous. Simulations that use discontinuous | | | +| | polynomials and cross ``Tsplit`` boundaries may fail when computing | | | +| | temperature from enthalpy. Therefore, this option should only be used | | | +| | when the simulation is not expected to cross over the transition | | | +| | temperature. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | None | +| | ``[fluid_name].viscosity.molecular`` = ``mixture``. | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` | | | +| | * ``Reid`` :cite:p:`reid87` | | | +| | | | | +| | See :ref:`viscosity model descriptions ` | | | +| | in fluid section. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ The following inputs are defined for each species using the prefix ``species.[species_name]``: -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+=======================================+========================================================================+==========+===========+ -| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | -| | model and when including chemical reactions. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | -| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | -| | specific heat model is ``constant`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Real | None | -| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | -| | and the specific heat model is ``NASA7-poly`` | | | -| | | | | -| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | -| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | -| | heat while all six coefficients are used to compute enthalpy. | | | -| | The seventh coefficient, ``a6``, may be provided for completeness | | | -| | but is not used as MFIX-Exa presently does not need to compute | | | -| | entropy. | | | -| | | | | -| | * By default, two sets of coefficients are required. The coefficients | | | -| | define two polynomials used to compute the specific heat and | | | -| | enthalpy across low and high temperatures. The split temperature | | | -| | (Tsplit) defines the transition from low-to-high polynomials. | | | -| | | | | -| | * A single polynomial can be used for all temperatures by setting the | | | -| | split temperature to -1. | | | -| | | | | -| | * Three or more polynomials can be used by defining multiple split | | | -| | temperatures. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.Tsplit | Defines the temperature transition between NASA-7 polynomials. | Real | 1000. | -| | | | | -| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | -| | the specific heat and enthalpy for all temperatures. | | | -| | | | | -| | * Specifying ``N`` split temperatures where ``N>1``, signifies ``N+1`` | | | -| | polynomials are used to define the specific heat and enthalpy. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | -| | the specific heat model is ``constant`` and chemical reactions | | | -| | are defined. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| density | Constant species density (only used for solid phase species). | Real | 0. | -| | | | | -| | Needed for all solid species if | | | -| | :ref:`solids density` is set to ``mixture`` or | | | -| | :ref:`density updates` is set to | | | -| | ``Constant_Species_Densities`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=======================================+=========================================================================+==========+===========+ +| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | +| | model and when including chemical reactions. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | +| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | +| | specific heat model is ``constant`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Reals | None | +| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | +| | and the specific heat model is ``NASA7-poly`` | | | +| | | | | +| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | +| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | +| | heat while all six coefficients are used to compute enthalpy. | | | +| | The seventh coefficient, ``a6``, may be provided for completeness | | | +| | but is not used as MFIX-Exa presently does not need to compute | | | +| | entropy. | | | +| | | | | +| | * By default, two sets of coefficients are required. The coefficients | | | +| | define two polynomials used to compute the specific heat and | | | +| | enthalpy across low and high temperatures. The transition temperature | | | +| | (``Tsplit``) defines the transition from low to high polynomials. | | | +| | | | | +| | * A single polynomial can be used for all temperatures by setting the | | | +| | transition temperature to -1. | | | +| | | | | +| | * Three or more polynomials can be used by defining multiple transition | | | +| | temperatures. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | +| | the specific heat and enthalpy for all temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` polynomials are | | | +| | used to define the specific heat and enthalpy. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | +| | the specific heat model is ``constant`` and chemical reactions | | | +| | are defined. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| density | Constant species density (only used for solid phase species). | Real | 0. | +| | | | | +| | Needed for all solid species if | | | +| | :ref:`solids density` is set to ``mixture`` or | | | +| | :ref:`density updates` is set to | | | +| | ``Constant_Species_Densities`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ Example inputs diff --git a/duplicate-keys.txt b/duplicate-keys.txt index efa155e..dd14b43 100644 --- a/duplicate-keys.txt +++ b/duplicate-keys.txt @@ -1,11 +1,4 @@ Duplicate keys: -ascent.actions: ascent.rst, plotting.rst -bc.[region_name].[fluid_name].massflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].velocity: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volfrac: boundary_conditions.rst, boundary_conditions.rst hypre.bamg_num_down_sweeps: multigrid_solvers.rst, multigrid_solvers.rst hypre.bamg_num_sweeps: multigrid_solvers.rst, multigrid_solvers.rst -mfix.ascent_int: ascent.rst, plotting.rst -mfix.ascent_per_approx: ascent.rst, plotting.rst mfix.restart: checkpointing.rst, initialization.rst diff --git a/nonexistent-keys.txt b/nonexistent-keys.txt index 72a209e..de176cc 100644 --- a/nonexistent-keys.txt +++ b/nonexistent-keys.txt @@ -1,4 +1,3 @@ -# This group may be incorrect bc.[region_name].[solid_name].bins bc.[region_name].[solid_name].density.custom bc.[region_name].[solid_name].density.max -- GitLab From 20d8271ba37428926a818c3c67e5cf7bd90c0388 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:19:05 -0600 Subject: [PATCH 11/16] spelling --- docs/source_docs/user_guide/inputs/solids_model.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index 9ad6278..53669a4 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -59,7 +59,7 @@ models. The following inputs are defined with the prefix ``solids``: | | A value is required for ``constant`` specific heat model. | | | +------------------------------------------+--------------------------------------------------------------+----------+----------+ | specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | -| | temperatures (``Tsplit``). More information is provied in | | | +| | temperatures (``Tsplit``). More information is provided in | | | | | :ref:`species specific heats`. | | | | | | | | +------------------------------------------+--------------------------------------------------------------+----------+----------+ -- GitLab From 2b062d589c18cfea63c290bb6bd6e38aa25ed5bb Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:22:50 -0600 Subject: [PATCH 12/16] bullet lists --- docs/source_docs/user_guide/inputs/fluid_model.rst | 9 +++++---- docs/source_docs/user_guide/inputs/solids_model.rst | 7 ++++--- docs/source_docs/user_guide/inputs/species.rst | 1 + 3 files changed, 10 insertions(+), 7 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 2b6d787..dd689c5 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -169,6 +169,11 @@ a single fluid; therefore, it is common to name the fluid ``fluid``. This is ill | | are not solved. | | | +--------------------------------------------+----------------------------------------------------------------+--------+----------+ | specific_heat | Specify which specific heat model to use for the fluid. | String | None | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | +| | | | | | | Options: | | | | | | | | | | * ``constant`` - the fluid has a constant specific heat | | | @@ -189,10 +194,6 @@ a single fluid; therefore, it is common to name the fluid ``fluid``. This is ill | | | | | | | :math:`c_{p,\mathrm{mixture}} = \sum_n X_n c_{p,n}` | | | | | | | | -| | A specific heat model is required if advecting enthalpy. | | | -| | Furthermore, the model must be ``mixture`` if species | | | -| | equations are solved. | | | -| | | | | +--------------------------------------------+----------------------------------------------------------------+--------+----------+ | specific_heat.constant | Constant fluid specific heat | Real | 0 | | | A value is required for ``constant`` specific heat model. | | | diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index 53669a4..bb68fcc 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -32,6 +32,10 @@ models. The following inputs are defined with the prefix ``solids``: +------------------------------------------+--------------------------------------------------------------+----------+----------+ | specific_heat | 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 | | | @@ -51,9 +55,6 @@ models. The following inputs are defined with the prefix ``solids``: | | | | | | | :math:`c_{p,\mathrm{mixture}} = \sum_n X_n c_{p,n}` | | | | | | | | -| | A specific heat model is required if advecting enthalpy. | | | -| | Furthermore, the model must be ``mixture`` if species | | | -| | equations are solved. | | | +------------------------------------------+--------------------------------------------------------------+----------+----------+ | specific_heat.constant | Constant solids specific heat value. | Real | 0 | | | A value is required for ``constant`` specific heat model. | | | diff --git a/docs/source_docs/user_guide/inputs/species.rst b/docs/source_docs/user_guide/inputs/species.rst index f02aa4e..56af28c 100644 --- a/docs/source_docs/user_guide/inputs/species.rst +++ b/docs/source_docs/user_guide/inputs/species.rst @@ -125,6 +125,7 @@ The following inputs are defined for each species using the prefix ``species.[sp | | | | | | | * Specifying ``N`` temperatures signifies that ``N+1`` polynomials are | | | | | used to define the specific heat and enthalpy. | | | +| | | | | +---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | | | the specific heat model is ``constant`` and chemical reactions | | | -- GitLab From b9fa99f7884a1b6e4c6cf6c299c751dd4927614f Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:26:27 -0600 Subject: [PATCH 13/16] bullet lists --- docs/source_docs/user_guide/inputs/species.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/species.rst b/docs/source_docs/user_guide/inputs/species.rst index 56af28c..ebd5c72 100644 --- a/docs/source_docs/user_guide/inputs/species.rst +++ b/docs/source_docs/user_guide/inputs/species.rst @@ -124,7 +124,7 @@ The following inputs are defined for each species using the prefix ``species.[sp | | the specific heat and enthalpy for all temperatures. | | | | | | | | | | * Specifying ``N`` temperatures signifies that ``N+1`` polynomials are | | | -| | used to define the specific heat and enthalpy. | | | +| | used to define the specific heat and enthalpy. | | | | | | | | +---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ | enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | -- GitLab From a5f4782d8b02d4f6882116500fbfb0fc53b090b8 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:07:45 -0600 Subject: [PATCH 14/16] move mfix.solve_species to species section --- docs/source_docs/user_guide/inputs/model_options.rst | 2 -- 1 file changed, 2 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 7983264..658b9c3 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -13,8 +13,6 @@ The following inputs are defined using the prefix ``mfix``: +------------------------+----------------------------------------------------------------------------+----------+---------------------+ | advect_enthalpy | Enable time evolution of fluid temperature and enthalpy. | Int | 0 | +------------------------+----------------------------------------------------------------------------+----------+---------------------+ -| solve_species | Enable the time evolution of fluid species mass fractions. | Int | 0 | -+------------------------+----------------------------------------------------------------------------+----------+---------------------+ | constraint | Select low Mach number constraint. | String | IncompressibleFluid | | | | | | | | Options (case-insensitive): | | | -- GitLab From 8b4722e3b3d7e24b2e94ee486a0738430b90f73d Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:11:56 -0600 Subject: [PATCH 15/16] update docs move mfix.solve_species to species section rename ignore_discontinuities to specific_heat.ignore_discontinuities --- .../user_guide/inputs/fluid_model.rst | 420 +++++++++--------- .../user_guide/inputs/initial_conditions.rst | 2 +- .../user_guide/inputs/model_options.rst | 79 ++-- .../user_guide/inputs/solids_model.rst | 172 ++++--- .../source_docs/user_guide/inputs/species.rst | 286 ++++++------ duplicate-keys.txt | 7 - nonexistent-keys.txt | 1 - 7 files changed, 516 insertions(+), 451 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 198894b..2b6d787 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -27,210 +27,226 @@ a single fluid; therefore, it is common to name the fluid ``fluid``. This is ill .. |Sato_Eq| replace:: :math:`\mu_{pit} = C d_s \rho \left|\boldsymbol{u} - \boldsymbol{u_s}\right|` .. |eff_visc| replace:: :math:`\mu_{eff} = \mu_{mol} + \mu_{eddy} + \mu_{susp} + \mu_{pit}` - .. |mix_Eq| replace:: :math:`\mu_{mix} = \sum_{\alpha=1}^{N} \frac{X_{\alpha} \mu_{\alpha}}{\sum_{\beta} X_{\beta} \phi_{\alpha \beta}}` -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| | Description | Type | Default | -+============================================+============================================================+========+==========+ -| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - no viscosity model. | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` - | | | -| | | | | -| | |Sutherland_Eq| | | | -| | | | | -| | * ``Reid`` :cite:p:`reid87` - | | | -| | | | | -| | |Reid_4parm_Eq| | | | -| | | | | -| | * ``mixture`` - a mixture viscosity is computed from | | | -| | :ref:`species viscosities` | | | -| | and local species mass fractions :cite:p:`bird06` | | | -| | | | | -| | |mix_Eq| | | | -| | | | | -| | A viscosity model is required if the fluid solver is | | | -| | enabled. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy | Specify eddy viscosity model. | | | -| | | | | -| | Options: | | | -| | |String | None | -| | * ``None`` - No eddy viscosity model | | | -| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | -| | * ``WALE`` :cite:p:`ducros98` | | | -| | * ``usr`` | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | -| | between 0.1 and 0.2. | | | -| | | | | -| | A value is required when using the ``Smagorinsky-Lilly`` | | | -| | eddy viscosity model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension | Specify suspension viscosity model of the form | String | None | -| | | | | -| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - No eddy suspension model | | | -| | | | | -| | :math:`\mu^* = 1` | | | -| | | | | -| | * ``Einstein`` :cite:p:`einstein11` | | | -| | | | | -| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | -| | | | | -| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | -| | | | | -| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | -| | | | | -| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | -| | | | | -| | * ``ChengLaw`` :cite:p:`cheng03` | | | -| | | | | -| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | -| | | | | -| | * ``Sato`` :cite:p:`sato81` | | | -| | | | | -| | |Sato_Eq| | | | -| | | | | -| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | -| | | | | -| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Brinkman`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``Roscoe`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | -| | | | | -| | A value is required when using the ``ChengLaw`` model. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | -| | molecular viscosity | | | -| | | | | -| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | -| | |eff_visc| | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| species | Specify the species that constitute the fluid. | Strings| None | -| | | | | -| | All listed species must be properly defined. See the | | | -| | :ref:`species definition documentation` | | | -| | for additional details. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| molecular_weight | Constant fluid molecular weight. | Real | 0 | -| | | | | -| | A molecular weight should only be defined when using one | | | -| | of the ``IdealGas`` constraints and the species equations | | | -| | are not solved. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat | Specify which specific heat model to use for the fluid. | String | None | -| | Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant specific heat | | | -| | | | | -| | * ``NASA7-poly`` the fluid 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}` | | | -| | | | | -| | A specific heat model is required if advecting enthalpy. | | | -| | Furthermore, the model must be ``mixture`` if species | | | -| | equations are solved. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.constant | Constant fluid specific heat | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| specific_heat.NASA7.a[i] | Specific heat polynomial coefficients. Polynomial | Real | 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`. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | -| | fluid. Options: | | | -| | | | | -| | * ``constant`` - the fluid has a constant thermal | | | -| | conductivity model | | | -| | | | | -| | A thermal conductivity model is required if advecting | | | -| | enthalpy. | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | -| | A value is required when using the | | | -| | ``IdealGasClosedSystem`` constraint. If using the | | | -| | ``IdealGasOpenSystem`` constraint, then the value must | | | -| | match all specified outflow pressures. | | | -| | | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| reference_temperature | Reference temperature used to compute specific | Real | 298. | -| | enthalpy | | | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| | Description | Type | Default | ++============================================+================================================================+========+==========+ +| viscosity.molecular | Specify which molecular viscosity model to use. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model. | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` - | | | +| | | | | +| | |Sutherland_Eq| | | | +| | | | | +| | * ``Reid`` :cite:p:`reid87` - | | | +| | | | | +| | |Reid_4parm_Eq| | | | +| | | | | +| | * ``mixture`` - a mixture viscosity is computed from | | | +| | :ref:`species viscosities` | | | +| | and local species mass fractions :cite:p:`bird06` | | | +| | | | | +| | |mix_Eq| | | | +| | | | | +| | A viscosity model is required if the fluid solver is | | | +| | enabled. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.constant | Constant fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.molecular.Reid.A | Reid model constants. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy | Specify eddy viscosity model. | | | +| | | | | +| | Options: | | | +| | |String | None | +| | * ``None`` - No eddy viscosity model | | | +| | * ``Smagorinsky-Lilly`` :cite:p:`Smag63,Lilly66` | | | +| | * ``WALE`` :cite:p:`ducros98` | | | +| | * ``usr`` | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.Smagorinsky-Lilly.constant | Smagorinsky-Lilly constant which usually has values | Real | None | +| | between 0.1 and 0.2. | | | +| | | | | +| | A value is required when using the ``Smagorinsky-Lilly`` | | | +| | eddy viscosity model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension | Specify suspension viscosity model of the form | String | None | +| | | | | +| | :math:`\mu_{susp} = \mu_{mol}(\mu^* - 1)` | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - No eddy suspension model | | | +| | | | | +| | :math:`\mu^* = 1` | | | +| | | | | +| | * ``Einstein`` :cite:p:`einstein11` | | | +| | | | | +| | :math:`\mu^* = 1 + 2.5\varepsilon_s` | | | +| | | | | +| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s)^{-c}` | | | +| | | | | +| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | +| | | | | +| | :math:`\mu^* = (1-\varepsilon_s/c_1)^{-c_2}` | | | +| | | | | +| | * ``ChengLaw`` :cite:p:`cheng03` | | | +| | | | | +| | :math:`\mu^* = e^{2.5(1/(1-\varepsilon_s)^c-1)/c}` | | | +| | | | | +| | * ``Sato`` :cite:p:`sato81` | | | +| | | | | +| | |Sato_Eq| | | | +| | | | | +| | * ``Subramaniam`` :cite:p:`mehrabadi15` | | | +| | | | | +| | :math:`\mu_{pit} = C_\mu k_f^2 / \epsilon_f` | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Brinkman`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``Roscoe`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | +| | | | | +| | A value is required when using the ``ChengLaw`` model. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Sato.constant | Constant for Sato suspension expression. | Real | 0.65 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Subramaniam.constant | Constant for Subramaniam suspension expression. | Real | 0.09 | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| viscosity.max_effective_factor | Max limit of the effective viscosity as a factor of the | Real | 1.e6 | +| | molecular viscosity | | | +| | | | | +| | :math:`\mu_{eff} < fac*\mu_{mol}`, where | | | +| | |eff_visc| | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| species | Specify the species that constitute the fluid. | Strings| None | +| | | | | +| | All listed species must be properly defined. See the | | | +| | :ref:`species definition documentation` | | | +| | for additional details. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| molecular_weight | Constant fluid molecular weight. | Real | 0 | +| | | | | +| | A molecular weight should only be defined when using one | | | +| | of the ``IdealGas`` constraints and the species equations | | | +| | are not solved. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat | Specify which specific heat model to use for the fluid. | String | None | +| | Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant specific heat | | | +| | | | | +| | * ``NASA7-poly`` the fluid 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.constant | Constant fluid specific heat | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provided | | | +| | in :ref:`species specific heats`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.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`. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single | | | +| | polynomial defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat | | | +| | and enthalpy. | | | +| | | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for | String | None | +| | fluid. Options: | | | +| | | | | +| | * ``constant`` - the fluid has a constant thermal | | | +| | conductivity model | | | +| | | | | +| | A thermal conductivity model is required if advecting | | | +| | enthalpy. | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermal_conductivity.constant | Constant fluid thermal conductivity. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| thermodynamic_pressure | Thermodynamic pressure. | Real | 0 | +| | A value is required when using the | | | +| | ``IdealGasClosedSystem`` constraint. If using the | | | +| | ``IdealGasOpenSystem`` constraint, then the value must | | | +| | match all specified outflow pressures. | | | +| | | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ +| reference_temperature | Reference temperature used to compute specific | Real | 298. | +| | enthalpy | | | ++--------------------------------------------+----------------------------------------------------------------+--------+----------+ Example inputs diff --git a/docs/source_docs/user_guide/inputs/initial_conditions.rst b/docs/source_docs/user_guide/inputs/initial_conditions.rst index 579419c..5c1397c 100644 --- a/docs/source_docs/user_guide/inputs/initial_conditions.rst +++ b/docs/source_docs/user_guide/inputs/initial_conditions.rst @@ -20,7 +20,7 @@ The following inputs are defined using the prefix ``ic``: | | * ``Volume`` - sort by volume of each IC region | | | | | * ``Priority`` - sort by specified priority value | | | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ -| [region_name].priority| Priority value for IC regions if ``ranking_type` == ``Priority`` | Int | Max | +| [region_name].priority| Priority value for IC regions if ``ranking_type`` == ``Priority`` | Int | Max | +-----------------------+------------------------------------------------------------------------+-------------+-----------+ diff --git a/docs/source_docs/user_guide/inputs/model_options.rst b/docs/source_docs/user_guide/inputs/model_options.rst index 658b9c3..f5c0338 100644 --- a/docs/source_docs/user_guide/inputs/model_options.rst +++ b/docs/source_docs/user_guide/inputs/model_options.rst @@ -308,45 +308,46 @@ The following input is defined using the prefix ``mfix.drag``: .. |VM_Nijssen_Eq| replace:: :math:`\begin{split}C_p,vm =& C_0,vm \\ &+ \sum_{k=1}^3 \left\left[a_{k,0} + a_{k,1} \ln\left(\frac{\rho_p}{\rho_f}\right)\right] \varepsilon_p^k \right\end{split}` -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| | Description | Type | Default | -+==========================+===============================================================================+=============+===========+ -| model | Fluid-particle drag model. Options: | String | None | -| | | | | -| | * ``WenYu`` :cite:p:`Wen66` | | | -| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | -| | * ``BVK2`` :cite:p:`Bee07,Bee071,tenneti11,tang15` | | | -| | * ``SyamOBrien`` -- Syamlal-O'Brien :cite:p:`syam88` | | | -| | * ``UserDrag`` | | | -| | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | -| model.SyamOBrien.d1 | | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | -| | fluid-particle drag force. The force is applied to the fluid by | | | -| | multiplying the shear stress by fluid volume fraction. | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | -| | The virtual mass force is not included by default. | | | -| | | | | -| | |VirtualMass_Eq| | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``None`` - virtual mass force calculation is skipped | | | -| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | -| | * ``Constant`` - specified constant coefficient | | | -| | * ``Zuber`` - :cite:p:`zuber64` | | | -| | | | | -| | |VM_Zuber_Eq| | | | -| | | | | -| | * ``Nijssen`` - :cite:p:`nijssen20` | | | -| | | | | -| | |VM_Nijssen_Eq| | | | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ -| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | -+--------------------------+-------------------------------------------------------------------------------+-------------+-----------+ ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| | Description | Type | Default | ++==========================+=========================================================================================+=============+===========+ +| model | Fluid-particle drag model. Options: | String | None | +| | | | | +| | * ``WenYu``- Wen-Yu :cite:p:`Wen66` | | | +| | * ``Gidaspow`` :cite:p:`ding90,Lat00` | | | +| | * ``BVK2`` - Beetstra, van der Hoef, and Kuipers :cite:p:`Bee07,Bee071,tenneti11,tang15`| | | +| | * ``SyamOBrien`` - Syamlal-O'Brien :cite:p:`syam88` | | | +| | * ``UserDrag`` - User-supplied drag model | | | +| | | | | +| | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| model.SyamOBrien.c1 | Fitting parameters for ``SyamOBrien`` drag model | Real | None | +| model.SyamOBrien.d1 | | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| include_divtau | Interpolate the fluid shear stress to particles and include in the | Int | 0 | +| | fluid-particle drag force. The force is applied to the fluid by | | | +| | multiplying the shear stress by fluid volume fraction. | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass | Include virtual mass force in fluid-particle momentum transfer. | String | None | +| | The virtual mass force is not included by default. | | | +| | | | | +| | |VirtualMass_Eq| | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - virtual mass force calculation is skipped | | | +| | * ``Null`` - zero virtual mass coefficient, :math:`C_p,vm = 0` (for testing) | | | +| | * ``Constant`` - specified constant coefficient | | | +| | * ``Zuber`` - :cite:p:`zuber64` | | | +| | | | | +| | |VM_Zuber_Eq| | | | +| | | | | +| | * ``Nijssen`` - :cite:p:`nijssen20` | | | +| | | | | +| | |VM_Nijssen_Eq| | | | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ +| virtual_mass.constant | Constant virtual-mass coefficient. | Real | 0.5 | ++--------------------------+-----------------------------------------------------------------------------------------+-------------+-----------+ .. note:: diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index 1436464..9ad6278 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -9,71 +9,113 @@ Solids settings 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| molecular_weight | Value of constant solid molecular weight. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| density | Specify which density model to use for solids. | String | None | -| | | | | -| | Options (case-insensive): | | | -| | | | | -| | * ``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. | | | -| | | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat | Specify which specific heat model to use for solids. | String | constant | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant specific heat | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| specific_heat.constant | Constant solids specific heat value. | Real | 0 | -| | A value is required for ``constant`` specific heat model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| species | Specify which species can constitute | Strings | None | -| | the solids phase. Defined species must | | | -| | be a subset of the species defined in ``species.solve`` | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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 | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | -| | Options: | | | -| | | | | -| | * ``constant`` - constant thermal conductivity | | | -| | | | | -| | A thermal conductivity model is not required for solids. | | | -| | Defining a thermal conductivity model enables conductive | | | -| | heat transfer for DEM. However, it has no effect for PIC | | | -| | solids heat transfer. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | -| | A value is required for ``constant`` thermal conductivity | | | -| | model. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ -| 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. | | | -+------------------------------------------+-------------------------------------------------------------+----------+----------+ ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| | 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| molecular_weight | Value of constant solid molecular weight. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| density | Specify which density model to use for solids. | String | None | +| | | | | +| | Options (case-insensive): | | | +| | | | | +| | * ``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. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat | Specify which specific heat model to use for solids. | String | constant | +| | | | | +| | 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}` | | | +| | | | | +| | A specific heat model is required if advecting enthalpy. | | | +| | Furthermore, the model must be ``mixture`` if species | | | +| | equations are solved. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.constant | Constant solids specific heat value. | Real | 0 | +| | A value is required for ``constant`` specific heat model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.ignore_discontinuities | Ignore discontinuities in NASA-7 polynomials at transition | Int | 0 | +| | temperatures (``Tsplit``). More information is provied in | | | +| | :ref:`species specific heats`. | | | +| | | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.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`. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 | Reals | 1000 | +| | polynomials. | | | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial | | | +| | defines the specific heat and enthalpy for all | | | +| | temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` | | | +| | polynomials are used to define the specific heat and | | | +| | enthalpy. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| reference_temperature | Reference temperature used for specific enthalpy. | Real | 0 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| species | Specify which species can constitute | Strings | None | +| | the solids phase. Defined species must | | | +| | be a subset of the species defined in ``species.solve`` | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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 | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity | Specify which thermal conductivity model to use for solids. | String | constant | +| | Options: | | | +| | | | | +| | * ``constant`` - constant thermal conductivity | | | +| | | | | +| | A thermal conductivity model is not required for solids. | | | +| | Defining a thermal conductivity model enables conductive | | | +| | heat transfer for DEM. However, it has no effect for PIC | | | +| | solids heat transfer. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| thermal_conductivity.constant | Constant solids thermal conductivity value. | Real | 0 | +| | A value is required for ``constant`` thermal conductivity | | | +| | model. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ +| 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. | | | ++------------------------------------------+--------------------------------------------------------------+----------+----------+ Below is an example for specifying the solids solver model options. diff --git a/docs/source_docs/user_guide/inputs/species.rst b/docs/source_docs/user_guide/inputs/species.rst index 17037ad..f02aa4e 100644 --- a/docs/source_docs/user_guide/inputs/species.rst +++ b/docs/source_docs/user_guide/inputs/species.rst @@ -1,149 +1,163 @@ .. _InputsSpeciesDefs: + +Species solver +============== + +Enabling the species solver +--------------------------- + +The following input is defined using the prefix ``mfix``: + ++-------------------------+------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=========================+========================================================================+==========+===========+ +| solve_species | Enable the MFIX species solver. | Int | 0 | ++-------------------------+------------------------------------------------------------------------+----------+-----------+ + Species definitions -=================== +------------------- The following inputs are defined using the prefix ``species``: -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+==========================+========================================================================+==========+===========+ -| solve | Specified name(s) of the species or None to disable the | Strings | None | -| | species solver. The name assigned to the species | | | -| | solver is used to specify species inputs. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity | Fluid species diffusivity model. | String | None | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - a constant diffusion coefficient is used for | | | -| | all fluid species. | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | -| | ``constant`` diffusivity model | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat | Species specific heat model. This setting only applies if either | String | | -| | ``fluid.specific_heat = mixture`` or | | | -| | ``solids.specific_heat = mixture`` | | | -| | | | | -| | Options (case-insensitive): | | | -| | | | | -| | * ``constant`` - a constant specific heat is defined for each species | | | -| | and a *mixture* specific heat is computed based on the fluid or | | | -| | particle composition. | | | -| | | | | -| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | -| | or more polynomials that are a function of temperature, then the | | | -| | *mixture* specific heat is computed based on the fluid or particle | | | -| | composition. | | | -| | | | | -| | NASA7 polynomial format: | | | -| | | | | -| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | -| | | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| ignore_discontinuities | MFIX-Exa asserts that the low and high NASA-7 polynomials are | Int | 0 | -| | continuous at the 1000K temperature split. However, some polynomials | | | -| | do not have valid high temperature coefficients, for example, liquid | | | -| | water. By setting this keyword to 1, the simulation will only warn | | | -| | that the polynomials are discontinuous. Simulations that use | | | -| | discontinuous polynomials and cross between the 1000K temperature | | | -| | split may fail when computing temperature from enthalpy. | | | -| | Therefore, this option should only be used when the simulation is not | | | -| | expected to cross over the transition temperature. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | | -| | ``fluid.viscosity.molecular = mixture``. | | | -| | | | | -| | Options: | | | -| | | | | -| | * ``constant`` - constant viscosity | | | -| | * ``Sutherland`` :cite:p:`suth1893` | | | -| | * ``Reid`` :cite:p:`reid87` | | | -| | | | | -| | See :ref:`viscosity model descriptions ` | | | -| | in fluid section. | | | -+--------------------------+------------------------------------------------------------------------+----------+-----------+ ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++======================================+=========================================================================+==========+===========+ +| solve | Specified name(s) of the species or None to disable the | Strings | None | +| | species solver. The name assigned to the species | | | +| | solver is used to specify species inputs. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity | Fluid species diffusivity model. | String | None | +| | | | | +| | Options: | | | +| | | | | +| | * ``constant`` - a constant diffusion coefficient is used for | | | +| | all fluid species. | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| diffusivity.constant | Constant species diffusivity - required for | Real | 0 | +| | ``constant`` diffusivity model | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat | Species specific heat model. This setting only applies if either | String | | +| | ``fluid.specific_heat = mixture`` or | | | +| | ``solids.specific_heat = mixture`` | | | +| | | | | +| | Options (case-insensitive): | | | +| | | | | +| | * ``constant`` - a constant specific heat is defined for each species | | | +| | and a *mixture* specific heat is computed based on the fluid or | | | +| | particle composition. | | | +| | | | | +| | * ``NASA7-poly`` - the specific heat of each species is defined by one | | | +| | or more polynomials that are a function of temperature, then the | | | +| | *mixture* specific heat is computed based on the fluid or particle | | | +| | composition. | | | +| | | | | +| | NASA7 polynomial format: | | | +| | | | | +| | :math:`c_p(T)/R = \sum_{i=0}^4 a_{i}T^{i}` | | | +| | | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.ignore_discontinuities | MFIX-Exa checks that NASA-7 polynomials are continuous at the specified | Int | 0 | +| | transition temperature(s)(``Tsplit``). However, some polynomials do not | | | +| | have valid high temperature coefficients, for example, liquid water. By | | | +| | setting this keyword to 1, the simulation will only warn that the | | | +| | polynomials are discontinuous. Simulations that use discontinuous | | | +| | polynomials and cross ``Tsplit`` boundaries may fail when computing | | | +| | temperature from enthalpy. Therefore, this option should only be used | | | +| | when the simulation is not expected to cross over the transition | | | +| | temperature. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular | Molecular viscosity model of species. This setting only applies when | String | None | +| | ``[fluid_name].viscosity.molecular`` = ``mixture``. | | | +| | | | | +| | Options: | | | +| | | | | +| | * ``None`` - no viscosity model | | | +| | * ``constant`` - constant viscosity | | | +| | * ``Sutherland`` :cite:p:`suth1893` | | | +| | * ``Reid`` :cite:p:`reid87` | | | +| | | | | +| | See :ref:`viscosity model descriptions ` | | | +| | in fluid section. | | | ++--------------------------------------+-------------------------------------------------------------------------+----------+-----------+ The following inputs are defined for each species using the prefix ``species.[species_name]``: -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| | Description | Type | Default | -+=======================================+========================================================================+==========+===========+ -| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | -| | model and when including chemical reactions. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | -| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | -| | specific heat model is ``constant`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Real | None | -| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | -| | and the specific heat model is ``NASA7-poly`` | | | -| | | | | -| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | -| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | -| | heat while all six coefficients are used to compute enthalpy. | | | -| | The seventh coefficient, ``a6``, may be provided for completeness | | | -| | but is not used as MFIX-Exa presently does not need to compute | | | -| | entropy. | | | -| | | | | -| | * By default, two sets of coefficients are required. The coefficients | | | -| | define two polynomials used to compute the specific heat and | | | -| | enthalpy across low and high temperatures. The split temperature | | | -| | (Tsplit) defines the transition from low-to-high polynomials. | | | -| | | | | -| | * A single polynomial can be used for all temperatures by setting the | | | -| | split temperature to -1. | | | -| | | | | -| | * Three or more polynomials can be used by defining multiple split | | | -| | temperatures. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| specific_heat.NASA7.Tsplit | Defines the temperature transition between NASA-7 polynomials. | Real | 1000. | -| | | | | -| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | -| | the specific heat and enthalpy for all temperatures. | | | -| | | | | -| | * Specifying ``N`` split temperatures where ``N>1``, signifies ``N+1`` | | | -| | polynomials are used to define the specific heat and enthalpy. | | | -| | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | -| | the specific heat model is ``constant`` and chemical reactions | | | -| | are defined. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | -| | | | | -| | A value is required for ``constant`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | -| | | | | -| | A value is required for ``Sutherland`` species viscosity model. | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | -| viscosity.molecular.Reid.B | | | | -| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | -| viscosity.molecular.Reid.D | | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ -| density | Constant species density (only used for solid phase species). | Real | 0. | -| | | | | -| | Needed for all solid species if | | | -| | :ref:`solids density` is set to ``mixture`` or | | | -| | :ref:`density updates` is set to | | | -| | ``Constant_Species_Densities`` | | | -+---------------------------------------+------------------------------------------------------------------------+----------+-----------+ ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| | Description | Type | Default | ++=======================================+=========================================================================+==========+===========+ +| molecular_weight | Molecular weight of species. Required for mixture molecular weight | Real | 0 | +| | model and when including chemical reactions. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.constant | Constant species specific heat. Required for all fluid (solids) | Real | 0 | +| | species if the fluid (solids) ``specific_heat = mixture`` and the | | | +| | specific heat model is ``constant`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.a[i] | Species specific heat polynomial coefficients. Required for all fluid | Reals | None | +| | (solids) species if the fluid (solids) ``specific_heat = mixture`` | | | +| | and the specific heat model is ``NASA7-poly`` | | | +| | | | | +| | * Each polynomial is defined by six coefficients (``a0`` ... ``a5``). | | | +| | Coefficients ``a0`` through ``a4`` are use for computing specific | | | +| | heat while all six coefficients are used to compute enthalpy. | | | +| | The seventh coefficient, ``a6``, may be provided for completeness | | | +| | but is not used as MFIX-Exa presently does not need to compute | | | +| | entropy. | | | +| | | | | +| | * By default, two sets of coefficients are required. The coefficients | | | +| | define two polynomials used to compute the specific heat and | | | +| | enthalpy across low and high temperatures. The transition temperature | | | +| | (``Tsplit``) defines the transition from low to high polynomials. | | | +| | | | | +| | * A single polynomial can be used for all temperatures by setting the | | | +| | transition temperature to -1. | | | +| | | | | +| | * Three or more polynomials can be used by defining multiple transition | | | +| | temperatures. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| specific_heat.NASA7.Tsplit | Defines the transition temperature between NASA-7 polynomials. | Reals | 1000 | +| | | | | +| | * Setting ``Tsplit = -1`` signifies that a single polynomial defines | | | +| | the specific heat and enthalpy for all temperatures. | | | +| | | | | +| | * Specifying ``N`` temperatures signifies that ``N+1`` polynomials are | | | +| | used to define the specific heat and enthalpy. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| enthalpy_of_formation | Enthalpy of formation of species. Required only when | Real | 0 | +| | the specific heat model is ``constant`` and chemical reactions | | | +| | are defined. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.constant | Constant species fluid viscosity. | Real | 0 | +| | | | | +| | A value is required for ``constant`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.T_ref | Sutherland model reference temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.mu_ref | Sutherland model reference viscosity at T_ref for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Sutherland.S | Sutherland model temperature for species. | Real | 0 | +| | | | | +| | A value is required for ``Sutherland`` species viscosity model. | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| viscosity.molecular.Reid.A | Reid model constants for species. | Real | 0 | +| viscosity.molecular.Reid.B | | | | +| viscosity.molecular.Reid.C | Values are required for ``Reid`` species viscosity model. | | | +| viscosity.molecular.Reid.D | | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ +| density | Constant species density (only used for solid phase species). | Real | 0. | +| | | | | +| | Needed for all solid species if | | | +| | :ref:`solids density` is set to ``mixture`` or | | | +| | :ref:`density updates` is set to | | | +| | ``Constant_Species_Densities`` | | | ++---------------------------------------+-------------------------------------------------------------------------+----------+-----------+ Example inputs diff --git a/duplicate-keys.txt b/duplicate-keys.txt index efa155e..dd14b43 100644 --- a/duplicate-keys.txt +++ b/duplicate-keys.txt @@ -1,11 +1,4 @@ Duplicate keys: -ascent.actions: ascent.rst, plotting.rst -bc.[region_name].[fluid_name].massflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].velocity: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volflow: boundary_conditions.rst, boundary_conditions.rst -bc.[region_name].[fluid_name].volfrac: boundary_conditions.rst, boundary_conditions.rst hypre.bamg_num_down_sweeps: multigrid_solvers.rst, multigrid_solvers.rst hypre.bamg_num_sweeps: multigrid_solvers.rst, multigrid_solvers.rst -mfix.ascent_int: ascent.rst, plotting.rst -mfix.ascent_per_approx: ascent.rst, plotting.rst mfix.restart: checkpointing.rst, initialization.rst diff --git a/nonexistent-keys.txt b/nonexistent-keys.txt index 72a209e..de176cc 100644 --- a/nonexistent-keys.txt +++ b/nonexistent-keys.txt @@ -1,4 +1,3 @@ -# This group may be incorrect bc.[region_name].[solid_name].bins bc.[region_name].[solid_name].density.custom bc.[region_name].[solid_name].density.max -- GitLab From b9e71013ea3c8809ffd8cee70ba89a043b4d1d35 Mon Sep 17 00:00:00 2001 From: Charles G Waldman Date: Fri, 19 Dec 2025 14:31:49 -0600 Subject: [PATCH 16/16] table fixup --- docs/source_docs/user_guide/inputs/solids_model.rst | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/solids_model.rst b/docs/source_docs/user_guide/inputs/solids_model.rst index bb68fcc..adfeaa0 100644 --- a/docs/source_docs/user_guide/inputs/solids_model.rst +++ b/docs/source_docs/user_guide/inputs/solids_model.rst @@ -40,14 +40,14 @@ models. The following inputs are defined with the prefix ``solids``: | | | | | | | * ``constant`` - constant solid specific heat | | | | | * ``NASA7-poly`` the solid specific heat is defined by | | | -| | NASA-7 polynomials. | | | +| | NASA-7 polynomials. | | | | | | | | -| | NASA7 polynomial format: | | | +| | 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`. | | | +| | 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` | | | -- GitLab