Modification of reactingFoam solver to FGMFoam solver

killua1915 · February 6, 2025, 14:07

Hello,

I am trying to follow tutorial by Michael Bertsch to compile FGMFoam from reactingFoam.

https://www.tfd.chalmers.se/~hani/ku...rt_FGMFoam.pdf

I am getting following error and I am not sure what is the issue here.

build_log.zip

Any idea on how to resolve this ?

HTML Code:

In file included from FGMFoam.C:88:PVEqn.H: In function ‘int main(int, char**)’:PVEqn.H:8:20: error: ‘class Foam::CombustionModel<Foam::psiReactionThermo>’ has no member named ‘SourcePV’    8 |          reaction->SourcePV()   //--- added by senbin, 2019-04-26      |                    ^~~~~~~~

HTML Code:

make: *** [/usr/lib/openfoam/openfoam2406/wmake/rules/General/transform:38: Make/linux64GccDPInt32Opt/FGMFoamPost.o] Error 1

The source code of file is given below.
https://drive.google.com/file/d/1ZL3...ew?usp=sharing

dlahaye · February 7, 2025, 05:02

I suggest to use the version of OpenFoam for which the code was developed.

See also post related to FGM on this forum.

killua1915 · February 7, 2025, 09:15

I am not sure in which version it was written. Most files have following comment at top.
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2015 OpenFOAM Foundation
------------------------------------------------------------------------------

dlahaye · February 7, 2025, 11:52

Not sure either. I believe something like OpenFoam-2.2 or 2.3.

The alternative approach is to persue a new implementation. A function object to compute the mixture fraction is already available. Next step its to extend to. computation of the mixture fraction variance.

killua1915 · February 10, 2025, 04:11

I have successfully compiled and run the FGMFoam for SandiaD flame tutorial. But I am getting following warnings.

Quote:

Duplicate entry cubeRootVol in runtime table LESdelta
[stack trace]

I have went through some threads in this forums but I am not sure which library linkage is causing this issue. Any idea on this ?

dlahaye · February 10, 2025, 06:43

Not sure.

Is the entry "cubeRootVol" by mistake repeated in some input file?

killua1915 · February 10, 2025, 07:04

Not really. I have tried using grep function to find and I am getting following results :

Quote:

./TurbulenceModels/turbulenceModels/RAS/kOmegaSSTSAS/kOmegaSSTSAS.H: delta cubeRootVol;
./TurbulenceModels/turbulenceModels/RAS/kOmegaSSTSAS/kOmegaSSTSAS.H: cubeRootVolCoeffs
./TurbulenceModels/turbulenceModels/LES/LESdeltas/vanDriestDelta/vanDriestDelta.H: delta cubeRootVol;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/vanDriestDelta/vanDriestDelta.H: cubeRootVolCoeffs
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: Foam::cubeRootVolDelta
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: cubeRootVolDelta.C
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H:#ifndef cubeRootVolDelta_H
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H:#define cubeRootVolDelta_H
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: Class cubeRootVolDelta Declaration
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H:class cubeRootVolDelta
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: cubeRootVolDelta(const cubeRootVolDelta&) = delete;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: void operator=(const cubeRootVolDelta&) = delete;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: TypeName("cubeRootVol");
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: cubeRootVolDelta
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H: virtual ~cubeRootVolDelta() = default;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C:#include "cubeRootVolDelta.H"
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C: defineTypeNameAndDebug(cubeRootVolDelta, 0);
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C: addToRunTimeSelectionTable(LESdelta, cubeRootVolDelta, dictionary);
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C:void Foam::LESModels::cubeRootVolDelta::calcDelta()
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C:Foam::LESModels::cubeRootVolDel ta::cubeRootVolDelta
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C:void Foam::LESModels::cubeRootVolDelta::read(const dictionary& dict)
./TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C:void Foam::LESModels::cubeRootVolDelta::correct()
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C: cubeRootVolDelta_.calcDelta();
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C: static_cast<const volScalarField&>(cubeRootVolDelta_)
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C: cubeRootVolDelta_
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C: name + "cubeRootVolDelta",
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C: cubeRootVolDelta_.read(dict.subDict(typeName + "Coeffs"));
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.H: Maximum delta between maxDeltaxyz and cubeRootVolDelta.
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.H:#include "cubeRootVolDelta.H"
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.H: //- cubeRootVolDelta model
./TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.H: cubeRootVolDelta cubeRootVolDelta_;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/PrandtlDelta/PrandtlDelta.H: delta cubeRootVol;
./TurbulenceModels/turbulenceModels/LES/LESdeltas/PrandtlDelta/PrandtlDelta.H: cubeRootVolCoeffs
./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/sourceFiles:$(LESdelta)/cubeRootVolDelta/cubeRootVolDelta.C \
./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C.dep:$(OBJECTS_DIR)/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C.dep: \
./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C.dep:LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C \
./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C.dep:LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.H \
grep: ./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.o: binary file matches
grep: ./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.o: binary file matches
./TurbulenceModels/turbulenceModels/Make/linux64GccDPInt32Opt/LES/LESdeltas/maxDeltaxyzCubeRootLESDelta/maxDeltaxyzCubeRootLESDelta.C.dep:lnInclude/cubeRootVolDelta.H \
./TurbulenceModels/turbulenceModels/Make/files:$(LESdelta)/cubeRootVolDelta/cubeRootVolDelta.C

