[muhammadlr95]

Thanks Hans, it worked!

[swsyoon]

Hello

It seems the last post about this model is almost a year ago. I recently tried to use this pemfc model in openfoam4 and was able to compile it according to the earlier post. However when I tried to run the main script by typing ' make run' in run directory, I got this error.
I ran blockMesh command to generate mesh file but I cannot see the 'cellRegionAddressing' file that the error message is referring to.

Can you help to fix this issue ?
Thanks,


( pemfcSinglephaseNonIsothermalSolver | tee log.run );
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.0
Exec : pemfcSinglephaseNonIsothermalSolver
Date : May 19 2020
Time : 08:22:14
Host : "ws-Satellite-L655"
PID : 6511
Case : /home/ws/OpenFOAM/ws-4.0/run/pemfcSinglePhaseModel-4.0/run
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading cell properties

galvanostatic run
V = V [1 2 -3 0 0 -1 0] 0.6
Ibar0 = Ibar0 [0 -2 0 0 0 1 0] 7000
Rhat = Rhat [1 4 -3 0 0 -2 0] 1e-05
Rhat0 = Rhat0 [1 4 -3 0 0 -2 0] 1e-05
CL electronic conductivity(kappaCl) = kappaCl [-1 -2 3 0 0 2 0] 1000
GDL electronic conductivity(kappaGdl) = kappaGdl [-1 -2 3 0 0 2 0] 5000
GDL - BP Contact resistance (rContact)= rContact [1 4 -3 0 0 -2 0] 2e-06
BP electronic conductivity(kappaBp) = kappaBp [-1 -2 3 0 0 2 0] 83000
CL thickness (tCl) = tCl [0 1 0 0 0 0 0] 3.7e-06
GDL thickness (tGdl) = tGdl [0 1 0 0 0 0 0] 0.00041
BP thickness (tBp) = tBp [0 1 0 0 0 0 0] 0.003
concentration overpotential constant = etaConCaConstant [0 0 0 0 0 0 0] 0.25

Tinit = Tinit [0 0 0 1 0 0 0] 353

anodeName = fuel_to_electrolyte
fuelInletName = fuelInlet
fuelOutletName = fuelOutlet
fuelAbpName = fuel_to_abp
cathodeName = air_to_electrolyte
airInletName = airInlet
airOutletName = airOutlet
airCbpName = air_to_cbp
electrolyteAnodeName = electrolyte_to_fuel
electrolyteCathodeName = electrolyte_to_air
abpFuelName = abp_to_fuel
cbpAirName = cbp_to_air

Reading cell field T

Reading cell field k



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/ws/OpenFOAM/ws-4.0/run/pemfcSinglePhaseModel-4.0/run/constant/abp/polyMesh/cellRegionAddressing at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting

[ramakant]

Originally Posted by swsyoon

Hello

It seems the last post about this model is almost a year ago. I recently tried to use this pemfc model in openfoam4 and was able to compile it according to the earlier post. However when I tried to run the main script by typing ' make run' in run directory, I got this error.
I ran blockMesh command to generate mesh file but I cannot see the 'cellRegionAddressing' file that the error message is referring to.

Can you help to fix this issue ?
Thanks,


( pemfcSinglephaseNonIsothermalSolver | tee log.run );
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.0
Exec : pemfcSinglephaseNonIsothermalSolver
Date : May 19 2020
Time : 08:22:14
Host : "ws-Satellite-L655"
PID : 6511
Case : /home/ws/OpenFOAM/ws-4.0/run/pemfcSinglePhaseModel-4.0/run
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading cell properties

