acuoci / laminarsmoke Goto Github PK
View Code? Open in Web Editor NEWCFD solver (based on OpenFOAM) for laminar reacting flow with detailed kinetic mechanisms based on OpenSMOKE++
CFD solver (based on OpenFOAM) for laminar reacting flow with detailed kinetic mechanisms based on OpenSMOKE++
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.
laminarSMOKE/solvers/laminarSMOKE/properties.H
Lines 232 to 235 in 65b29ff
My trial solutions)
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
<> 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?
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.
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!
Hello,
The links to the open source solvers (except Radau) no longer work.
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
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]
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
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
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?
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 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
The DRG must be adjusted because of the recent changes in the interfaces provided by OpenSMOKE++ library.
DVODE is not included in ODEPACK, and could not be applied, how to link it with opensmoke?
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.