Changing Polynomial Properties of species from 8 to 16 coefficients

Kenna · October 9, 2013, 12:25

P { margin-bottom: 0.21cm; } I have a case in which I need the implement thermophysical properties for a wide temperature range of fluids from 300 K to 30000 K. I got the data but have trouble to let OpenFoam use them.

My first choice for this problem would be copying the existing library about polynomial functions for the thermophysicalProperties-File and make out of 8 coefficients 16. I made a custom solver already which is directed to the copied libary and the solver works with it.

I changed in the make/options file this:

Code:

EXE_INC = \
    -I$(WM_PROJECT_DIR)/src/thermophysicalModels/basic/lnInclude \

by this

Code:

EXE_INC = \
    -I$(WM_PROJECT_USER_DIR)/src/thermophysicalModels/basic/lnInclude \

Is it enough to change in all relevant files the 8 to a 16? Relevant files might be for example:

thermophysicalModels/specie/include/thermoPhysicsTyes.H //i.e icoPolynomial<specie, 16>,

For icoPolynomial:
src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomial.C
src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomial.H //i.e. PolySize=16
src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomialI.H

And the same files for:
thermo/hPolynomial
transport/polynomial

I tried to change only the icoPolynomial part. Used ./Allwmake on /src/thermophysicalModels
and a got a Error message from the terminal.

Code:

 kennametal@kennametal-Precision-WorkStation-T3500:~/OpenFOAM/kennametal-2.2.0/src/thermophysicalModels$ ./Allwmake 
 + wmake libso specie 
 Making dependency list for source file reaction/reactions/makeLangmuirHinshelwoodReactions.C 
 Making dependency list for source file reaction/reactions/makeReactions.C 
 SOURCE=reaction/reactions/makeReactions.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100  -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/makeReactions.o 
 SOURCE=reaction/reactions/makeLangmuirHinshelwoodReactions.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100  -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/makeLangmuirHinshelwoodReactions.o 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libspecie.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libspecie.so] Fehler 1 
 + wmake libso solidSpecie 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidSpecie.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidSpecie.so] Fehler 1 
 + wmake libso thermophysicalFunctions 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libthermophysicalFunctions.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libthermophysicalFunctions.so] Fehler 1 
 + ./properties/Allwmake 
 + wmake libso liquidProperties 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libliquidProperties.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libliquidProperties.so] Fehler 1 
 + wmake libso liquidMixtureProperties 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libliquidMixtureProperties.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libliquidMixtureProperties.so] Fehler 1 
 + wmake libso solidProperties 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidProperties.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidProperties.so] Fehler 1 
 + wmake libso solidMixtureProperties 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidMixtureProperties.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidMixtureProperties.so] Fehler 1 
 + wmake libso basic 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so] Fehler 1 
 + wmake libso reactionThermo 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so] Fehler 1 
 + wmake libso laminarFlameSpeed 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/liblaminarFlameSpeedModels.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/liblaminarFlameSpeedModels.so] Fehler 1 
 + wmake libso chemistryModel 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libchemistryModel.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libchemistryModel.so] Fehler 1 
 + wmake libso barotropicCompressibilityModel 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libbarotropicCompressibilityModel.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libbarotropicCompressibilityModel.so] Fehler 1 
 + wmake libso SLGThermo 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libSLGThermo.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libSLGThermo.so] Fehler 1 
 + wmake libso solidThermo 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidThermo.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidThermo.so] Fehler 1 
 + wmake libso solidChemistryModel 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidChemistryModel.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libsolidChemistryModel.so] Fehler 1 
 + wmake libso radiationModels 
 /usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libradiationModels.so: Permission denied 
 collect2: ld gab 1 als Ende-Status zurück 
 make: *** [/opt/openfoam220/platforms/linuxGccDPOpt/lib/libradiationModels.so] Fehler 1 
 kennametal@kennametal-Precision-WorkStation-T3500:~/OpenFOAM/kennametal-2.2.0/src/thermophysicalModels$

In short I think the first problem starts here
/usr/bin/ld: cannot open output file /opt/openfoam220/platforms/linuxGccDPOpt/lib/libspecie.so: Permission denied
I guess I need to adjust somewhere an other option file or an #include command.

I don't know any further. I would appreciate your help.

Thanks
Kenna

Also this post was helpful to me so far:
http://www.cfd-online.com/Forums/ope...olynomial.html

I am using OpenFOAM 2.2.0
and especially the solver compressibleTwoPhasesEulerFoam

Kenna · November 20, 2013, 04:44

So, after some time I tried it again and it works. Maybe the problem was something else in OpenFOAM.

I tested it with the PolynomialTransport.H
I changed the PolySize=8 to PolySize=16

In the thermoPhysicsTypes.H

I changed in the one specific thermoType description with a polynomial Transport the last 8 to a 16.

I tested a case with compressibleTwoPhaseEulerFoam. The Error Messages demandes 16 coefficients. I put zeros at the mueCoeffs till there were 16 numbers. And it calculated.

I also guess it calculates correct. But I didn't check this with other Cases jet.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Natural convection- properties calculation by polynomial eqn	ajayiitkgp	FLUENT	1	June 3, 2014 18:45
Ethanol Gas Polynomial Coefficients	Trustydelux	Siemens	0	March 8, 2012 09:54
Polynomial thermophysical properties and laminar porous Zone	vemps	OpenFOAM	0	June 30, 2011 11:13
species properties of coalchemistryfoam	yingkun	OpenFOAM Running, Solving & CFD	0	May 14, 2011 10:31
polynomial thermophysical properties	jason.ryon	OpenFOAM	2	May 11, 2011 07:16

November 20, 2013, 04:44		#2
Kenna New Member Join Date: Jul 2013 Posts: 14 Rep Power: 13	So, after some time I tried it again and it works. Maybe the problem was something else in OpenFOAM. I tested it with the PolynomialTransport.H I changed the PolySize=8 to PolySize=16 In the thermoPhysicsTypes.H I changed in the one specific thermoType description with a polynomial Transport the last 8 to a 16. I tested a case with compressibleTwoPhaseEulerFoam. The Error Messages demandes 16 coefficients. I put zeros at the mueCoeffs till there were 16 numbers. And it calculated. I also guess it calculates correct. But I didn't check this with other Cases jet.