as the title suggests, I am having some issues running a real fluid simulation with CoolProp (CP from now on). I am using CS-6.1.1. If I remember correctly, CP is now officially supported since CS-6.0 (correct me if I am wrong). Before anything else, I want to point out that I am a fairly new person to CFD and very new to CS.
The error that I get when I try to run with CP with the compressible model is always the same:
Code: Select all
Local case configuration:
Date: mié 23 sep 2020 21:25:37
System: Linux 5.8.7-1-default (Welcome to openSUSE Tumbleweed 20200919 - Kernel)
Machine: localhost.localdomain
Processor: model name : Intel(R) Core(TM) i5-6300HQ CPU @ 2.30GHz
Memory: 15894 MB
User: fernando (Fernando Oleo Blanco)
Directory: /home/fernando/Dirt/cfd/nozzle/co2nozz/RESU/result_compress_correct2
OpenMP threads: 4
Processors/node: 4
Reading metadata from file: "mesh_input.csm"
===============================================================
CALCULATION PREPARATION
=======================
===========================================================
../../../code_saturne-6.1.1/src/base/cs_coolprop.cxx:196: Fatal error.
CoolProp was unable to compute some fluid properties with
input variable names: "H", "P" and backend "HEOS".
First of 1 fluid(s) considered: "CarbonDioxide".
Call stack:
1: 0x7fb5db51bcc9 <cs_phys_prop_coolprop+0x659> (libsaturne-6.1.so)
2: 0x7fb5db43a5a4 <cs_phys_prop_compute+0x2b4> (libsaturne-6.1.so)
3: 0x7fb5db9a22be <csphys_+0x35e> (libsaturne-6.1.so)
4: 0x7fb5db4bd7d4 <iniusi_+0x199> (libsaturne-6.1.so)
5: 0x7fb5db4bd22b <initi1_+0x84> (libsaturne-6.1.so)
6: 0x7fb5dc864ede <main+0x54e> (libcs_solver-6.1.so)
7: 0x7fb5d9f4ecca <__libc_start_main+0xea> (libc.so.6)
8: 0x40107a <_start+0x2a> (cs_solver)
End of stack
The build that I have, had the following configuration (in reality, this is the configuration for CS-6.2.0, but it is the same):
Code: Select all
Configuration options:
use debugging code: no
use malloc hooks: no
use graphical user interface: yes
use long global numbers: yes
Zlib (gzipped file) support: yes
MPI (Message Passing Interface) support: no
OpenMP support: yes
OpenMP Fortran support: yes
CUDA support: no
BLAS (Basic Linear Algebra Subprograms) support: no
ParMETIS (Parallel Graph Partitioning) support: no
METIS (Graph Partitioning) support: no
PT-SCOTCH (Parallel Graph Partitioning) support: no
SCOTCH (Graph Partitioning) support: no
CCM support: no
HDF (Hierarchical Data Format) support: yes
CGNS (CFD General Notation System) support: yes
MED (Model for Exchange of Data) support: yes
MED MPI I/O support: no
MEDCoupling support: no
Catalyst (ParaView co-processing) support: no
Melissa (in-situ statistics) support: no
EOS support: no
freesteam support: no
CoolProp support: yes
PETSc support: no
MUMPS support: no
Dynamic loader support: dlopen
Auto load environment modules: no
PdfLaTeX (document typesetting) available: yes
Doxygen (document generation) available: no
The package has been configured. Type:
make
make install
As it can be seen in the error message above, the simulation fails at the very very beginning.
By taking a look at the cs_coolprop.cxx file, it shows that this message is produced when the number of outputs is not the same as the inputs. While that is in itself fine, I do not think this is my case.
The initialization parameters are all within CP's CO2 range (other fluids also fail, in compressible cases). I have tried different initialization values, and it always fails with the same error. I want to point out that all my tests were generated with the GUI, and that all incompressible simulations run fine when selecting CP and "material laws" for the density and viscosity (however, I have the feeling the values are not correct).
Is this error expected? As in: "this error appears when X is not setup correctly". How could it be prevented?
-------------------------------------
BACKGROUND:
I am currently tasked with designing SCO2 turbomachinery (preliminar design). While doing so, we were pointed to an article (still not published) where it showed that the severe non-linear properties near the critical point could lead to very unintuitive results. I am testing CS to see if it could run real fluid simulations with SCO2 and repeat the results presented in the aforementioned paper.
------------------------------------
FURTHER QUESTIONS:
Is the CP backend supposed to be used with the GUI or does it require manual intervention with some of the user scripts?
Why does CS GUI not allow compressible fluid with varying gamma? This is very important to me since SCO2 does have quite a change near de critical point, and I would like to model it as close as possible.
Is it a bug that when I select compressible fluid model, in the fluid tab, the fluid selection dropdown menu is not shown? And is it expected that the density does not allow to select "material law"? I suppose that in this case, the density is computed from the perfect gas equations and therefore CS does not "allow" for the density to follow a "material law".
Reading the cs_coolprop file, it shows that asking CP for Mach number (as well as some other properties) generates an error. Why is that? I mean, CP does support most (if not all) the properties for which an error is being generated as far as I know.
Could third-party codes be added to the python build script for easy setup? The setup file has no entry for CP, FreeSteam, EOS, etc.
And finally: is CP well tested and used? I have found close to no information about it in the forums nor anywhere else (tutorials, manuals, etc). And the fact that some important properties of fluids are "blacklisted" from the interface, does not look very reassuring.
Kind regards,
Fernando Oleo Blanco