[vanbonguyen]

Originally Posted by lx882211

i am very interested in DDT, i want to simulated this amazing phenomenon, compared with a C3H8 fuel. However, i don't know how to generate the ignition delay and other information used during the simulation. Could you tell the method or some tips about that?

Thank you very much
Best Regards

The website (https://sourceforge.net/projects/ddtfoam/) that you can download both source code for ddtfoam and two files in matlab (
lookupY_05.m and


lookuptIgn_06.m)

the first matlab code is used to generate the species at different temperature and pressure for lookup table. The second matlab code is used to generate the ignition time. In order to do this, you need to install the STDtoolbox for matlab (http://shepherd.caltech.edu/EDL/publ...ml/SD_Toolbox/) where you can find both cantera and python.
After installation you just find your own chemical kinetic model and generate the file.cti. Put the file.cti in the data base and run the two matlab codes above to get the ignition time and species table.


Good luck,
NGUYEN VAN BO

[lx882211]

Originally Posted by vanbonguyen

The website (https://sourceforge.net/projects/ddtfoam/) that you can download both source code for ddtfoam and two files in matlab (
lookupY_05.m and


lookuptIgn_06.m)

the first matlab code is used to generate the species at different temperature and pressure for lookup table. The second matlab code is used to generate the ignition time. In order to do this, you need to install the STDtoolbox for matlab (http://shepherd.caltech.edu/EDL/publ...ml/SD_Toolbox/) where you can find both cantera and python.
After installation you just find your own chemical kinetic model and generate the file.cti. Put the file.cti in the data base and run the two matlab codes above to get the ignition time and species table.


Good luck,
NGUYEN VAN BO

Thank you very much for your reply, i will try it. Thank you.
If i want to change H2 to C3H8, where should i take care in these two matlab scripts?

[pankajonline2]

Originally Posted by dohnie

Dear Anh,
you should get the best result by using the waveTransmissive BC.
However, this will require some extra programming work. The standard wave transmissive BC in OpenFOAM is based on the assumption of a constant ratio of specific heats (gamma) and a flux field which is called phi in most OpenFOAM solvers. In ddtFoam, where the fluxes are computed using a Riemann solver, this flux corresponds to Riemann.rhoUFlux() divided by Riemann.rhoFlux().
I hope that helps!
Florian

Hi Florian,

Thanks for your ddt FOAM, I am using now.
I need to use wave-transmissive BC.

The openFOAM defines "phi= mass-flux=pho*velocity*Area" during calculation of advection velocity. Also in ddt Foam , the Riemen solver provides Riemann.rhoFlux()= pho*velocity*Area ; So, I set directly in BC dictionary phi=rhoFlux. I am bit confused regarding the above suggestion from you that "flux corresponds to Riemann.rhoUFlux() divided by Riemann.rhoFlux()."

Could you please, let me know whether "phi " will be = rhoFlux() only or Riemann.rhoUFlux() divided by Riemann.rhoFlux()." ??

Thanks
Pankaj

[sy2516]

Did anyone successfully compile the DDTFoam Solver in OpenFOAM4.0?

[pankajonline2]

Hi Yuan Sorry I am running ddtFoam on OF-2.1. R u working on ddtFOAM? Then I have a question for you. Did you use waveTransmissive B.C with ddtFoam? I am always getting different outlet pressure when I change the "lref" in waveTransmissive parameter . For coldFlow ddtFOAM +WaveTransmissive B.C give consistent outlet pressure.

But, with strong detonation it does give accurate outlet pressure using WaveTransmissive B.C.

Anyone has any comment.

[sy2516]

Originally Posted by ecbmxer

Gloq, that is exactly what I did. I installed OF 2.1.1 in parallel with 2.3 and created aliases to switch between the two. Surprisingly, it worked out perfectly and I can run all of the ddt routines. The parafoam integration isn't perfect between the two versions, but I just use regular paraview with 2.1.1

How to install OF 2.1.1 in parallel with OF4.0? Thank you.

[jvkoch]

I'm also interested in porting ddtFoam to current versions (4.x and 5.0) of OpenFOAM. From these comments, it appears that the solver only works with version 2.1.1, which is unfortunate. The compiling of the solvers and dependencies fail for current versions of OpenFOAM... that's expected. Any people interested in diving into sorting this out with me?

[Richardk]

Hi formers,

I tried to change the fuel from hydrogen to methane, so it requires to change the ignition tables with using (lookupY_05.m and lookuptIgn_06.m))as well as some extera changes in the code, such as:
Changing the UnburnedThermo.C and eUnburnedThermo.C, as well as laminar flame speed in the compute_sl.

I have done these changes but still I cannot run a test case and even in the first time step, it goes out of janaf limit!
And I used Gulder's equation for Methane's SL.

I will appreciate if anyone can help me with these issue.


Thanks.

[cunconbkhp]

Hello ddtFoamers!

I have installed the ddtFoam solver with OpenFoam 2.1.1. I'm on ubuntu 16.

When i run the ddtFoam tutorial i have this error :

problem while reading header for object cTable

Does anyone knows the solution for this issue?

Thanks in advance!

[vitorgabriel]

Hello,

Thanks for ddtfoam, it's very good solver. But I'm having a problem to generate the igntion time:

Error using ThermoPhase/setTemperature (line 5)
the temperature must be positive
Error in Solution/set (line 152)
setTemperature(a,tval);
Error in uvsys (line 30)
set(gas, 'T', y(1), 'Rho', density(gas), 'Y', y(2:end));
Error in ode15s (line 579)
rhs = hinvGak*feval(odeFcn,tnew,ynew,odeArgs{:}) - Mtnew*(psi+difkp1);
Error in explosion (line 42)
out = ode15s(@uvsys,tel,y0,options,gas,mw);
Error in lookuptIgn_06CH4 (line 64)
[out] = explosion(gas0,fig_num); % constant volume explosion

I installed cantera and SDToolbox, everything is working. I generated the ddtfoamcTable_fp for methane base in cti for methane.

But the ignition time I trying all the ways and I dont solve.

Any help is welcome.

Thanks

[koken]

To compile, run the Allwmake script in the main directory. But，there are many errors, like this:

fatal error: primitivePatch.H: No such file or directory
compilation terminated.

I do not why，who can help me?

[wyldckat]

Originally Posted by koken

I do not why？who can help me?

Quick answer: Anyone can potentially help you, as long as you give enough information for people to help you
See for example the instructions given on this page: How to give enough info to get help


More specifically, in your post/message is missing the critical information of which OpenFOAM version/variant/fork you are using.
Knowing which Operating System in which you are using OpenFOAM can also help, along with knowing which installation instructions you've followed.

[sanjeev_adhikari]

I was trying to compile the ddtfoam in Openfoam V7,6, and 4 but there was an compilation error? Which version of Openfoam is compatible for ddtfoam?


+ wmake libso basic
./src/thermophysicalModels/Allwmake: 6: ./src/thermophysicalModels/Allwmake: wmake: not found
+ wmake libso reactionThermo
./src/thermophysicalModels/Allwmake: 7: ./src/thermophysicalModels/Allwmake: wmake: not found
+ wmake libso chemistryModel
./src/thermophysicalModels/Allwmake: 8: ./src/thermophysicalModels/Allwmake: wmake: not found
./Allwmake: 5: ./Allwmake: wmake: not found
./Allwmake: 6: ./Allwmake: wmake: not found
./Allwmake: 7: ./Allwmake: wmake: not found
./Allwmake: 8: ./Allwmake: wmake: not found

Originally Posted by dohnie

Dear Foamers,

Originally Posted by dohnie

I developed an OpenFOAM solver for the simulation of the deflagration-to-detonation transition in gases.
Some features:

• accurate shock capturing with a density based solver (largely based on Oliver Borm's development densityBasedTurbo)
• deflagration modelling via a reaction progress variable (based on XiFoam)
• auto-ignition sub-grid model separating each computational cell into a "shocked" and an "unshocked" part
• auto-ignition from tabulated ignition delay times, gained from a detailed chemical mechanism
• low computational cost due to applicability to coarse grids and usage of pre-tabulated chemistry

Even if you're not interested in detonations/combustion, it might be a nice alternative to existing OpenFOAM solvers when it comes to shock resolution in transonic flows.

The code incl. a short documentation and tutorials can be downloaded from
http://sourceforge.net/projects/ddtfoam/

[forwardyang]

Hi friends,

Does ddtFoam accept 3D unstructured mesh?

[forwardyang]

according to the readme file, try of211

[sanjeev_adhikari]

Does anybody tell me which version of openfoam is ddtfoam compatible? As i tired compiling in v4,6,7. However, had save error during compiling.

[li siye]

Dear Florian, it's a amazing work. can i ask how many number of grids for the case? and can the ddtFoam be used with dynamicRefineFvMesh? with Adaptive Refinement mesh can get lower cost for computation.

[li siye]

Dear Florian, it's a amazing work. can i ask how many number of grids for the case? and can the ddtFoam be used with dynamicRefineFvMesh? with Adaptive Refinement mesh can get lower cost for computation

[Ruixuan]

Originally Posted by Alvin

Hello, I installed ddtFoam onto OpenFoam 2.1.1 (The recommended version) in a SUSE Linux environment and trying to run the tutorial cases, but segmentation faults shows up and it stopped computing. (Attached picture. This is an error for the shocktube tutorial.)
According to the previous posts in this thread, it looks like everyone is using without any issues. Is there anyone had the same issues as mine? If you had and solved it, could you share how to fix it? Any information you can provide me would be greatly appreciated.

Hi Alvin, may I ask if you have solved the segmentation error and how? I am facing a similar problem.

[umbertosomma]

Hello! I am trying to run the tutorial simulations and I have a problem. I am following the short documentation "Go to the directory tutorials/pddtFoam Tutorial, run the Setup script and examine the initialconditions. Start running the case by invoking pddtFoam ". Basically when I try to run the case by invoking pddtFoam I get like '' command not found ''. I am a beginner of OpenFoam and CFD . Maybe you can give me a tip.
Thanks a lot.