CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

New solver: ddtFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree19Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 16, 2016, 04:06
Default
  #41
New Member
 
NGUYENVB
Join Date: Aug 2015
Posts: 11
Rep Power: 8
vanbonguyen is on a distinguished road
Quote:
Originally Posted by lx882211 View Post
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 likes this.
vanbonguyen is offline   Reply With Quote

Old   December 16, 2016, 04:23
Default
  #42
New Member
 
Mr.liu
Join Date: Sep 2012
Posts: 27
Rep Power: 11
lx882211 is on a distinguished road
Quote:
Originally Posted by vanbonguyen View Post
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?
lx882211 is offline   Reply With Quote

Old   January 6, 2017, 17:39
Default
  #43
New Member
 
Pankaj
Join Date: Nov 2010
Location: USA
Posts: 5
Rep Power: 13
pankajonline2 is on a distinguished road
Quote:
Originally Posted by dohnie View Post
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
pankajonline2 is offline   Reply With Quote

Old   April 7, 2017, 17:31
Default
  #44
New Member
 
Shuai Yuan
Join Date: Nov 2016
Posts: 29
Rep Power: 7
sy2516 is on a distinguished road
Did anyone successfully compile the DDTFoam Solver in OpenFOAM4.0?
sy2516 is offline   Reply With Quote

Old   April 7, 2017, 19:35
Default Wave-Transmissive Bc with ddtFoam
  #45
New Member
 
Pankaj
Join Date: Nov 2010
Location: USA
Posts: 5
Rep Power: 13
pankajonline2 is on a distinguished road
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.
pankajonline2 is offline   Reply With Quote

Old   April 10, 2017, 16:59
Default
  #46
New Member
 
Shuai Yuan
Join Date: Nov 2016
Posts: 29
Rep Power: 7
sy2516 is on a distinguished road
Quote:
Originally Posted by ecbmxer View Post
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.
sy2516 is offline   Reply With Quote

Old   October 8, 2017, 18:19
Default Update for current versions
  #47
New Member
 
James
Join Date: Oct 2017
Posts: 2
Rep Power: 0
jvkoch is on a distinguished road
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?
jvkoch is offline   Reply With Quote

Old   October 26, 2017, 15:44
Default Changing hydrogen fuel to methane in ddtFoam
  #48
New Member
 
Richard
Join Date: Oct 2017
Location: U.K.
Posts: 2
Rep Power: 0
Richardk is on a distinguished road
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.
Richardk is offline   Reply With Quote

Old   December 5, 2017, 08:23
Default problem while reading header for object cTable
  #49
New Member
 
Huy Quang Dong
Join Date: Nov 2017
Posts: 9
Rep Power: 6
cunconbkhp is on a distinguished road
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!
cunconbkhp is offline   Reply With Quote

Old   December 5, 2017, 14:41
Default Error lookuptIgn_06
  #50
New Member
 
vitor gabriel silva
Join Date: Nov 2015
Posts: 3
Rep Power: 8
vitorgabriel is on a distinguished road
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
Attached Files
File Type: zip methane.cti.zip (13.5 KB, 16 views)
File Type: zip lookupfiles.zip (5.0 KB, 12 views)
__________________
Vitor Gabriel Silva, Msc. Eng
Chemical and HSE Engineer

E-mail: vitorgabriel.eng@hotmail.com
vitorgabriel is offline   Reply With Quote

Old   November 12, 2018, 02:33
Default errors during installation of ddtfoam on the OpenFOAMŪ for Windows
  #51
New Member
 
Jian Gao
Join Date: Nov 2018
Posts: 1
Rep Power: 0
koken is on a distinguished road
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?

Last edited by koken; November 13, 2018 at 02:32.
koken is offline   Reply With Quote

Old   November 12, 2018, 16:56
Default
  #52
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,965
Blog Entries: 45
Rep Power: 126
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quote:
Originally Posted by koken View Post
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.
__________________
wyldckat is offline   Reply With Quote

Old   October 28, 2019, 00:36
Default Compilation Error
  #53
New Member
 
sanjeev adhikari
Join Date: May 2018
Posts: 13
Rep Power: 6
sanjeev_adhikari is on a distinguished road
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











Quote:
Originally Posted by dohnie View Post
Dear Foamers,
Quote:
Originally Posted by dohnie View Post
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/
sanjeev_adhikari is offline   Reply With Quote

Old   October 28, 2019, 04:06
Default
  #54
New Member
 
anonymous
Join Date: Nov 2010
Posts: 9
Rep Power: 13
forwardyang is on a distinguished road
Hi friends,

Does ddtFoam accept 3D unstructured mesh?
forwardyang is offline   Reply With Quote

Old   October 28, 2019, 22:24
Default
  #55
New Member
 
anonymous
Join Date: Nov 2010
Posts: 9
Rep Power: 13
forwardyang is on a distinguished road
according to the readme file, try of211
forwardyang is offline   Reply With Quote

Old   October 29, 2019, 03:18
Default combatibility
  #56
New Member
 
sanjeev adhikari
Join Date: May 2018
Posts: 13
Rep Power: 6
sanjeev_adhikari is on a distinguished road
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.
sanjeev_adhikari is offline   Reply With Quote

Old   December 23, 2020, 05:59
Default
  #57
New Member
 
Join Date: Oct 2020
Posts: 10
Rep Power: 3
li siye is on a distinguished road
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 is offline   Reply With Quote

Old   December 23, 2020, 06:01
Default
  #58
New Member
 
Join Date: Oct 2020
Posts: 10
Rep Power: 3
li siye is on a distinguished road
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 is offline   Reply With Quote

Old   November 20, 2021, 07:02
Default
  #59
New Member
 
Join Date: Nov 2021
Posts: 1
Rep Power: 0
Ruixuan is on a distinguished road
Quote:
Originally Posted by Alvin View Post
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.
Ruixuan is offline   Reply With Quote

Old   March 14, 2022, 08:20
Default
  #60
New Member
 
Umberto Somma
Join Date: Mar 2022
Posts: 1
Rep Power: 0
umbertosomma is on a distinguished road
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.
umbertosomma is offline   Reply With Quote

Reply

Tags
combustion, compressible, density based, detonation, transonic

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluent Adjoint Solver? ex10148 FLUENT 16 September 28, 2018 08:11
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
Divergence problem Smaras FLUENT 13 February 21, 2013 05:03
3d vof Smaras FLUENT 2 February 19, 2013 06:58
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 23:59.