2D granular flows with the μ(I)\mu(I) rheology and side walls friction: a well balanced multilayer discretization

We present here numerical modelling of granular flows with the $\mu(I)$ rheology in confined channels. The contribution is twofold: (i) a model to approximate the Navier-Stokes equations with the $\mu(I)$ rheology through an asymptotic analysis. Under the hypothesis of a one-dimensional flow, this model takes into account side walls friction; (ii) a multilayer discretization following Fern\'andez-Nieto et al. (J. Fluid Mech., vol. 798, 2016, pp. 643-681). In this new numerical scheme, we propose an appropriate treatment of the rheological terms through a hydrostatic reconstruction which allows this scheme to be well-balanced and therefore to deal with dry areas. Based on academic tests, we first evaluate the influence of the width of the channel on the normal profiles of the downslope velocity thanks to the multilayer approach that is intrinsically able to describe changes from Bagnold to S-shaped (and vice versa) velocity profiles. We also check the well balance property of the proposed numerical scheme. We show that approximating side walls friction using single-layer models may lead to strong errors. Secondly, we compare the numerical results with experimental data on granular collapses. We show that the proposed scheme allows us to qualitatively reproduce the deposit in the case of a rigid bed (i. e. dry area) and that the error made by replacing the dry area by a small layer of material may be large if this layer is not thin enough. The proposed model is also able to reproduce the time evolution of the free surface and of the flow/no-flow interface. In addition, it reproduces the effect of erosion for granular flows over initially static material lying on the bed. This is possible when using a variable friction coefficient $\mu(I)$ but not with a constant friction coefficient

Latitudinal distribution of the solar wind properties in the low- and high-pressure regimes: Wind observations

International audienceThe solar wind properties depend on ?, the heliomagnetic latitude with respect to the heliospheric current sheet (HCS), more than on the heliographic latitude. We analyse the wind properties observed by Wind at 1 AU during about 2.5 solar rotations in 1995, a period close to the last minimum of solar activity. To determine ?, we use a model of the HCS which we fit to the magnetic sector boundary crossings observed by Wind. We find that the solar wind properties mainly depend on the modulus |?|. But they also depend on a local parameter, the total pressure (magnetic pressure plus electron and proton thermal pressure). Furthermore, whatever the total pressure, we observe that the plasma properties also depend on the time: the latitudinal gradients of the wind speed and of the proton temperature are not the same before and after the closest HCS crossing. This is a consequence of the dynamical stream interactions. In the low pressure wind, at low |?|, we find a clear maximum of the density, a clear minimum of the wind speed and of the proton temperature, a weak minimum of the average magnetic field strength, a weak maximum of the average thermal pressure, and a weak maximum of the average ß factor. This overdense sheet is embedded in a density halo. The latitudinal thickness is about 5° for the overdense sheet, and 20° for the density halo. The HCS is thus wrapped in an overdense sheet surrounded by a halo, even in the non-compressed solar wind. In the high-pressure wind, the plasma properties are less well ordered as functions of the latitude than in the low-pressure wind; the minimum of the average speed is seen before the HCS crossing. The latitudinal thickness of the high-pressure region is about 20°. Our observations are qualitatively consistent with the numerical model of Pizzo for the deformation of the heliospheric current sheet and plasma sheet

Simulation and optimisation of terahertz emission from InGaAs and InP photoconductive switches

We simulate the terahertz emission from laterally-biased InGaAs and InP using a three-dimensional carrier dynamics model in order to optimise the semiconductor material. Incident pump-pulse parameters of current Ti:Sapphire and Er:fibre lasers are chosen, and the simulation models the semiconductor's bandstructure using parabolic Gamma, L and X valleys, and heavy holes. The emitted terahertz radiation is propagated within the semiconductor and into free space using a model based on the Drude-Lorentz dielectric function. As the InGaAs alloy approaches InAs an increase in the emitted power is observed, and this is attributed to a greater electron mobility. Additionally, low-temperature grown and ion-implanted InGaAs are modelled using a finite carrier trapping time. At sub-picosecond trapping times the terahertz bandwidth is found to increase significantly at the cost of a reduced emission power.Comment: 9 pages, 7 figure

Monolithic echo-less photoconductive switches for high-resolution terahertz time-domain spectroscopy

Interdigitated photoconductive (IPC) switches are convenient sources and detectors for terahertz (THz) time domain spectroscopy. However, reflection of the emitted or detected radiation within the device substrate can lead to echoes that inherently limits the spectroscopic resolution achievable. In this work, we design and realize low-temperature-grown-GaAs (LT-GaAs) IPC switches for THz pulse generation and detection that suppresses such unwanted echoes. This is realized through a monolithic geometry of an IPC switch with a metal plane buried at a subwavelength depth below the LT-GaAs surface. Using this device as a detector, and coupling it to an echo-less IPC source, enables echo-free THz-TDS and high-resolution spectroscopy, with a resolution limited only by the temporal length of the measurement governed by the mechanical delay line used