The non-equilibrium Marshak wave problem in non-homogeneous media

Nitay Derei,Shmuel Balberg,Shay I. Heizler,Elad Steinberg,Ryan G. McClarren,Menahem Krief

2024-11-22

Abstract:We derive a family of similarity solutions to the nonlinear non-equilibrium Marshak wave problem for an inhomogeneous planar medium which is coupled to a time dependent radiation driving source. We employ the non-equilibrium gray diffusion approximation in the supersonic regime. The solutions constitute a generalization of the non-equilibrium nonlinear solutions that were developed recently for homogeneous media. Self-similar solutions are constructed for a power law time dependent surface temperature, a spatial power law density profile and a material model with power law temperature and density dependent opacities and specific energy density. The extension of the problem to non-homogeneous media enables the existence of similarity solutions for a general power law specific material energy. It is shown that the solutions exist for specific values of the temporal temperature drive and spatial density exponents, which depend on the material exponents. We also illustrate how the similarity solutions take various qualitatively different forms which are analyzed with respect to various parameters. Based on the solutions, we define a set of non-trivial benchmarks for supersonic non-equilibrium radiative heat transfer. The similarity solutions are compared to gray diffusion simulations as well as to detailed implicit Monte-Carlo and discrete-ordinate transport simulations in the optically-thick regime, showing a great agreement, which highlights the benefit of these solutions as a code verification test problem.

High Energy Astrophysical Phenomena,Mathematical Physics,Computational Physics,Plasma Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the propagation of non - equilibrium Marshak waves in non - uniform media. Specifically, the paper derives a series of self - similar solutions by introducing a more general setting, that is, in a planar medium, the medium has a spatial power - law density distribution in the form of \(\rho(x)=\rho_0x^{-\omega}\), the surface temperature drive varies with time in the form of a power - law \(T_s(t) = T_0t^\tau\), and the material model follows a power - law relationship between temperature and density, including total opacity \(k_t(T,\rho)=\frac{1}{G}T^{-\alpha}\rho^{1 + \lambda}\), absorption opacity \(k_a(T,\rho)=\frac{1}{G'}T^{-\alpha'}\rho^{1+\lambda'}\), and material energy density \(u(T,\rho)=FT^\beta\rho^{1-\mu}\).

These solutions allow self - similarity when \(\beta\neq4\) and \(\omega\neq0\), thus directly generalizing the solutions in previous works that assume \(\beta = 4\) and \(\omega = 0\). The paper shows that this new family of nonlinear self - similar solutions exists only when the time index \(\tau\) and the spatial index \(\omega\) satisfy specific values related to the material model indices \(\alpha,\lambda,\alpha',\lambda',\beta,\mu\), and analyzes in detail the properties and behaviors of these solutions under different parameters. Finally, the paper uses these generalized solutions to define a set of non - trivial benchmark tests for supersonic non - equilibrium radiative heat conduction, and compares these solutions with detailed numerical stochastic and deterministic radiation transfer simulations as well as gray - diffusion simulations.