From 00e17023226fbe22bcaded85ae21acfaab4841cc Mon Sep 17 00:00:00 2001
From: fsonnino <31736430+fsonnino@users.noreply.github.com>
Date: Thu, 19 Oct 2017 16:29:45 +0200
Subject: [PATCH] Create CoolProp_Properties_of_Water.c (#1569)

Fluent udf to call the following thermo physical properties of water :
-Density
-Thermal Conductivity
-Viscosity
-Specific Heat

If the user want's to add other properties it should be done in priority with the following entry variables Temperature and Pressure which gives in the udf the generic following statement :
property=PropsSI((char*)"corresponding_property", (char*)"T", temperature, (char*)"P", pressure, (char*)FLUID);
---
 .../Fluent/CoolProp_Properties_of_Water.c     | 83 +++++++++++++++++++
 1 file changed, 83 insertions(+)
 create mode 100644 wrappers/Fluent/CoolProp_Properties_of_Water.c

diff --git a/wrappers/Fluent/CoolProp_Properties_of_Water.c b/wrappers/Fluent/CoolProp_Properties_of_Water.c
new file mode 100644
index 00000000..cf835952
--- /dev/null
+++ b/wrappers/Fluent/CoolProp_Properties_of_Water.c
@@ -0,0 +1,83 @@
+#include "udf.h"
+#include <string.h>
+struct string;
+const char FLUID[] = "Water";
+const real gauge = 101325; /*operating pressure in pascal (as defined in fluent) */
+double PropsSI (char*, char*, double, char*, double, char*);
+
+/*Density of the FLUID[]*/
+DEFINE_PROPERTY(water_density, c, t)
+{
+         real temperature = C_T(c, t);
+         real pressure = C_P(c, t)+gauge;
+         real density;
+         density = PropsSI((char*)"D",(char*)"T", temperature, (char*)"P", pressure, (char*)FLUID);
+         return density;
+}
+/*Thermal Conductivity of the FLUID[]*/
+DEFINE_PROPERTY(water_thermalConductivity, c, t)
+{
+	real thermalConductivity;
+	real temperature = C_T(c, t);
+	real pressure = C_P(c, t)+gauge;
+	thermalConductivity = PropsSI((char*)"L", (char*)"T", temperature, (char*)"P", pressure, (char*)FLUID);
+	return thermalConductivity; 
+}
+/*Viscosity of the FLUID[]*/
+DEFINE_PROPERTY(water_viscosity, c, t)
+{
+	real viscosity;
+	real temperature = C_T(c, t);
+	real pressure = C_P(c, t)+gauge;
+	viscosity = PropsSI((char*)"V", (char*)"T", temperature, (char*)"P", pressure, (char*)FLUID);
+	return viscosity;
+}
+/*Specific heat of the FLUID[]*/
+DEFINE_SPECIFIC_HEAT(water_specificHeat, temperature, Tref, enthalpy, yi)
+{	
+        real pressure ;
+	/*
+        density = 1.7730;
+	*/
+	/* The following commented code is supposed to get the pressure
+	 from the cell to use with Coolprop. Left commented because
+	 specific heat depends very little on pressure. Will increase
+	 computational time significantly. */
+	 
+	
+	 Domain *domain = Get_Domain(1);
+	 Thread *t;
+	 cell_t c;
+	 thread_loop_c(t, domain)
+	 {
+	 	begin_c_loop(c, t)
+	 	{
+	 		real pressure = C_P(c, t)+gauge;
+                        temperature = C_T(c,t);
+                        real specificHeat;
+	                specificHeat = PropsSI((char*)"C",(char*)"T", temperature, (char*)"P", pressure, (char*)FLUID);
+	                *enthalpy = specificHeat*( temperature - Tref );
+			return specificHeat;
+	 	}end_c_loop(c, t)
+	 }
+}
+/* test coolprop integration */
+DEFINE_ON_DEMAND(call_coolprop_water)
+{
+	real p, t, density, specificHeat, thermalConductivity, enthalpy, viscosity;
+
+	p = 101325.0;
+	t = 298.15;
+
+	density = PropsSI((char*)"D",(char*)"T",t,(char*)"P",p,(char*)FLUID);
+	specificHeat = PropsSI((char*)"C",(char*)"T",t,(char*)"P",p,(char*)FLUID);
+	viscosity = PropsSI((char*)"V",(char*)"T",t,(char*)"P",p,(char*)FLUID);
+	thermalConductivity = PropsSI((char*)"L",(char*)"T",t,(char*)"P",p,(char*)FLUID);
+	enthalpy = PropsSI((char*)"H",(char*)"T",t,(char*)"P",p,(char*)FLUID);
+
+	Message("p = %lf, T = %lf => density = %lf\n", p, t, density);
+	Message("p = %lf, T = %lf => specific heat = %lf\n", p, t, specificHeat);
+	Message("p = %lf, T = %lf => viscosity = %lf\n", p, t, viscosity);
+	Message("p = %lf, T = %lf => thermal conductivity = %lf\n", p, t, thermalConductivity);
+	Message("p = %lf, T = %lf => enthalpy = %lf\n", p, t, enthalpy);
+}