From e10e48d20d9bce4309ba36e41f1f182541706d02 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Tue, 27 Aug 2024 11:52:00 -0400 Subject: [PATCH 1/6] Add doc for suspension viscosity models --- docs/source_docs/refs.bib | 53 +++++++++++++++++++ .../user_guide/inputs/fluid_model.rst | 36 ++++++++++++- 2 files changed, 88 insertions(+), 1 deletion(-) diff --git a/docs/source_docs/refs.bib b/docs/source_docs/refs.bib index 0693753..92e5c82 100644 --- a/docs/source_docs/refs.bib +++ b/docs/source_docs/refs.bib @@ -82,6 +82,15 @@ pages = {379-402} url = {https://www.sciencedirect.com/science/article/pii/S0021999120305945}, } +@article{brinkman52, + author = {Brinkman, H.C.}, + title = {The viscosity of concentrated suspensions and solutions}, + journal = {The Journal of Chemical Physics}, + year = {1952}, + volume = {20}, + pages = {571} +} + @book{Burden10, title={Numerical Analysis}, author={Burden, R.L. and Faires, J.D.}, @@ -101,6 +110,14 @@ pages = {379-402} pages = {1-31} } +@article{cheng03, + author = {Nian-Sheng Cheng and Adrian Wing-Keung Law}, + title = {Exponential formula for computing effective viscosity}, + journal = {Powder Technology}, + year = {2003}, + volume = {129}, + pages = {156--160} +} @techreport{chern87, author = {Chern, I.-L., and Colella, P. }, @@ -190,6 +207,15 @@ pages = {379-402} publisher={Oxford University Computing Laboratory Oxford, UK} } +@article{einstein11, + title={Eine neue bestimmung der moleküldimensionen}, + author={Einstein, A}, + journal={Ann. Physik}, + volume={19}, + pages={289--306}, + year={1911}, +} + @article{elhilo12, author = {M. El-Hilo and R.W. Chantrell}, title = {Rationalisation of distribution functions for models of nanoparticle magnetism}, @@ -238,6 +264,15 @@ url = {https://www.sciencedirect.com/science/article/pii/S0304885312002132}, publisher={Elsevier} } +@article{gibilaro07, + author = {Gibilaro L.G., Gallucci K., Felice R. Di. and Pagliai P.}, + title = {On the apparent viscosity of a fluidized bed}, + journal = {Chemical Engineering Science}, + year = {2007}, + volume = {62}, + pages = {294--300} +} + @article{Gid83, author = {Gidaspow, D. and Ettehadieh, B. }, title = {Fluidization in Two-Dimensional Beds with a Jet; 2. Hydrodynamic Modeling}, @@ -290,6 +325,15 @@ doi = {https://doi.org/10.1016/j.jcp.2022.111305}, pages = {843-855} } +@article{krieger59, + author = {Krieger, I.M. and Dougherty, T.J. }, + title = {A mechanism for non-newtonian flow in suspensions of rigid spheres}, + journal = {Trans. Soc. Rheol}, + year = {1959}, + volume = {3}, + pages = {137--152} +} + @article{Lat011, author = {Lathouwers, D. and Bellan, J. }, title = {Modeling of dense gas-solid reative mixtures applied to biomass pyrolysis in a fluidized bed}, @@ -384,6 +428,15 @@ doi = {https://doi.org/10.1016/j.jcp.2022.111305}, volume = {32} } +@article{roscoe52, + author = {Roscoe, R.}, + title = {The viscosity of suspensions of rigid spheres}, + journal = {British Journal of Applied Physics}, + year = {1952}, + volume = {3}, + pages = {267} +} + @article{Sch33, author = {Schiller, L. and Naumann, A. }, title = {A drag coefficient correlation}, diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 374b671..7f1ef78 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -35,11 +35,45 @@ The following inputs are defined using the ``fluid`` prefix. | 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 ``Samgorinsky-Lilly`` | | | +| | 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 viscosty model. | String | None | +| | | | | +| | Available options include: | | | +| | | | | +| | * ``None`` No eddy viscosity model | | | +| | * ``Einstein`` :cite:p:`einstein11` | | | +| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | | | | +| | :math:`\mu_{susp}=\mu_{mol} (1-\varepsilon_s)^{-n}` | | | +| | | | | +| | * ``Roscoe`` :cite:p:`roscoe52,krieger59` | | | +| | | | | +| | :math:`\mu_{susp}=\mu_{mol} \mu^*`, where | | | +| | :math:`\mu^*=(max(c1-\varepsilon_s,e^{-9})/c1)^{-c2}` | | | +| | | | | +| | * ``ChengLaw`` :cite:p:`cheng03` | | | +| | | | | ++--------------------------------------------+------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | +| | | | | +| | A value is required when using the model ``Brinkman``. | | | ++--------------------------------------------+------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c1 | Constant for max packing in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the model ``Roscoe``. | | | ++--------------------------------------------+------------------------------------------------------------+--------+----------+ +| viscosity.suspension.Roscoe.c2 | Constant for exponent in Roscoe suspension expression. | Real | None | +| | | | | +| | A value is required when using the model ``Roscoe``. | | | ++--------------------------------------------+------------------------------------------------------------+--------+----------+ +| viscosity.suspension.ChengLaw.constant | Constant for exponent in ChengLaw suspension expression. | Real | None | +| | | | | +| | A value is required when using the model ``ChengLaw``. | | | ++--------------------------------------------+------------------------------------------------------------+--------+----------+ | species | Specify the species that constitute the fluid. | String | None | | | | | | | | All listed species must be properly defined. See the | | | -- GitLab From 010db8b93359481042fe590175e8cfbc63fe4863 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Tue, 27 Aug 2024 11:56:53 -0400 Subject: [PATCH 2/6] fix --- docs/source_docs/user_guide/inputs/fluid_model.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 7f1ef78..237dac7 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -53,7 +53,7 @@ The following inputs are defined using the ``fluid`` prefix. | | * ``Roscoe`` :cite:p:`roscoe52,krieger59` | | | | | | | | | | :math:`\mu_{susp}=\mu_{mol} \mu^*`, where | | | -| | :math:`\mu^*=(max(c1-\varepsilon_s,e^{-9})/c1)^{-c2}` | | | +| | :math:`\mu^*=(max(c1-\varepsilon_s,10^{-9})/c1)^{-c2}`| | | | | | | | | | * ``ChengLaw`` :cite:p:`cheng03` | | | | | | | | -- GitLab From 6831df4659a33c99a87817fe2933160813cd99a2 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Tue, 27 Aug 2024 13:58:56 -0400 Subject: [PATCH 3/6] add another reference --- docs/source_docs/refs.bib | 8 ++++++++ docs/source_docs/user_guide/inputs/fluid_model.rst | 2 +- 2 files changed, 9 insertions(+), 1 deletion(-) diff --git a/docs/source_docs/refs.bib b/docs/source_docs/refs.bib index 92e5c82..a250d50 100644 --- a/docs/source_docs/refs.bib +++ b/docs/source_docs/refs.bib @@ -374,6 +374,14 @@ doi = {https://doi.org/10.1016/j.jcp.2022.111305}, publisher = {NOAA/NCAR Boulder Atmospheric Observatory}, } +@article{maron56, + author = {Samuel H Maron and Percy E Pierce}, + title = {Application of ree-eyring generalized flow theory to suspensions of spherical particles}, + journal = {Journal of Colloid Science}, + year = {1956}, + volume = {11}, + pages = {80--95} +} @article{Massoudi03, author = {Mehrdad Massoudi}, diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 237dac7..a8f2a71 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -50,7 +50,7 @@ The following inputs are defined using the ``fluid`` prefix. | | | | | | | :math:`\mu_{susp}=\mu_{mol} (1-\varepsilon_s)^{-n}` | | | | | | | | -| | * ``Roscoe`` :cite:p:`roscoe52,krieger59` | | | +| | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | | | | | | | | :math:`\mu_{susp}=\mu_{mol} \mu^*`, where | | | | | :math:`\mu^*=(max(c1-\varepsilon_s,10^{-9})/c1)^{-c2}`| | | -- GitLab From f6d778129dd5b4e2123ab59eb80e8733cb1def81 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Wed, 28 Aug 2024 10:58:52 -0400 Subject: [PATCH 4/6] update readme --- README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/README.md b/README.md index 585f536..8a510fe 100644 --- a/README.md +++ b/README.md @@ -31,7 +31,7 @@ can be built locally by: 1. Install sphinx: ```shell -> pip install sphinx sphinx_rtd_theme sphinx-tabs +> pip install sphinx sphinx_rtd_theme sphinx-tabs sphinxcontrib-bibtex ``` 2. Build the documentation: -- GitLab From 035d45627b67de1c42e4b7592fa921bb6689d758 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Wed, 28 Aug 2024 14:04:57 -0400 Subject: [PATCH 5/6] apply review changes --- .../user_guide/inputs/fluid_model.rst | 19 +++++++++++-------- 1 file changed, 11 insertions(+), 8 deletions(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index a8f2a71..9d5581a 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -39,23 +39,26 @@ The following inputs are defined using the ``fluid`` prefix. | | eddy viscosity model. | | | +--------------------------------------------+------------------------------------------------------------+--------+----------+ | viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | -+--------------------------------------------+------------------------------------------------------------+--------+----------+ -| viscosity.suspension | Specify suspension viscosty model. | String | None | +--------------------------------------------+------------------------------------------------------------+--------+----------+ +| viscosity.suspension | Specify suspension viscosty model of the form | String | None | +| | :math:`\mu_{susp}=\mu_{mol}(\mu^* - 1)` | | | | | | | | | | Available options include: | | | | | | | | -| | * ``None`` No eddy viscosity model | | | +| | * ``None`` No eddy suspension model | | | +| | :math:`\mu^*=1` | | | +| | | | | | | * ``Einstein`` :cite:p:`einstein11` | | | -| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | :math:`\mu^*=1 + 2.5\varepsilon_s` | | | | | | | | -| | :math:`\mu_{susp}=\mu_{mol} (1-\varepsilon_s)^{-n}` | | | +| | * ``Brinkman`` :cite:p:`brinkman52,gibilaro07` | | | +| | :math:`\mu^*=(1-\varepsilon_s)^{-c}` | | | | | | | | | | * ``Roscoe`` :cite:p:`roscoe52,krieger59,maron56` | | | -| | | | | -| | :math:`\mu_{susp}=\mu_{mol} \mu^*`, where | | | -| | :math:`\mu^*=(max(c1-\varepsilon_s,10^{-9})/c1)^{-c2}`| | | +| | :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}` | | | | | | | | +--------------------------------------------+------------------------------------------------------------+--------+----------+ | viscosity.suspension.Brinkman.constant | Constant for exponent in Brinkman suspension expression. | Real | None | -- GitLab From 7c91212f0e809c684cdc2b433df80b5b6947b2e8 Mon Sep 17 00:00:00 2001 From: Deepak Rangarajan Date: Wed, 28 Aug 2024 14:06:36 -0400 Subject: [PATCH 6/6] fix format --- docs/source_docs/user_guide/inputs/fluid_model.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source_docs/user_guide/inputs/fluid_model.rst b/docs/source_docs/user_guide/inputs/fluid_model.rst index 9d5581a..7e1430d 100644 --- a/docs/source_docs/user_guide/inputs/fluid_model.rst +++ b/docs/source_docs/user_guide/inputs/fluid_model.rst @@ -39,7 +39,7 @@ The following inputs are defined using the ``fluid`` prefix. | | eddy viscosity model. | | | +--------------------------------------------+------------------------------------------------------------+--------+----------+ | viscosity.eddy.WALE.constant | WALE eddy viscosity model constant. | Real | 0.325 | ---------------------------------------------+------------------------------------------------------------+--------+----------+ +---------------------------------------------+------------------------------------------------------------+--------+----------+ | viscosity.suspension | Specify suspension viscosty model of the form | String | None | | | :math:`\mu_{susp}=\mu_{mol}(\mu^* - 1)` | | | | | | | | -- GitLab