In this paper, we discuss dissipation process of the binary mixture gas in the thermally relativistic flow \textcolor{red}{by focusing on the characteristics of the diffusion flux}. As an analytical object, we consider the relativistic rarefied-shock layer problem around the triangle prism. Numerical results of the diffusion flux are compared with the Navier-Stokes-Fourier (NSF) order approximation of the diffusion flux, which is calculated using the diffusion and thermal-diffusion coefficients by Kox \textit{et al}. [Physica A, 84, 1, pp.165-174 (1976)]. In the case of the uniform flow with the small Lorentz contraction, the diffusion flux, which is obtained by calculating the relativistic Boltzmann equation, is roughly approximated by the NSF order approximation inside the shock wave, whereas the diffusion flux in the vicinity of the wall is markedly different from the NSF order approximation. The magnitude of the diffusion flux, which is obtained by calculating the relativistic Boltzmann equation, is similar to that of the NSF order approximation inside the shock wave, unlike the pressure deviator, dynamic pressure and heat flux, even when the Lorentz contraction in the uniform flow becomes large, because the diffusion flux does not depend on the generic Knudsen number from its definition in Eckart's frame. \textcolor{red}{Finally, the author concludes that the accurate diffusion flux must be calculated from the particle four flow, which is formulated using the four velocity distinguished by each species of particles.

We discuss the thermal conduction and convection of thermally relativistic fluids between two parallel walls under the gravitational force, both theoretically and numerically. In the theoretical discussion, we assume that the Lorentz contraction is ignored and spacetime is flat. For understanding of the thermal conduction and convection of thermally relativistic fluids between two parallel walls under the gravitational force, we solve the relativistic Boltzmann equation using the direct simulation Monte Carlo method. Numerical results indicate that strongly nonequilibrium states are formed in vicinities of two walls, which do not allow us to discuss the transition of the thermal conduction to the thermal convection of thermally relativistic fluids under the gravitational force in the framework of the relativistic Navier-Stokes-Fourier equation, when the flow-field is under the transition regime between the rarefied and continuum regimes, whereas such strongly nonequilibrium states are not formed in vicinities of two walls under the nonrelativistic limit.

The previously developed Flux-Vector-Splitting (FVS) method was formulated here for the first-order version of Z4 formalism. Then, the characteristics of this method formulated for Z4 formalism were studied by numerically analyzing the evolution of two types of black holes (free and stuffed). Finally, these numerical results from the FVS method were compared with those from the Local Lax Friedrichs and Modified Local Lax Friedrichs methods to reveal the dependency of numerical solution by Z4 formalism on the choice of numerical scheme.

The relativistic shock layer problem was numerically analyzed by using two relativistic Boltzmann-kinetic equations. One is Marle model, and the other is Anderson-Witting model. As with Marle model, the temperature of the gain term was determined from its relation with the dynamic pressure in the framework of 14-moments theory. From numerical results of the relativistic shock layer problem, behaviors of projected moments in the nonequilibrium region were clarified. Profiles of the heat flux given by Marle model and Anderson-Witting model were quite adverse to the profile of the heat flux approximated by Navier-Stokes-Fourier law. On the other hand, profiles of the heat flux given by Marle model and Anderson-Witting model were similar to the profile approximated by Navier-Stokes-Fourier law. Additionally we discuss the differences between Anderson-Witting model and Marle model by focusing on the fact that the relaxational rate of the distribution function depends on both flow velocity and molecular velocity for Anderson-Witting model, while it depends only on the molecular velocity for Marle model.