Hi
I encounterd a problem when trying to make pimple branches with sundials package.
There is no problem with Simple packages.
Please help me with this problem.
Best regards

Dear developers and contributors,

I try to simulate a 1D premixed flame to get the laminar flame speed.
So I made 1D shape as below.

Boundary Conditions
Inlet: fixedValue for mass fractions
Outlet: InletOutlet for mass fractions
Sides: empty (to achieve 1D simulation)

Here is a problem. When I run this case, this error appears.

I checked ~laminarSMOKE/solvers/laminarSMOKE/properties.H, and it seems to be caused by some boundary faces that have zero mass fractions initially. “MoleFractions_From_MassFractions” does not work if all input mass fractions are zero.

My trial solutions)

1. Run reactingFoam by one timestep and restart with a new initial solution, so that zero mass fractions does not appear at boundary faces.
2. Or edit properties.H; (ex, If all mass fractions are zero, omit above lines.)

Two ways seem to work.
If there are better ways, please let me know.
It would be grateful if properties.H could work for this kind of situation.

Thank you

Best regards,
Daeyoung Jun

Hello Professor Antonio,
Sir, I am working on laminar methane flames. I am using ubuntu 14..04 and OpenFOAM 2.4. I had followed all the instructions given in Readme file of Laminar SMOKE. I am following mybashrc.minimalist. When I am compiling I am getting these errors as shown in the screenshot below. These errors are obtained during compilation of laminarBuoyantPimpleFoam. Is it a bug or error due to non-installation of libraries like(INTEL MKL).

Thanxs in Advance

Hi Alberto:

I am having difficulties compiling laminarSMOKE with OpenFOAM 4.1. The issue pertains to both laminarPimpleSMOKE and laminarSimpleSMOKE. I discuss the latter only.

Note that the rest of the solvers (i.e. laminarBuoyantPimpleSMOKE, laminarBuoyantSimpleSMOKE, and laminarSMOKEpostProcessor) compile correctly.

The compilation fails on solvers/laminarSimpleSMOKE as below.

It seems that laminarSimpleSMOKE.C requires pressureControl.H, which is available in OpenFoam-dev, but not in 4.1

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 -DOPENFOAM_VERSION=40 -w -DOPENSMOKE_USE_MKL=0 -DOPENSMOKE_USE_ISAT=0 -DDEVVERSION=0 -I../laminarSMOKE -I../laminarSMOKE/steady -I../../libs/radiationOpenSMOKE++/lnInclude -I../openSMOKEpp4laminarSMOKE/ -I/usr/lib64/include -I/usr/include/eigen3 -I/home/fbisetti/OpenFOAM/thirdparty-polimi/rapidxml-1.13 -I/include -I/include -I/ -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/finiteVolume/lnInclude -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/fvOptions/lnInclude -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/meshTools/lnInclude -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/sampling/lnInclude -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/thermophysicalModels/basic/lnInclude -IlnInclude -I. -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/OpenFOAM/lnInclude -I/apps/openfoam/4.1/gcc/system/openmpi-2.1.0/OpenFOAM-4.1/src/OSspecific/POSIX/lnInclude -fPIC -c laminarSimpleSMOKE.C -o Make/linux64GccDPInt32Opt/laminarSimpleSMOKE.o
laminarSimpleSMOKE.C:66:30: fatal error: pressureControl.H: No such file or directory
#include "pressureControl.H"
^
compilation terminated.
make: *** [Make/linux64GccDPInt32Opt/laminarSimpleSMOKE.o] Error 1

I am not sure how to fix this!

Thanks
Fabrizio

Dear Prof. acuoci,

I have been asking for Opensmoke++ libraries for testing current version of Laminarsmoke but it looks like you are modifying both at same time. Kindly upload the opensmoke++ libraries so that I can work on current version of Laminarsmoke.

I want to modify the current version to add conjugate heat transfer for solids.

Regards
Ratna Kishore

Hi,
I was able to get the code to compile successfully with the minimalist settings.
I'm trying to get the packaged tutorial case "MethaneCoflowFlame" to work for a CH4 global mechanism. At the first time step, the code exits with the following error:

Starting time loop

Courant Number mean: 6.41183184969e-05 max: 0.000150895373491
deltaT = 1.19999040008e-06
Time = 1.1999904e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 9.59486949027e-12, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.0201082144937, Final residual = 4.2849627155e-11, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.999999999844, Final residual = 3.69112647038e-10, No Iterations 1
 Umin: (-3.3731055518e-06 -6.95967895832e-06 0) Umax: (0.000159554715889 6.95967895832e-06 0.15)
