Null-collision meshless Monte-Carlo-Application to the validation of fast radiative transfer solvers embedded in combustion simulators

The Monte-Carlo method is often presented as a reference method for radiative transfer simulation when dealing with participating, inhomogeneous media. The reason is that numerical uncertainties are only of a statistical nature and are accurately evaluated by measuring the standard deviation of the Monte Carlo weight. But classical Monte-Carlo algorithms first sample optical thicknesses and then determine absorption or scattering locations by inverting the formal integral definition of optical thickness as an increasing function of path length. This function is only seldom analytically invertible and numerical inversion procedures are required. Most commonly, a volumic grid is introduced and optical properties within each cell are replaced by approximate homogeneous or linear fields. Simulation results are then sensitive to the grid and can no longer be considered as references. We propose a new algorithmic formulation based on the use of null-collisions that eliminate the need for numerical inversion: no volumic grid is required. Benchmark configurations are first considered in order to evaluate the effect of two free parameters: the amount of null-collisions, and the criterion used to decide at which stage a Russian Roulette is used to exit the path tracking process. Then the corresponding algorithm is implemented using a development environment allowing to deal with complex geometries (thanks to computer graphics techniques), leading to a Monte Carlo code that can be easily used for validation of fast radiative transfer solvers embedded in combustion simulators. ``Easily'' means here that the way the Monte Carlo algorithm deals with both the geometry and the temperature/pressure/concentration fields is independent of the choices made inside the combustion solver: there is no need for the design of a new path-tracking procedure adapted to each new CFD grid. The Monte Carlo simulator is ready for use as soon as combustion specialists provide a localization/interpolation tool defining what they consider as the continuous input fields best suiting their numerical assumptions. The radiation validation tool introduces no grid in itself.

[1]  Monte-Carlo algorithm for line-by-line calculations of thermal radiation in multiple scattering layered atmospheres , 2006 .

[2]  M. Pinar Mengüç,et al.  Thermal Radiation Heat Transfer , 2020 .

[3]  Jin Zhang.  Radiation Monte Carlo approcah dedicated to the coupling with LES reactive simulation. , 2011 .

[4]  M. Modest,et al.  An adaptive emission model for Monte Carlo simulations in highly inhomogeneous media represented by stochastic particle fields , 2007 .

[5]  Bénédicte Cuenot,et al.  Discrete Ordinates and Monte Carlo Methods for Radiative Transfer Simulation Applied to Computational Fluid Dynamics Combustion Modeling , 2009 .

[6]  J. Gautrais,et al.  Integral formulation of null-collision Monte Carlo algorithms , 2013 .

[7]  Frédéric André,et al.  Generalization of the k-moment method using the maximum entropy principle. Application to the NBKM and full spectrum SLMB gas radiation models , 2012 .

[8]  Bénédicte Cuenot,et al.  Analysis of Radiation Modeling for Turbulent Combustion: Development of a Methodology to Couple Turbulent Combustion and Radiative Heat Transfer in LES , 2011 .

[9]  Optimization of photon tracking in GATE , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[10]  Anouar Soufiani,et al.  High temperature gas radiative property parameters of statistical narrow-band model for H2O, CO2 and CO, and correlated-K model for H2O and CO2 , 1997 .

[11]  Greg Humphreys,et al.  Physically Based Rendering: From Theory to Implementation , 2004 .

[12]  J. Dufresne,et al.  A boundary-based net-exchange Monte Carlo method for absorbing and scattering thick media , 2005, 1902.09966.

[13]  D. Veynante,et al.  Monte Carlo method of radiative transfer applied to a turbulent flame modeling with LES , 2009 .

[14]  J. Dufresne,et al.  Net exchange parameterization of thermal infrared radiative transfer in Venus' atmosphere , 2009 .

[15]  Ömer L. Gülder,et al.  Application of the statistical narrow-band correlated-k method to non-grey gas radiation in CO , 2001 .

[16]  Michael F. Modest,et al.  Spectral Monte Carlo models for nongray radiation analyses in inhomogeneous participating media , 2007 .

[17]  Jorge Amaya,et al.  Analysis of the interaction between turbulent combustion and thermal radiation using unsteady coupled LES/DOM simulations , 2012 .

[18]  S. Tashkun,et al.  CDSD-4000: High-resolution, high-temperature carbon dioxide spectroscopic databank , 2011 .

[19]  Jonathan Tennyson,et al.  HITEMP, the high-temperature molecular spectroscopic database , 2010 .

[20]  Ephraim M Sparrow,et al.  Monte carlo radiation solutions—effect of energy partitioning and number of rays , 1973 .

[21]  Ö. Gülder,et al.  Application of the statistical narrow-band correlated-k method to low-resolution spectral intensity and radiative heat transfer calculations — effects of the quadrature scheme , 2000 .

[22]  Takayuki Tanaka,et al.  ARTS: Accelerated Ray-Tracing System , 1986, IEEE Computer Graphics and Applications.

[23]  Charles Meneveau,et al.  A priori testing of a similarity model for large eddysimulations of turbulent premixed combustion , 2002 .

[24]  Anouar Soufiani,et al.  Gas IR Radiative Properties: From Spectroscopic Data to Approximate Models , 1999 .

[25]  John R. Howell,et al.  Comparison of Monte Carlo Strategies for Radiative Transfer in Participating Media , 1998 .

[26]  A net-exchange Monte Carlo approach to radiation in optically thick systems , 2002 .

[27]  M. Modest.  Radiative heat transfer , 1993 .

[28]  Jean-François Cornet,et al.  Monte Carlo advances and concentrated solar applications , 2014 .

[29]  Bénédicte Cuenot,et al.  Comparison of three spatial differencing schemes in discrete ordinates method using three-dimensional unstructured meshes , 2005 .

[30]  Aldo Badano,et al.  Monte Carlo simulation of X-ray imaging using a graphics processing unit , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[31]  S. Blanco,et al.  The practice of recent radiative transfer Monte Carlo advances and its contribution to the field of microorganisms cultivation in photobioreactors , 2013 .

[32]  D. Veynante,et al.  Large eddy simulations of an acoustically excited turbulent premixed flame , 2000 .

[33]  Long-wave radiative analysis of cloudy scattering atmospheres using a net exchange formulation , 2004, 1902.08515.

[34]  D. Poitou.  Modélisation du rayonnement dans la simulation aux grandes échelles de la combustion turbulente , 2009 .