25–29 Mar 2019
Campus Puerta de Toledo of the Universidad Carlos III, Madrid, Spain
Europe/Amsterdam timezone

Radiation Modeling in the PIC-DSMC code PICLAS

26 Mar 2019, 10:00
30m
Campus Puerta de Toledo of the Universidad Carlos III, Madrid, Spain

Campus Puerta de Toledo of the Universidad Carlos III, Madrid, Spain

Puerta de Toledo Campus Ronda de Toledo, 1 28005 Madrid, Spain GPS coordinates: 40º24´30,24” N 3º42´39,59” O Metro: Puerta de Toledo Station (Line 5) Suburban train: Embajadores Station (Line C5) or Pirámides Station (Lines C1, C7 y C10)
Numerical Simulations Numerical Simulations

Speaker

Julian Beyer (University of Stuttgart)

Description

Abstract

Radiation processes are important for the description and understanding of plasma phenomena, since the radiative heat flux on the capsule can be substantial [1]. Experimental measurements are often complicated for several reasons. An alternative is the numerical simulation. Here, the common approach is the coupling of a CFD method with a radiation solver and a radiation transport solver. Especially the efficiency of these methods in dense flow regions, where radiation transfer becomes most important is an advantage. Nevertheless, the treatment of non-equilibrium effects in the flow field, which can strongly influence the radiative heat transfer, becomes problematic with CFD methods. To overcome this problem, a gas kinetic description of the plasma is necessary. This is often done by using the Direct Simulation Monte Carlo Method (DSMC) [2]. Previous approaches of coupling a particle code with a radiation solver [3, 4] revealed indeed promising results. Using DSMC, it is possible to produce detailed information about flow species, which can vary widely from the average. Detailed information about rotational, vibrational, and electronic excitation temperatures as well as the density of each species can be directly used as input variables for each cell of the radiation calculation, which would cause significant difficulties using CFD flow field data. A critical point of the radiation modeling is the solution of the Radiation Transport Equation (RTE). Here, different algorithms exist with different levels of accuracy and computational effort. In the Monte Carlo Method [5], the energy is divided into an integer number of particles. Their properties like wavelength, direction, and position are randomly assigned. The way of each beam through the cells is traced and the optical path is calculated. Subsequently, the RTE is solved in the direction of the beam. In this work, first results of solving the RTE for radiation data, which are calculated within the PIC-DSMC code PICLas, are presented.

Acknowledgement

The authors gratefully acknowledge funding provided by Airbus Defence and Space, by ArianeGroup and by the Deutsche Forschungsgemeinschaft (DFG) within the project “Partikelverfahren mit Strahlungslöser zur Simulation hochenthalper Nichtgleichgewichts-Plasmen” (project number 393159129).

References

  1. T. Ozawa, Improved Chemistry Models for DSMC Simulations of Ionized Rarefied Hypersonic Flows, Ph.D. Thesis, Department of Aerospace Engineering, Pennsylvania State University, University Park, PA, 2007.

  2. G. A. Bird, Molecular Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, New York 1994.

  3. I. D. Boyd, J. Zhong, D. A. Levin, and P. Jenniskens, Flow and radiation analysis for stardust entry at high altitude, in 46th AIAA Aerospace Sciences Meeting and Exhibit, page 1215, 2008.

  4. T. Ozawa, A. Wang, D. A. Levin, and M. Modest, Development of a Coupled DSMC-Particle Photon Monte Carlo Method for Simulating Atomic Radiation in Hypersonic Reentry Flows, in 40th Thermophysics Conference, page 3916, 2010.

  5. J. R. Howell, The Monte Carlo Method in Radiative Heat Transfer, ASME J. Heat Transfer, 120, pp. 547-560, 1998.

Primary author

Julian Beyer (University of Stuttgart)

Co-authors

Dr Marcel Pfeiffer (University of Stuttgart) Prof. Stefanos Fasoulas (University of Stuttgart)

Presentation materials