galvanostatic run
V = V [1 2 -3 0 0 -1 0] 0.6
Ibar0 = Ibar0 [0 -2 0 0 0 1 0] 7000
Rhat = Rhat [1 4 -3 0 0 -2 0] 1e-05
Rhat0 = Rhat0 [1 4 -3 0 0 -2 0] 1e-05
CL electronic conductivity(kappaCl) = kappaCl [-1 -2 3 0 0 2 0] 1000
GDL electronic conductivity(kappaGdl) = kappaGdl [-1 -2 3 0 0 2 0] 5000
GDL - BP Contact resistance (rContact)= rContact [1 4 -3 0 0 -2 0] 2e-06
BP electronic conductivity(kappaBp) = kappaBp [-1 -2 3 0 0 2 0] 83000
CL thickness (tCl) = tCl [0 1 0 0 0 0 0] 3.7e-06
GDL thickness (tGdl) = tGdl [0 1 0 0 0 0 0] 0.00041
BP thickness (tBp) = tBp [0 1 0 0 0 0 0] 0.003
concentration overpotential constant = etaConCaConstant [0 0 0 0 0 0 0] 0.25

Tinit = Tinit [0 0 0 1 0 0 0] 353

anodeName = fuel_to_electrolyte
fuelInletName = fuelInlet
fuelOutletName = fuelOutlet
fuelAbpName = fuel_to_abp
cathodeName = air_to_electrolyte
airInletName = airInlet
airOutletName = airOutlet
airCbpName = air_to_cbp
electrolyteAnodeName = electrolyte_to_fuel
electrolyteCathodeName = electrolyte_to_air
abpFuelName = abp_to_fuel
cbpAirName = cbp_to_air

Reading cell field T

Reading cell field k



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/ws/OpenFOAM/ws-4.0/run/pemfcSinglePhaseModel-4.0/run/constant/abp/polyMesh/cellRegionAddressing at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting

Yes...I can help you
please follow the steps to run the program and it is working fine for me.
My email id:ramakantgadhewal@gmail.com

[mariduake]

I'm also checking this code. I'm new to OpenFoam, I used the Fluent module, but I wanted to check this open source software. My question is, do you upload a geometry? because I can't find it in the files. I am reading about meshing here https://www.cfd.at/downloads/FoamTut...ExampleTwo.pdf
So, if I want to run this code, d I have to make my own geometry? Sorry I was accustomed with Fluent. Thank you in advance with your response.

[ramakant]

Originally Posted by mariduake

I'm also checking this code. I'm new to OpenFoam, I used the Fluent module, but I wanted to check this open source software. My question is, do you upload a geometry? because I can't find it in the files. I am reading about meshing here https://www.cfd.at/downloads/FoamTut...ExampleTwo.pdf
So, if I want to run this code, d I have to make my own geometry? Sorry I was accustomed with Fluent. Thank you in advance with your response.

Please check it Geometry is already there.

[mat16]

Originally Posted by wyldckat

Greetings to all,

I've split off from the OpenFuelCell thread, the posts that are referring to the package provided here: http://dx.doi.org/10.17632/3gz7pxznzn.1

Apparently that ZIP file was somehow corrupted by someone in the middle... what I mean is that the file "lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H" was replaced by the content from "lib/myPatchToPatchInterpolation/myPatchToPatchInterpolation.H", which is why this problem is occurring. It's possible to recover the broken file, by using the original code "patchToPatchInterpolation.H" from OpenFOAM 4.0 as a basis.

After downloading the ZIP file, I've ran the following commands, in order to be able to compile the code with OpenFOAM 4.0, along with the necessary corrections:

unzip pemfcSinglePhaseModel-4.0.zip
cd pemfcSinglePhaseModel-4.0
cd pemfcSinglePhaseNonIsothermalSolver
find -name "All*" | xargs chmod +x
cp $FOAM_SRC/OpenFOAM/interpolations/patchToPatchInterpolation/patchToPatchInterpolation.H lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H
sed -i -e 's=patchToPatchInterpolation=mypatchToPatchInterpolation=g'  -e 's=PatchToPatchInterpolation=myPatchToPatchInterpolation=g' lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H
./Allwmake > log.make 2>&1
./Allwmake > log.make2 2>&1

The last 2 commands will create 2 log files: "log.make" and "log.make2", where the second file should have something like this inside:

...
..../libSinglephasePEMFC.so' is up to date.
...
.../pemfcSinglephaseNonIsothermalSolver' is up to date.

Best regards,
Bruno

Hey everyone,

I downloaded the version of OpenFoam 4.0 on ubuntu 16.04.

