OpenFOAM Low Mach Number Formulation for Combustion
Hi Dear Friends,
I want to know if anybody tried to simulate combustion using low mach number model in OpenFOAM. I would be appreciate if anybody can help me Are there any solver for modeling reacting flow with Low mach number Formulation? 
fireFoam 1.6.0 is the one you wantLow ma number. Others are not as far as I know. They are incompressible or compressible, but obtaining the Low ma number solver only needs very small modification based on the compressible codes.

Thank you for response
I want to simulate combustion in tube (3D) and channel (2D). and I think reactingFOAM is a suitable solver but I dont know How to change this solver for Low mach number flow. Can you give me some help that I can develop this solver? thank you so much 
If your thermodynamic pressure can be considered constant throughout the simulation (i.e. not a piston engine), you'll find what you need to change here:
http://web.student.chalmers.se/group...SlidesOFW5.pdf This was for OF 1.5, but the change of reactingfoam2.1 is straight forward. The solver works good, and has generated a number of papers: Duwig, Christophe, Sébastien Ducruix, and Denis Veynante. "Studying the Stabilization Dynamics of Swirling Partially Premixed Flames by Proper Orthogonal Decomposition." Journal of Engineering for Gas Turbines and Power 134 (2012): 101501. Duwig, Christophe, et al. "Large Eddy Simulations of a piloted lean premix jet flame using finiterate chemistry." Combustion Theory and Modelling 15.4 (2011): 537568. K 
Hello.
You may see if it helps: http://code.google.com/p/lowmachspra...lowmachSolver/ I've adapted dieselFoam 1.7 for a spray simulation in low Mach number. "pd" is the hydrodynamic pressure and "p" is the thermodynamic one. I considered the case of "p" being constant only. 
Thank you for your response.

Boundary Condition Low mach number in combustion
Hi Friends
I am working on simulation combustion in channel. I change reactingfoam for solving low mach formulation. now I have some question about p. 1. I have a open system, and I set preff = 101325 (1 atm), Is it true? 2. for pd (pdynamics), I set Inlet boundary : zero gradient outlet boundary: constant 101235 but I think this is wrong because pd is calculated between 101325 , 101225 and according relation P=Pd+Preff so P is 2 atm Whereas in reference are said Pd<<preff which boundary condition must I use for pdynamic? Initial condition for pd is important? thanks a lot 
All times are GMT 4. The time now is 15:28. 