DILUPBiCGStab:  Solving for CH4, Initial residual = 1, Final residual = 3.23802485131e-16, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 0.999932492575, Final residual = 1.33812383028e-16, No Iterations 1
DILUPBiCGStab:  Solving for CO2:  solution singularity
DILUPBiCGStab:  Solving for H2O:  solution singularity
Inert species is N2 with local index equal to 4
Transport equations of species solved in 0.004862 s 
DILUPBiCGStab:  Solving for T, Initial residual = 0.542046994922, Final residual = 2.76226339329e-10, No Iterations 1
 * T gas min/max (after transport)   = 300, 300
 * Solving homogeneous chemistry (OpenSMOKE++ solver, Direct integration)... 
terminate called after throwing an instance of 'rapidxml::parse_error'
  what():  unexpected end of data

Any idea about what this is due to? I am using the latest version of rapidxml and gcc on OpenFOAM-6.

Thanks!

The DRG must be adjusted because of the recent changes in the interfaces provided by OpenSMOKE++ library.

Hello Prof.
I am trying to cmpile the code with odepack and encounter with dscal_ redifinition.
The error is:
/home/ali/Desktop/odepack/libodepack.a(opkda2.o): In function dscal_': opkda2.f:(.text+0x86b): multiple definition of dscal_'
/opt/intel/mkl/lib/intel64/libmkl_intel_lp64.a(_dscal_lp64.o):_dscal.c:(.text+0x0): first defined here
collect2: error: ld returned 1 exit status
make: *** [/home/ali/OpenFOAM/ali-4.1/platforms/linux64GccDPInt32Opt/bin/laminarPimpleSMOKE] Error 1

What is the solution for this problem?
Best regards

The current implementation of ISAT does not work, because the ISATLib was recently updated by Mauro Bracconi.
We need to update the interfaces to correctly use the most recent version of ISATLib.

After the introduction of the customized radiation library, in case of wrong input data from the user, the solvers report the error message but an additional, annoying, error message about memory corruption.
This is not a big issue, because in any case the solver would abort, but it is better to remove it.

FOAM exiting

*** glibc detected *** laminarBuoyantPimpleSMOKE: double free or corruption (!prev): 0x0000000001bcde60 ***
======= Backtrace: =========
/lib64/libc.so.6[0x3b71276166]
/lib64/libc.so.6[0x3b71278c93]
/lib64/libc.so.6(__cxa_finalize+0x9d)[0x3b712361bd]
/software/OpenFOAM/OpenFOAM-2.4.x/platforms/linux64GccDPOpt/lib/libradiationModels.so(+0x4ca06)[0x7fac8b997a06]
======= Memory map: ========
00400000-00ec0000 r-xp 00000000 00:13 17218158037 /home/chimica2/cuoci/OpenFOAM/cuoci-2.4.x/platforms/linux64GccDPOpt/bin/laminarBuoyantPimpleSMOKE
010c0000-010e2000 rw-p 00ac0000 00:13 17218158037 /home/chimica2/cuoci/OpenFOAM/cuoci-2.4.x/platforms/linux64GccDPOpt/bin/laminarBuoyantPimpleSMOKE
010e2000-01126000 rw-p 00000000 00:00 0
01a21000-021b2000 rw-p 00000000 00:00 0 [heap]

DVODE is not included in ODEPACK, and could not be applied, how to link it with opensmoke?

<> Is it possible to simulate a case in a similar way to SandiaD (OpenFoam/combustion/reactingFoam), in which the flow develops at an initial time and then, the chemical part begins?

Hello,
The links to the open source solvers (except Radau) no longer work.

Except the already known declaration problems, there are also some problems which seems from interface.


I tried several version of sundials package, none of them worked.
Current set up for sundials-4.1.0.

Prof Cuoci,

There is a problem when I tried to compile "solver/laminarBuoyantPimpleSMOKE" according to the README.md in section 4, I got errors regarding undefined variables:

parcels
and
surfaceFilm

I thought maybe this is a bug/problem from the original git source repository, I commented the lines 75-76 regarding the two variable "parcels" and "surfaceFilm" in the file /home/pawans76/newver/laminarSMOKE/solvers/laminarBuoyantPimpleSMOKE/pEqn.22.H
[pawans76@mike5 laminarBuoyantPimpleSMOKE]$ pwd
/home/pawans76/newver/laminarSMOKE/solvers/laminarBuoyantPimpleSMOKE
[pawans76@mike5 laminarBuoyantPimpleSMOKE]$ diff pEqn.22.H pEqn.22.H.orig 75,78c75,77
< // parcels.Srho()
< // + surfaceFilm.Srho()
< // + fvOptions(psi, p_rgh, rho.name())
< fvOptions(psi, p_rgh, rho.name())
--- (in place of)

    parcels.Srho()
  + surfaceFilm.Srho()
  + fvOptions(psi, p_rgh, rho.name())

The above changes allow the solver laminarBuoyantPimpleSMOKE to compile.

Is this a bug?