EOSViscoDiffH2O
Author
Frederic Aitken, Ferdinand Volino
Title
EOSViscoDiffH2O
Description
EOSViscoDiffH2O is a Mathematica package which allows to calculate the Dynamic Viscosity and Self-Diffusion Coefficient for all fluid phases of Water from 200 K to 1800 K and pressures from very low values until 6 GPa
Category
Academic Articles & Supplements
Keywords
equations of state, dynamic viscosity, self-diffusion coefficient, thermodynamic equilibrium, phase transitions, liquids, water, steam, supercooled, IAPWS, correlation length, fractional derivative, Frenkel/Widom line, latent heat of vaporization
URL
http://www.notebookarchive.org/2022-04-9od0wto/
DOI
https://notebookarchive.org/2022-04-9od0wto
Date Added
2022-04-21
Date Last Modified
2022-04-21
File Size
19.21 kilobytes
Supplements
Rights
CC BY 4.0
Download
Open in Wolfram Cloud
EOSViscoDiffH2O
EOSViscoDiffH2O
Authors : F . Aitken and F . Volino
Description
EOSViscoDiffH2O is a Mathematica package which allows to calculate the Dynamic Viscosity and Self-Diffusion Coefficient for all fluid phases of Water from 200 K to 1800 K and pressures from very low values until 6 GPa. The “EOSViscoDiffH2O” package can be combined with the “ThermodynamicData” package included in Mathematica to determine for example the density of water (except for the supercooled phase) as shown in the examples at the bottom of this Notebook.
The equations of state have been published in “A New Single Equation of State to Describe the Dynamic Viscosity and Self-Diffusion Coefficient for all Fluid Phases of Water from 200 K to 1800 K Based on a New Original Microscopic Model,” Physics of Fluids, 33(11), 117112 (2021); https://doi.org/10.1063/5.0069488.
A preprint version of the published paper can also be accessed from the arXiv open-source archive at: https://arxiv.org/abs/2108.10666.
Installation
To download the package from the Notebook Archive open the "[+]" tab and select "Supplements 2022-04-9od0wto_Supplements.zip". It contains the file “EOSViscoDiffH2O. zip”. The installation of the package in Mathematica is described in the ReadMe.txt file included in EOSViscoDiffH2O. zip.
If the installation was successful then you should see a new palette in the “Palettes” menu called “EOSViscoDiffH2O”. This palette allows you to access some equations more quickly.
If the installation was successful then you should see a new palette in the “Palettes” menu called “EOSViscoDiffH2O”. This palette allows you to access some equations more quickly.
Documentation
IIf the installation was successful, you can then access the documentation through the “Help” menu in a Mathematica Notebook: click on “Wolfram Documentation” and then at the bottom of the documentation page click on “Add-ons and Packages”; in the list of “Installed Add-Ons” you should see a cell labeled “EOSViscoDiffH2O”; click on this cell to have access to the documentation of the program with many examples.
The online documentation is also available at the following address :
http://mathematica.g2elab.grenoble-inp.fr/guide/EOSViscoDiffH2O-HTML/html/guide/EOSViscoDiffH2O.html
http://mathematica.g2elab.grenoble-inp.fr/guide/EOSViscoDiffH2O-HTML/html/guide/EOSViscoDiffH2O.html
How to use the package
To be able to use the functions of the package, it is necessary to load it before any new use in a Notebook. To do this, you must write the following command and validate it:
In[]:=
Needs["EOSViscoDiffH2O`"]
You can then calculate any function of the package. For example, the viscosity and self-diffusion coefficient of liquid water at atmospheric pressure for a temperature of 298.15 K (i.e. 25 °C) by choosing or not to use a specific unit system:
(*Densityequationofstateofliquidwateratatmosphericpressurefromthe"ThermodynamicData"package*)ρH20IAPWS[t_,unit_]:=UnitConvert[ThermodynamicData["Water","Density",{"Pressure"Quantity[1,"Atmospheres"],"Temperature"Quantity[t,"Kelvins"]}],unit];
(*Thecalculateddensitiesforatemperatureof298.15Kwithdifferentunits*)ρH20IAPWS[298.15,"Grams"]ρH20IAPWS298.15,"Kilograms"
3
"Centimeters"
3
"Meters"
Out[]=
0.997048
g/
3
cm
Out[]=
997.048
kg/
3
m
(*Thecalculatedviscosityfromthepresent"EOSViscoDiffH2O"packageforatemperatureof298.15Katatmosphericpressure*)etaH2O298.15,QuantityMagnitude@ρH20IAPWS298.15,"Grams"etaH2OQuantity[298.15,"Kelvins"],ρH20IAPWS298.15,"Grams"etaH2OQuantity[25,"DegreesCelsius"],ρH20IAPWS298.15,"Kilograms"
3
"Centimeters"
3
"Centimeters"
3
"Meters"
Out[]=
0.89018
Out[]=
0.89018
cP
Out[]=
0.89018
cP
(*Thecalculatedself-diffusioncoeffcientfromthepresent"EOSViscoDiffH2O"packageforatemperatureof298.15KatatmosphericpressureusingMills'parameters*)DtH2O298.15,QuantityMagnitude@ρH20IAPWS298.15,"Grams",Coefd0.873248,CoefN0.266544DtH2OQuantity[298.15,"Kelvins"],ρH20IAPWS298.15,"Grams",Coefd0.873248,CoefN0.266544DtH2OQuantity[25,"DegreesCelsius"],ρH20IAPWS298.15,"Kilograms",Coefd0.873248,CoefN0.266544
3
"Centimeters"
3
"Centimeters"
3
"Meters"
Out[]=
0.0000231252
Out[]=
0.0000231252/s
2
cm
Out[]=
0.0000231252/s
2
cm
Before using any of the functions, it is strongly recommended to study the corresponding documentation . The documentation can be accessed more quickly by going to the package' s corresponding palette .
Cite this as: Frederic Aitken, Ferdinand Volino, "EOSViscoDiffH2O" from the Notebook Archive (2022), https://notebookarchive.org/2022-04-9od0wto
Download