CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   ANSYS (https://www.cfd-online.com/Forums/ansys/)
-   -   RGP Files (https://www.cfd-online.com/Forums/ansys/249898-rgp-files.html)

TahminehAd May 15, 2023 08:52

RGP Files
 
Hello everyone, I am using Python to create a RGP file, I am using the following commands that have been suggested by ANSYS to create such file. But I get the error: No module named pyExt.RefProp. I have installed all the pip packages in Python but I still get this error. Can anyone help me to understand how I can solve this problem?

Python commands:

import pyExt.RefProp as arp
fluidsPath = arp.getFluidsPath()
mat = arp.RefPropLib()

fluidList = ['fluids/water.fld']
mat.setup(fluidsPath, fluidList)

interpError = 1.0e-3
Tmin = 300.0
Tmax = 700.0
Pmin = 100000.0
Pmax = 25000000.0

# Generate both the liquid and vapor components, independent settings possible
h2ollow = arp.RGPLiquid("H2OL", Tmin, Tmax, Pmin, Pmax, interpError)
h2ollow.setMetastable(False)
h2ovlow = arp.RGPVapor("H2OV", Tmin, Tmax, Pmin, Pmax, interpError)
h2ovlow.setMetastable(False)
arp.generateRGPfile("H2O-low.rgp", h2ollow, h2ovlow)

# Generate both the liquid and vapor components in the RGP file in one go
h2omed = arp.RGPSettings("H2O", Tmin, Tmax, Pmin, Pmax, 1.0e-4, arp.ADAPT_AUTO_TP)
arp.generateRGPfile("H2O-med.rgp", h2omed)

# Generate just the vapor component in the RGP file, for entire range
h2ovapor = arp.RGPVapor("H2OV", interpError)
arp.generateRGPfile("H2O-vapor.rgp", h2ovapor)

evcelica May 15, 2023 10:07

What ANSYS Version are you using? I have this working for v21.2, and I believe it may be different for different ANSYS versions.
This output a 1399 kB .rgp file names Water.rgp. My version includes no metastable properties.

Paste this into command prompt: Make sure you change the paths and file names in the second part to match where your python script is, saved as a .py file.
My python file name is "genRGPnometa.py" your may be different.
"C:\Program Files\ANSYS Inc\v212\SystemCoupling\bin\systemcoupling" -R "C:\My Programs\RGP_generator\genRGPnometa\genRGPnometa.p y"

Here is what is in my python file, with the stuff I don't need commented out:
Change the last uncommented line to the path where you want your rgp file written.

import pyExt.RefProp as arp
fluidsPath = arp.getFluidsPath()
mat = arp.RefPropLib()

fluidList = ['fluids/Water.fld'] # Change the REFPROP fluid here~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
mat.setup(fluidsPath, fluidList)

interpError = 1e-3 # Change accuracy/file size here, units for T,P are [K] and [Pa]~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Tmin = 300
Tmax = 700
Pmin = 100000.0
Pmax = 25000000.0
fluid_name = 'Water' #This is the Fluid Name, which will be followed by L or V for the Phase~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~
# Component names for Liq, Vap and Mix are NameL, NameV, and NameL

# Generate both the liquid and vapor components in the RGP file in one go

fluid = arp.RGPSettings(fluid_name, Tmin, Tmax, Pmin, Pmax, interpError, arp.ADAPT_AUTO_TP)
print(dir(fluid))
# fluid.setMetastable(False) # For version 202 and before

# Uncomment below for 212 version
fluid_liq = fluid.getLiquid()
fluid_liq.setMetastable(False)
fluid_vap = fluid.getVapor()
fluid_vap.setMetastable(False)

arp.generateRGPfile("C:\My Programs\RGP_generator\genRGPnometa\Water.rgp", fluid) # Provide Output path and filename here~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~

# Generate both the liquid and vapor components, independent settings possible
##h2ollow = arp.RGPLiquid("H2OL", Tmin, Tmax, Pmin, Pmax, interpError)
##h2ollow.setMetastable(False)
##h2ovlow = arp.RGPVapor("H2OV", Tmin, Tmax, Pmin, Pmax, interpError)
##h2ovlow.setMetastable(False)
##arp.generateRGPfile("H2O-low.rgp", h2ollow, h2ovlow)


# Generate just the vapor component in the RGP file, for entire range
##h2ovapor = arp.RGPVapor("H2OV", interpError)
##arp.generateRGPfile("H2O-vapor.rgp", h2ovapor)

TahminehAd May 16, 2023 04:52

Thank you for your response. I am using ANSYS V22.2 and I do exactly what you said. But the issue is that the command "import pyExt.RefProp as arp" gives me the error. Even though I installed all the packages but still it says " No module named 'pyExt' ". I am wondering if I need to install something else or if the folder RefProp needs to be changed? or anything else I should do?

TahminehAd May 16, 2023 05:51

Or if someone could just give me the RGP table for water liquid in normal pressure (ambient pressure 1-1.5Mpa), I would appreciate it a lot.

evcelica May 16, 2023 14:04

For water, you can just use the built in IAWPS library data (If you are using CFX) I don't know if fluent has one?
Otherwise, PM me your email and I'll send you the file I generated using your inputs.

Actually, I just did this with v222, worked perfectly.
Change the path and name file to where the python script is on your machine.
Your path and python name are definitely different than what I have in Bold. Also get rid of the space between p and y, I don't know why it does that here?
"C:\Program Files\ANSYS Inc\v222\SystemCoupling\bin\systemcoupling" -R "C:\My Programs\RGP_generator\genRGPnometa\genRGPnometa.p y"

Make sure you change all the paths to folders on your machine, including where to write the output file.

Bak16 August 1, 2023 12:38

Thank you very much for your help. I tried running the code in Command Console as ANSYS User guide says. Specifically on Command prompt I type
> cd C:\Program Files\ANSYS Inc\v222\SystemCoupling\bin
> systemcoupling.bat
After that I call evcelica's code in a file which I have named RGPGeneration.py with the following command:
exec(open('C:/Users/Documents/RGPGeneration.py').read())

They following are printed on the screen and no .rgp file is generated inside the chosen folder. How can I fix this?

Selecting adapation option for component WaterL
Relative error is 0.001000
Adapting using the automatic option.
Adapting including the saturation boundary.
Number of temperature seed points = 3
Number of pressure seed points = 3
Building tables for component WaterL
Tmin = 3.000000e+02 K
Tmax = 7.000000e+02 K
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Viscosity Tmin = 3.000000e+02 K
Viscosity Tmax = 7.000000e+02 K
Viscosity Pmax = 2.500000e+07 Pa
Conductivity Tmin = 3.000000e+02 K
Conductivity Tmax = 7.000000e+02 K
Conductivity Pmax = 2.500000e+07 Pa
Metastable extrapolation deactivated
Property nT nP Spinodal Failures (7.00e-01) Superheat Failures
Enthalpy 24 81 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Sound speed 21 21 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Volume 73 75 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cv 18 25 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cp 41 45 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
dP/dV_T 36 47 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Entropy 25 85 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Viscosity 29 28 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Conductivity 27 26 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Saturation
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Tmin = 3.727559e+02 K
Tmax = 6.470960e+02 K
Surface tension
Liquid saturation properties.
Vapor saturation properties.

Selecting adapation option for component WaterV
Relative error is 0.001000
Adapting using the automatic option.
Adapting including the saturation boundary.
Number of temperature seed points = 3
Number of pressure seed points = 3
Building tables for component WaterV
Tmin = 3.000000e+02 K
Tmax = 7.000000e+02 K
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Viscosity Tmin = 3.000000e+02 K
Viscosity Tmax = 7.000000e+02 K
Viscosity Pmax = 2.500000e+07 Pa
Conductivity Tmin = 3.000000e+02 K
Conductivity Tmax = 7.000000e+02 K
Conductivity Pmax = 2.500000e+07 Pa
Metastable extrapolation deactivated
Property nT nP Spinodal Failures (6.00e-01) Superheat Failures
Enthalpy 51 127 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Sound speed 19 22 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Volume 187 97 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cv 21 23 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cp 48 52 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
dP/dV_T 90 47 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Entropy 30 84 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Viscosity 17 23 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Conductivity 35 40 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Saturation
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Tmin = 3.727559e+02 K
Tmax = 6.470960e+02 K
Surface tension
Liquid saturation properties.
Vapor saturation properties.

Writing file...
Writing header information...
Writing tables...
Writing WaterL 2D tables...
Writing WaterL saturation tables...
Writing WaterV 2D tables...
Writing WaterV saturation tables...

evcelica August 1, 2023 16:38

Did you change the output path to a valid path on your machine:
It's on this line:
arp.generateRGPfile("C:\My Programs\RGP_generator\genRGPnometa\Water.rgp", fluid) # Provide Output path and filename here~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~


I just open command prompt, change to C:, then paste in:
"C:\Program Files\ANSYS Inc\v222\SystemCoupling\bin\systemcoupling" -R "C:\My Programs\RGP_generator\genRGPnometa\genRGPnometa.p y" with the quotes, and get rid of the space between p and y. Change the path and file names of course in this as well as the .py file.



CMD Outputs:

Microsoft Windows [Version 10.0.19045.3208]
(c) Microsoft Corporation. All rights reserved.

U:\>C:

C:\>"C:\Program Files\ANSYS Inc\v222\SystemCoupling\bin\systemcoupling" -R "C:\My Programs\RGP_generator\genRGPnometa\genRGPnometa.p y"

Ansys(R) System Coupling

Executing from: C:\PROGRA~1\ANSYSI~1\v222\SYSTEM~1\bin\systemcoupl ing.bat

2022 R2

Point Releases and Patches installed:

ANSYS, Inc. License Manager 2022 R2
Autodyn 2022 R2
SpaceClaim 2022 R2
LS-DYNA 2022 R2
CFD-Post only 2022 R2
CFX (includes CFD-Post) 2022 R2
FENSAP-ICE 2022 R2
Fluent (includes CFD-Post) 2022 R2
Polyflow (includes CFD-Post) 2022 R2
Forte (includes EnSight) 2022 R2
TurboGrid 2022 R2
ICEM CFD 2022 R2
Aqwa 2022 R2
Customization Files for User Programmable Features 2022 R2
Mechanical Products 2022 R2
Remote Solve Manager Standalone Services 2022 R2
Viewer 2022 R2
ACIS Geometry Interface 2022 R2
AutoCAD Geometry Interface 2022 R2
Catia, Version 4 Geometry Interface 2022 R2
Catia, Version 5 Geometry Interface 2022 R2
Catia, Version 6 Geometry Interface 2022 R2
Creo Elements/Direct Modeling Geometry Interface 2022 R2
Creo Parametric Geometry Interface 2022 R2
Inventor Geometry Interface 2022 R2
JTOpen Geometry Interface 2022 R2
NX Geometry Interface 2022 R2
Parasolid Geometry Interface 2022 R2
Solid Edge Geometry Interface 2022 R2
SOLIDWORKS Geometry Interface 2022 R2

(c) 2014-2022 ANSYS, Inc. All rights reserved. Unauthorized use, distribution
or duplication is prohibited. This product is subject to U.S. laws governing
export and re-export. For full Legal Notice, see documentation.

executing script 'C:\My Programs\RGP_generator\genRGPnometa\genRGPnometa.p y'
Found REFPROP fluids files at C:\PROGRA~1\ANSYSI~1\v222\SYSTEM~1\runTime\winx64\ etc\RefProp\refprop10
['__class__', '__delattr__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__', '__gt__', '__hash__', '__init__', '__init_subclass__', '__le__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__swig_destroy__', '__weakref__', 'getLiquid', 'getVapor', 'limitsProvided', 'matName', 'this', 'thisown']
Generating RGP file C:\My Programs\RGP_generator\genRGPnometa\Water.rgp
Selecting adapation option for component WaterL
Relative error is 0.001000
Adapting using the automatic option.
Adapting including the saturation boundary.
Number of temperature seed points = 3
Number of pressure seed points = 3
Building tables for component WaterL
Tmin = 3.000000e+02 K
Tmax = 7.000000e+02 K
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Viscosity Tmin = 3.000000e+02 K
Viscosity Tmax = 7.000000e+02 K
Viscosity Pmax = 2.500000e+07 Pa
Conductivity Tmin = 3.000000e+02 K
Conductivity Tmax = 7.000000e+02 K
Conductivity Pmax = 2.500000e+07 Pa
Metastable extrapolation deactivated
Property nT nP Spinodal Failures (7.00e-01) Superheat Failures
Enthalpy 24 81 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Sound speed 21 21 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Volume 73 75 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cv 18 25 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cp 41 45 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
dP/dV_T 36 47 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Entropy 25 85 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Viscosity 29 28 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Conductivity 27 26 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Saturation
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Tmin = 3.727559e+02 K
Tmax = 6.470960e+02 K
Surface tension
Liquid saturation properties.
Vapor saturation properties.

Selecting adapation option for component WaterV
Relative error is 0.001000
Adapting using the automatic option.
Adapting including the saturation boundary.
Number of temperature seed points = 3
Number of pressure seed points = 3
Building tables for component WaterV
Tmin = 3.000000e+02 K
Tmax = 7.000000e+02 K
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Viscosity Tmin = 3.000000e+02 K
Viscosity Tmax = 7.000000e+02 K
Viscosity Pmax = 2.500000e+07 Pa
Conductivity Tmin = 3.000000e+02 K
Conductivity Tmax = 7.000000e+02 K
Conductivity Pmax = 2.500000e+07 Pa
Metastable extrapolation deactivated
Property nT nP Spinodal Failures (6.00e-01) Superheat Failures
Enthalpy 51 127 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Sound speed 19 22 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Volume 187 97 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cv 21 23 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Cp 48 52 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
dP/dV_T 90 47 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Entropy 30 84 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Viscosity 17 23 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Conductivity 35 40 0 0.00% 0 0.00% 0 0.00% | 0 0.00% 0 0.00% 0 0.00%
Saturation
Pmin = 1.000000e+05 Pa
Pmax = 2.500000e+07 Pa
Tmin = 3.727559e+02 K
Tmax = 6.470960e+02 K
Surface tension
Liquid saturation properties.
Vapor saturation properties.

Writing file...
Writing header information...
Writing tables...
Writing WaterL 2D tables...
Writing WaterL saturation tables...
Writing WaterV 2D tables...
Writing WaterV saturation tables...

evcelica August 1, 2023 17:01

Why don't you just make a folder C:/temp and put everything there, including your py file and write your output there as well. Then we can see that your paths are set up correctly.
Right now we can't see what you are doing and where you are writing to.

Bak16 August 2, 2023 10:33

It worked! Thank you so much!


All times are GMT -4. The time now is 13:34.