https://openfoamwiki.net/index.php/I...OAM-4.0/Ubuntu

I have followed this procedure but it gives me the same error.





Also during the unzip of the file it asks me:



I do yes and it opens correctly.

Someone help me?

Thanks in advance

[mariduake]

I installed the same and i didn't have any problem compiling. I copied wyldckat instructions line by line and it worked out fine. I got the logs correctly. This is what I did:


1. After instalation I ran mkdir -p $FOAM_RUN
2. I copied the zip to the user4.0/run folder
3. I ran the commands line by line and it compiled just fine.

[mariduake]

Do you have any idea why i have this error. I followed instructions, I even installed OF4.0 instead of 4.1 but I got the same error.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.0
Exec : pemfcSinglephaseNonIsothermalSolver
Date : Dec 02 2020
Time : 05:01:29
Host : "fat-HP-625"
PID : 23372
Case : /home/fat/OpenFOAM/fat-4.0/run/pemfcSinglePhaseModel-4.0/run
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading cell properties

galvanostatic run
V = V [1 2 -3 0 0 -1 0] 0.6
Ibar0 = Ibar0 [0 -2 0 0 0 1 0] 7000
Rhat = Rhat [1 4 -3 0 0 -2 0] 1e-05
Rhat0 = Rhat0 [1 4 -3 0 0 -2 0] 1e-05
CL electronic conductivity(kappaCl) = kappaCl [-1 -2 3 0 0 2 0] 1000
GDL electronic conductivity(kappaGdl) = kappaGdl [-1 -2 3 0 0 2 0] 5000
GDL - BP Contact resistance (rContact)= rContact [1 4 -3 0 0 -2 0] 2e-06
BP electronic conductivity(kappaBp) = kappaBp [-1 -2 3 0 0 2 0] 83000
CL thickness (tCl) = tCl [0 1 0 0 0 0 0] 3.7e-06
GDL thickness (tGdl) = tGdl [0 1 0 0 0 0 0] 0.00041
BP thickness (tBp) = tBp [0 1 0 0 0 0 0] 0.003
concentration overpotential constant = etaConCaConstant [0 0 0 0 0 0 0] 0.25

Tinit = Tinit [0 0 0 1 0 0 0] 353

anodeName = fuel_to_electrolyte
fuelInletName = fuelInlet
fuelOutletName = fuelOutlet
fuelAbpName = fuel_to_abp
cathodeName = air_to_electrolyte
airInletName = airInlet
airOutletName = airOutlet
airCbpName = air_to_cbp
electrolyteAnodeName = electrolyte_to_fuel
electrolyteCathodeName = electrolyte_to_air
abpFuelName = abp_to_fuel
cbpAirName = cbp_to_air

Reading cell field T

Reading cell field k



--> FOAM FATAL ERROR:
Cannot find file "points" in directory "abp/polyMesh" in times 0 down to constant

From function Foam::word Foam::Time::findInstance(const Foam::fileName&, const Foam::word&, Foam::IOobject::readOption, const Foam::word&) const
in file db/Time/findInstance.C at line 182.

FOAM exiting

[mat16]

Hi Maricarmen,

It worked!

Many thanks for your support!

Best regards

[mariduake]

Has anyone tried to modify the fuel? What happens if I change H2 for CH3OH? I think for Methanol model the cathode remains the same, but the anode equations are modified. Has anyone tried to do so?

[ramakant]

Originally Posted by mariduake

Has anyone tried to modify the fuel? What happens if I change H2 for CH3OH? I think for Methanol model the cathode remains the same, but the anode equations are modified. Has anyone tried to do so?

Are you writing code for DMFC in OpenFOAM? or trying to modify the PEM fuel cell code?

[mariduake]

Originally Posted by ramakant

Are you writing code for DMFC in OpenFOAM? or trying to modify the PEM fuel cell code?

I'm trying to modify it. I know there are other phenomena but I just want to check how this can be done. I saw someone modified the Fluent pem fuel cell module but it didn't say how they did it.

[ramakant]

Originally Posted by mariduake