dlahaye · February 10, 2025, 07:34

1/ Suggest to check the case files (system/ and constant/) using grep;

2/ Suggest to share entire log file (allowing to check how and where the error message appears);

killua1915 · February 10, 2025, 09:18

log.zip

Please find attached log files.

dlahaye · February 10, 2025, 09:27

It is possible that errors are due to v2406 being used?

killua1915 · February 10, 2025, 09:36

Not sure. I have tried compiling it in version v2206 as you recommended previously but I was getting same warnings. The solver is running fine but it gives these warnings which I am not sure why.

dlahaye · February 10, 2025, 09:45

I was recommending going back way earlier and compile against even older versions of openfoam. Apologies for not being sufficiently clear.

What are the simulation results like?

killua1915 · February 10, 2025, 10:22

Haven't had the chance to look at it properly but based on initial expectation, contour plots for temperature and OH radicals are quite close as given in this tutorial.

https://www.tfd.chalmers.se/~hani/ku...rt_FGMFoam.pdf

killua1915 · February 12, 2025, 06:20

For anyone wondering how I got rid of these warnings, I renamed all the classes which were mentioned in the warning. Crude approach but it worked. I have tried copying library from src to my local and recompiling it with different name but not sure why it didn't work.

February 6, 2025, 14:07	Modification of reactingFoam solver to FGMFoam solver	#1
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	Hello, I am trying to follow tutorial by Michael Bertsch to compile FGMFoam from reactingFoam. https://www.tfd.chalmers.se/~hani/ku...rt_FGMFoam.pdf I am getting following error and I am not sure what is the issue here. build_log.zip Any idea on how to resolve this ? HTML Code: In file included from FGMFoam.C:88:PVEqn.H: In function ‘int main(int, char)’:PVEqn.H:8:20: error: ‘class Foam::CombustionModel<Foam::psiReactionThermo>’ has no member named ‘SourcePV’ 8 \| reaction->SourcePV() //--- added by senbin, 2019-04-26 \| ^~~~~~~~ HTML Code: make: * [/usr/lib/openfoam/openfoam2406/wmake/rules/General/transform:38: Make/linux64GccDPInt32Opt/FGMFoamPost.o] Error 1 The source code of file is given below. https://drive.google.com/file/d/1ZL3...ew?usp=sharing

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
temp dependent density	Linas	Fluent UDF and Scheme Programming	2	February 13, 2025 19:22
Get the rate of reactions from the reactingFoam solver	Mirza8	OpenFOAM Post-Processing	10	June 24, 2024 03:37
Error: WorkBench Error: Could not handle event: SolutionStatusUpdate	Kieyo	Fluent Multiphase	0	November 10, 2022 00:58
Error SIGSEGV using VOF and UDF	JERC_UTFSM	Fluent UDF and Scheme Programming	14	November 8, 2021 00:17
Is it possible to accelerate the ode solver of reactingFoam?	pajofego	OpenFOAM Running, Solving & CFD	2	August 10, 2014 06:38

February 7, 2025, 05:02		#2
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	I suggest to use the version of OpenFoam for which the code was developed. See also post related to FGM on this forum.

February 7, 2025, 09:15		#3
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	I am not sure in which version it was written. Most files have following comment at top. /---------------------------------------------------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| \\ / A nd \| www.openfoam.com \\/ M anipulation \| ------------------------------------------------------------------------------- Copyright (C) 2011-2015 OpenFOAM Foundation ------------------------------------------------------------------------------

February 7, 2025, 11:52		#4
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	Not sure either. I believe something like OpenFoam-2.2 or 2.3. The alternative approach is to persue a new implementation. A function object to compute the mixture fraction is already available. Next step its to extend to. computation of the mixture fraction variance.

February 10, 2025, 06:43		#6
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	Not sure. Is the entry "cubeRootVol" by mistake repeated in some input file?

February 10, 2025, 07:34		#8
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	1/ Suggest to check the case files (system/ and constant/) using grep; 2/ Suggest to share entire log file (allowing to check how and where the error message appears);

February 10, 2025, 09:18		#9
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	log.zip Please find attached log files.

February 10, 2025, 09:27		#10
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	It is possible that errors are due to v2406 being used?

February 10, 2025, 09:36		#11
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	Not sure. I have tried compiling it in version v2206 as you recommended previously but I was getting same warnings. The solver is running fine but it gives these warnings which I am not sure why.

February 10, 2025, 09:45		#12
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 834 Blog Entries: 1 Rep Power: 19	I was recommending going back way earlier and compile against even older versions of openfoam. Apologies for not being sufficiently clear. What are the simulation results like?

February 10, 2025, 10:22		#13
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	Haven't had the chance to look at it properly but based on initial expectation, contour plots for temperature and OH radicals are quite close as given in this tutorial. https://www.tfd.chalmers.se/~hani/ku...rt_FGMFoam.pdf

February 12, 2025, 06:20		#14
killua1915 New Member Killua Zoldyck Join Date: Apr 2023 Posts: 15 Rep Power: 3	For anyone wondering how I got rid of these warnings, I renamed all the classes which were mentioned in the warning. Crude approach but it worked. I have tried copying library from src to my local and recompiling it with different name but not sure why it didn't work.