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 want--Low 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:
This was for OF 1.5, but the change of reactingfoam-2.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 finite-rate chemistry." Combustion Theory and Modelling 15.4 (2011): 537-568.
You may see if it helps:
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
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 05:31.|