Back to Computational Fluid Dynamics  |  Previous paper  |  Next paper

Predicting the Retention Time of Nuclear Reaction Products in the PSI Recoil Chamber Using COMSOL Multiphysics


R. Dressler1 R. Eichler1
1Paul Scherrer Institute, Villigen, Switzerland


CAD drawing of the PSI recoil chamber and mesh of the COMSOL Multiphysics model.

Introduction: The chemical properties of the heaviest elements (atomic number Z > 103) depend on the influence of the high nuclear charge to their electronic structure. Enhanced chemical stability of copernicium (Cn, Z = 112) and flerovium (Fl, Z = 114) was predicted already 40 year ago by Pitzer [1]. The challenge of chemical investigations of these elements is the tiny production rates of few atoms per week using 48Ca induced nuclear fusion reactions obfuscated with short half-lifes in the order of seconds and below (see e.g. [2]). The formed evaporation residues (EVRs) recoil out of the target and are thermalized in a recoil chamber (RC). Their retention time in the RC is comparable to the half-life of the produced isotopes. Recently performed experiments with Cn and Fl [3,4] showed that gas phase chemistry is suitable to meet all these requirements. In Figure 1 our experimental setup is schematically depicted. Detailed knowledge of the flow pattern and transport regime in the RC is essential to maximize the transport efficiency and the number of atoms available for chemical investigations. In contrast to the capillary transport, the dwell time prediction of reaction products out of the RC needs more efforts. Computational Methods: The steady state flow pattern of the carrier gas was calculated in a full 3D simulation using the CFD Module of COMSOL Multiphysics. The used geometry of the gas volume is based on a CAD model of the PSI-RC (Figure 2). Particular attention was paid to keep the original structure up to small details like imprints of skews or quartz columns. The final model consists of 504,298 mesh elements, and 466,377 degrees of freedom. In a following step the time dependence of the transport of a diluted species out of the RC was calculated, where the initial distribution of the EVRs was determined using the stopping range program SRIM-2011 [5]. Results: Due to the Z-limitations of the SRIM code uranium was used as target material and radon as EVR taken into account the known deviations from experimental results [6]. The obtained ion distribution was fitted using a Gaussian distribution in perpendicular direction to the recoiling EVRs and a Rayleigh type distribution along the direction of the recoiling EVRs, respectively.
The flow calculations were performed with a constant total carrier gas flow of 1.5 l/min, where the ratio of the flows entering inflow tube 1 (in1) and inflow tube 2 (in2) was varied. Figure 3 shows as an example of the isotope distributions after 75 ms if 20% and 80% flow passing inlet in1, respectively. Conclusions: These results were used to estimate the yield for interesting EVRs, in particular for the isotopes 287Fl (t1/2 ~ 0.5 s), 288Fl (t1/2 ~ 0.8 s), and 289Fl (t1/2 ~ 2.6 s). The results are depicted in figure 4. A flat yield maximum is predicted if 60% - 90% of the gas flow entering the inlet in1.

Share it on Social Media: