Investigating in the laboratory the process of matter accretion onto forming stars through scaled experiments is important in order to better understand star and planetary systems formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A first step has been made in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528]. in allowing such investigations. It revealed the existence of a two components stream: a hot shell surrounding a cooler inner stream. The shell was formed by matter laterally ejected upon impact and refocused by the local magnetic field. That laboratory investigation was limited tonormal incidence impacts. However, in young stellar objects, complex structure of magnetic fields causes variability of the incidentangles of the accretion columns. This led us to undertake an investigation, using laboratory plasmas, of the consequence of having a slanted accretion impacting a young star. Here we use high power laser interactions and strong magnetic field generation in the laboratory, complemented by numerical simulations, to study the asymmetry induced upon accretion structures when columns of matter impact the surface of young stars with an oblique angle. Compared to the scenario where matter accretes normal to the star surface, we observe strongly asymmetric plasma structure, strong lateral ejecta of matter, poor confinement of the accreted material and reduced heating compared to the normal incidence case. Thus, slanted accretion is a configuration that seems to be capable of inducing perturbations of the chromosphere and hence possibly influence the level of activity of the corona.

Models and observations indicate that the impact of matter accreting onto the surface of young stars produces regions at the base of accretion columns, in which optically thin and thick plasma components coexist. Thus an accurate description of these impacts requires to account for the effects of absorption and emission of radiation. We study the effects of radiation emerging from shock-heated plasma in impact regions on the structure of the pre-shock downfalling material. We investigate if a significant absorption of radiation occurs and if it leads to a pre-shock heating of the accreting gas. We developed a radiation hydrodynamics model describing an accretion column impacting onto the surface of a Classical T Tauri Star. The model takes into account the stellar gravity, the thermal conduction, and the effects of both radiative losses and absorption of radiation by matter in the non local thermodynamic equilibrium regime. After the impact, a hot slab of post-shock plasma develops at the base of the accretion column. Part of radiation emerging from the slab is absorbed by the pre-shock accreting material. As a result, the pre-shock accretion column gradually heats up to temperatures of $10^5$ K, forming a radiative precursor of the shock. The precursor has a thermal structure with the hottest part at $T \approx 10^5$ K, with size comparable to that of the hot slab, above the post-shock region. At larger distances the temperature gradually decreases to $T \approx 10^4$ K.

Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those among them that use either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to the local thermodynamic equilibrium (LTE). This assumption may be not valid in some astrophysical cases. We present an upgraded version of the LTE radiation-hydrodynamics module implemented in the PLUTO code, originally developed by Kolb et al. (2013), which we have extended to handle non-LTE regimes. Starting from the general frequency-integrated comoving-frame equations of radiation hydrodynamics (RHD), we have justified all the assumptions made to obtain the non-LTE equations actually implemented in the module, under the FLD approximation. An operator-split method is employed, with two substeps: the hydrodynamic part is solved with an explicit method by the solvers already available in PLUTO, the non-LTE radiation diffusion and energy exchange part is solved with an implicit method. The module is implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: the radiative power loss, the Planck and Rosseland mean opacities. Our implementation has been validated through different tests, in particular radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that non-LTE approach is of paramount importance to properly model accretion shock structures. Our radiation FLD module represents a step toward the general non-LTE RHD modeling. The module is available, under request, for the community.