### Abstract

The equations of the Shell model are reexamined with a view to their more effective implementation into a computational fluid dynamics code. The simplification of the solution procedure without compromising accuracy is achieved by replacing time as an independent variable with the fuel depletion, which is the difference between the initial fuel concentration and the current one. All the other variables used in this model, including temperature, concentration of oxygen, radicals, intermediate and branching agents are expressed as functions of fuel depletion. Equations for the temperature and concentration of the intermediate agent are of the first order and allow analytical solutions. The concentrations of oxygen and fuel are related via an algebraic equation which is solved in a straightforward way. In this case the numerical solution of five coupled first-order ordinary differential equations is reduced to the solution of only two coupled first- order differential equations for the concentration of radicals and branching agent. It is possible to rearrange these equations even further so that the equation for the concentration of the radicals is uncoupled from the equation for the branching agent. In this case the equation for the concentration of radicals becomes a second-order ordinary differential equation. This equation is solved analytically in two limiting cases and numerically in the general case. The solution of the first-order ordinary differential equation for the concentration of the branching agent and the solution of the first-order differential equation for time are presented in the form of integrals containing the concentration of the radicals obtained earlier. This approach allows the central processing unit (CPU) time to be more than halved and makes the calculation of the autoignition process using the Shell model considerably more effective.

Original language | English |
---|---|

Pages (from-to) | 529-540 |

Number of pages | 12 |

Journal | Combustion and Flame |

Volume | 117 |

Issue number | 3 |

DOIs | |

Publication status | Published - May 1999 |

Externally published | Yes |

### Fingerprint

### ASJC Scopus subject areas

- Energy Engineering and Power Technology
- Fuel Technology
- Mechanical Engineering

### Cite this

*Combustion and Flame*,

*117*(3), 529-540. https://doi.org/10.1016/S0010-2180(98)00072-8

**The Shell autoignition model : A new mathematical formulation.** / Sazhin, S. S.; Sazhina, E. M.; Heikal, M. R.; Marooney, C.; Mikhalovsky, S. V.

Research output: Contribution to journal › Article

*Combustion and Flame*, vol. 117, no. 3, pp. 529-540. https://doi.org/10.1016/S0010-2180(98)00072-8

}

TY - JOUR

T1 - The Shell autoignition model

T2 - A new mathematical formulation

AU - Sazhin, S. S.

AU - Sazhina, E. M.

AU - Heikal, M. R.

AU - Marooney, C.

AU - Mikhalovsky, S. V.

PY - 1999/5

Y1 - 1999/5

N2 - The equations of the Shell model are reexamined with a view to their more effective implementation into a computational fluid dynamics code. The simplification of the solution procedure without compromising accuracy is achieved by replacing time as an independent variable with the fuel depletion, which is the difference between the initial fuel concentration and the current one. All the other variables used in this model, including temperature, concentration of oxygen, radicals, intermediate and branching agents are expressed as functions of fuel depletion. Equations for the temperature and concentration of the intermediate agent are of the first order and allow analytical solutions. The concentrations of oxygen and fuel are related via an algebraic equation which is solved in a straightforward way. In this case the numerical solution of five coupled first-order ordinary differential equations is reduced to the solution of only two coupled first- order differential equations for the concentration of radicals and branching agent. It is possible to rearrange these equations even further so that the equation for the concentration of the radicals is uncoupled from the equation for the branching agent. In this case the equation for the concentration of radicals becomes a second-order ordinary differential equation. This equation is solved analytically in two limiting cases and numerically in the general case. The solution of the first-order ordinary differential equation for the concentration of the branching agent and the solution of the first-order differential equation for time are presented in the form of integrals containing the concentration of the radicals obtained earlier. This approach allows the central processing unit (CPU) time to be more than halved and makes the calculation of the autoignition process using the Shell model considerably more effective.

AB - The equations of the Shell model are reexamined with a view to their more effective implementation into a computational fluid dynamics code. The simplification of the solution procedure without compromising accuracy is achieved by replacing time as an independent variable with the fuel depletion, which is the difference between the initial fuel concentration and the current one. All the other variables used in this model, including temperature, concentration of oxygen, radicals, intermediate and branching agents are expressed as functions of fuel depletion. Equations for the temperature and concentration of the intermediate agent are of the first order and allow analytical solutions. The concentrations of oxygen and fuel are related via an algebraic equation which is solved in a straightforward way. In this case the numerical solution of five coupled first-order ordinary differential equations is reduced to the solution of only two coupled first- order differential equations for the concentration of radicals and branching agent. It is possible to rearrange these equations even further so that the equation for the concentration of the radicals is uncoupled from the equation for the branching agent. In this case the equation for the concentration of radicals becomes a second-order ordinary differential equation. This equation is solved analytically in two limiting cases and numerically in the general case. The solution of the first-order ordinary differential equation for the concentration of the branching agent and the solution of the first-order differential equation for time are presented in the form of integrals containing the concentration of the radicals obtained earlier. This approach allows the central processing unit (CPU) time to be more than halved and makes the calculation of the autoignition process using the Shell model considerably more effective.

UR - http://www.scopus.com/inward/record.url?scp=0033136125&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0033136125&partnerID=8YFLogxK

U2 - 10.1016/S0010-2180(98)00072-8

DO - 10.1016/S0010-2180(98)00072-8

M3 - Article

VL - 117

SP - 529

EP - 540

JO - Combustion and Flame

JF - Combustion and Flame

SN - 0010-2180

IS - 3

ER -