I'm trying to modify it. I know there are other phenomena but I just want to check how this can be done. I saw someone modified the Fluent pem fuel cell module but it didn't say how they did it.

Could you please share the Fluent file?

[mat16]

Originally Posted by mariduake

Do you have any idea why i have this error. I followed instructions, I even installed OF4.0 instead of 4.1 but I got the same error.

Hi,
I changed the geometry and I have the same error as you.




Can someone help me?

Thank you.

[jose.ceballos]

Originally Posted by wyldckat

Greetings to all,

I've split off from the OpenFuelCell thread, the posts that are referring to the package provided here: http://dx.doi.org/10.17632/3gz7pxznzn.1

Apparently that ZIP file was somehow corrupted by someone in the middle... what I mean is that the file "lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H" was replaced by the content from "lib/myPatchToPatchInterpolation/myPatchToPatchInterpolation.H", which is why this problem is occurring. It's possible to recover the broken file, by using the original code "patchToPatchInterpolation.H" from OpenFOAM 4.0 as a basis.

After downloading the ZIP file, I've ran the following commands, in order to be able to compile the code with OpenFOAM 4.0, along with the necessary corrections:

unzip pemfcSinglePhaseModel-4.0.zip
cd pemfcSinglePhaseModel-4.0
cd pemfcSinglePhaseNonIsothermalSolver
find -name "All*" | xargs chmod +x
cp $FOAM_SRC/OpenFOAM/interpolations/patchToPatchInterpolation/patchToPatchInterpolation.H lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H
sed -i -e 's=patchToPatchInterpolation=mypatchToPatchInterpolation=g'  -e 's=PatchToPatchInterpolation=myPatchToPatchInterpolation=g' lib/myPatchToPatchInterpolation/mypatchToPatchInterpolation.H
./Allwmake > log.make 2>&1
./Allwmake > log.make2 2>&1

The last 2 commands will create 2 log files: "log.make" and "log.make2", where the second file should have something like this inside:

...
..../libSinglephasePEMFC.so' is up to date.
...
.../pemfcSinglephaseNonIsothermalSolver' is up to date.

Best regards,
Bruno

I've followed the instructions, and it seems the previous errror is fixed. But I can't obtaint the up date message, instead I got this message:

+ cd lib
+ ./Allwmake
+ wmake libso
'/home/cicy/OpenFOAM/cicy-4.x/platforms/linux64GccDPInt32Opt/lib/libSinglephasePEMFC.so' is up to date.
+ pwd
/home/cicy/OpenFOAM/cicy-4.x/run/pemfcSinglePhaseModel-4.0/pemfcSinglePhaseNonIsothermalSolver
+ cd applications
+ wmake
g++ -std=c++0x -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3  -DNoRepository -ftemplate-depth-100 -g -I../src/constants -I../src/createFields -I../src/createMesh -I../src/createSpecies -I../src/diffusivity -I../src/electrochemistry -I../src/energyTransport -I../src/mapping -I../src/momentumTransport -I../src/multiSpeciesTransport -I../src/patch -I../src/physicalProperties -I../src/readParameters -I../src/readProperties -I../src/tools/functions -I../src/tools/functionIntegration -I../src/tools/ReynoldsNumber -I../src/tools/RiddersRoot -I../src/tools/varInit -I../lib/lnInclude -I/home/cicy/OpenFOAM/OpenFOAM-4.x/src/finiteVolume/lnInclude -I/home/cicy/OpenFOAM/OpenFOAM-4.x/src/meshTools/lnInclude -I/home/cicy/OpenFOAM/OpenFOAM-4.x/src/thermophysicalModels/specie/lnInclude  -IlnInclude -I. -I/home/cicy/OpenFOAM/OpenFOAM-4.x/src/OpenFOAM/lnInclude -I/home/cicy/OpenFOAM/OpenFOAM-4.x/src/OSspecific/POSIX/lnInclude   -fPIC -c pemfcSinglephaseNonIsothermalSolver.C -o Make/linux64GccDPInt32Opt/pemfcSinglephaseNonIsothermalSolver.o
In file included from pemfcSinglephaseNonIsothermalSolver.C:104:0:
../src/patch/createPatchToPatchInterpolation.H: In function ‘int main(int, char**)’:
../src/patch/createPatchToPatchInterpolation.H:10:33: error: missing template arguments before ‘anodeToCathode’
     myPatchToPatchInterpolation anodeToCathode
                                 ^~~~~~~~~~~~~~
