We report on two- and three-dimensional numerical simulations of Rayleigh–Taylor instabilities in immiscible fluids. A diffuse-interface model that combines the Cahn–Hilliard equation, governing the evolution of the volume fraction of one fluid, and the Navier–Stokes equations, governing the bulk velocity and pressure, is used. The study is limited to low Atwood numbers owing to the use of the Boussinesq approximation. The code is based on a pseudo-spectral method. A linear analysis is first performed in a two-dimensional case of Rayleigh–Taylor instability to confirm that the model very well captures this phenomenon in the case of inviscid or viscid fluids. One key aspect of this work is that the influence of the thermodynamic parameters related to the Cahn–Hilliard equation (interface thickness and mobility) is quantitively studied. Three-dimensional results of Rayleigh–Taylor instabilities in viscous fluids are then presented to show the possibilities of this modeling. We observe the effect of the viscosity and the wavelength of an initial single-mode perturbation on the mass transport during the nonlinear regime.

The future PLINIUS-2 platform of CEA Cadarache will be dedicated to the study of corium interactions in severe nuclear accidents, and will host innovative large-scale experiments. The Nuclear Measurement Laboratory of CEA Cadarache is in charge of real-time high-energy X-ray imaging set-ups, for the study of the corium-water and corium-sodium interaction, and of the corium stratification process. Imaging such large and high-density objects requires a 15 MeV linear electron accelerator coupled to a tungsten target creating a high-energy Bremsstrahlung X-ray flux, with corresponding dose rate about 100 Gy/min at 1 m. The signal is detected by phosphor screens coupled to high-framerate scientific CMOS cameras. The imaging set-up is established using an experimentally-validated home-made simulation software (MODHERATO). The code computes quantitative radiographic signals from the description of the source, object geometry and composition, detector, and geometrical configuration (magnification factor, etc.). It accounts for several noise sources (photonic and electronic noises, swank and readout noise), and for image blur due to the source spot-size and to the detector unsharpness. In a view to PLINIUS-2, the simulation has been improved to account for the scattered flux, which is expected to be significant. The paper presents the scattered flux calculation using the MCNP transport code, and its integration into the MODHERATO simulation. Then the validation of the improved simulation is presented, through confrontation to real measurement images taken on a small-scale equivalent set-up on the PLINIUS platform. Excellent agreement is achieved. This improved simulation is therefore being used to design the PLINIUS-2 imaging set-ups (source, detectors, cameras, etc.). Simulation will be presented for different configurations, as well as preliminary tests of alternative detectors considered for PLINIUS-2.