../src/patch/createPatchToPatchInterpolation.H:16:33: error: missing template arguments before ‘cathodeToAnode’
     myPatchToPatchInterpolation cathodeToAnode
                                 ^~~~~~~~~~~~~~
../src/patch/createPatchToPatchInterpolation.H:58:33: error: missing template arguments before ‘fuelAnodeToElectrolyteAnode’
     myPatchToPatchInterpolation fuelAnodeToElectrolyteAnode
                                 ^~~~~~~~~~~~~~~~~~~~~~~~~~~
../src/patch/createPatchToPatchInterpolation.H:64:33: error: missing template arguments before ‘airCathodeToElectrolyteCathode’
     myPatchToPatchInterpolation airCathodeToElectrolyteCathode
                                 ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
In file included from pemfcSinglephaseNonIsothermalSolver.C:141:0:
../src/electrochemistry/solveElectrochemistry.H:75:18: error: ‘anodeToCathode’ was not declared in this scope
                  anodeToCathode.faceInterpolate
                  ^~~~~~~~~~~~~~
../src/electrochemistry/solveElectrochemistry.H:75:18: note: suggested alternative: ‘anodePatch’
                  anodeToCathode.faceInterpolate
                  ^~~~~~~~~~~~~~
                  anodePatch
In file included from ../src/electrochemistry/ASRfunction.H:8:0,
                 from ../src/electrochemistry/solveElectrochemistry.H:158,
                 from pemfcSinglephaseNonIsothermalSolver.C:141:
../src/electrochemistry/sigmaElectrolyte.H:78:17: error: ‘anodeToCathode’ was not declared in this scope
     wActivity = anodeToCathode.faceInterpolate(wActivityAnode) + wActivityCathode;
                 ^~~~~~~~~~~~~~
../src/electrochemistry/sigmaElectrolyte.H:78:17: note: suggested alternative: ‘pSatCathode’
     wActivity = anodeToCathode.faceInterpolate(wActivityAnode) + wActivityCathode;
                 ^~~~~~~~~~~~~~
                 pSatCathode
In file included from pemfcSinglephaseNonIsothermalSolver.C:141:0:
../src/electrochemistry/solveElectrochemistry.H:251:18: error: ‘cathodeToAnode’ was not declared in this scope
                  cathodeToAnode.faceInterpolate(mflux[s])*(1.0 - YsBC);
                  ^~~~~~~~~~~~~~
../src/electrochemistry/solveElectrochemistry.H:251:18: note: suggested alternative: ‘cathodeZones’
                  cathodeToAnode.faceInterpolate(mflux[s])*(1.0 - YsBC);
                  ^~~~~~~~~~~~~~
                  cathodeZones
../src/electrochemistry/solveElectrochemistry.H:274:11: error: ‘cathodeToAnode’ was not declared in this scope
          -cathodeToAnode.faceInterpolate(mfluxSum)
           ^~~~~~~~~~~~~~
../src/electrochemistry/solveElectrochemistry.H:274:11: note: suggested alternative: ‘cathodeZones’
          -cathodeToAnode.faceInterpolate(mfluxSum)
           ^~~~~~~~~~~~~~
           cathodeZones
/home/cicy/OpenFOAM/OpenFOAM-4.x/wmake/rules/General/transform:8: recipe for target 'Make/linux64GccDPInt32Opt/pemfcSinglephaseNonIsothermalSolver.o' failed
make: *** [Make/linux64GccDPInt32Opt/pemfcSinglephaseNonIsothermalSolver.o] Error 1

I've now installed OF4.X in Ubuntu 18.04 Windows Bash....

[jose.ceballos]

Someone could upload the pemfcSinglePhaseModel-4.0 fixed to downloand it?

[Pratik_21]

Did anyone try running this code in parallel ?