Search

Efficient numerical models are derived for problems of natural convection and material solidification in a horizontal differentially heated slender cavity. These 2D1/2 models are obtained by averaging the equations of momentum, heat, and mass conservation along the transverse direction assuming both a constant temperature and a well defined velocity profile in this direction. Based on our former works, the transverse velocity profile is assumed to be either a Poiseuille profile ( 2 D p 1 / 2 model), or Hartmann-type profiles featuring two boundary layers on the sides of a uniform bulk ( 2 D H 1 / 2 model). For this 2 D H 1 / 2 model, however, a parameter δ (giving the boundary layer thickness) has to be adjusted: optimal values have been found in a large range of the control parameters and expressed as a reliable fitted function of G r . The ability of the model to reproduce 3D results in a 2D framework is investigated in a large range of the control parameters (Prandtl number P r and Grashof number G r ); the validity domain of the model in this parameter space is also clarified and rigorously defined. A good precision is obtained for natural convection problems (intensity of the flow, temperature field) as well as for solid-liquid phase change problems (shape, position, and evolution of the front). A comparison with unpublished experimental data of solidification of pure tin is also conducted. The boundary conditions for the simulation are first defined after a post-treatment of the time-dependent experimental data in order for them to be representative of the experimental process despite a significant and time dependent thermal resistance between the walls of the crucible and the liquid. A very good agreement is observed between the 2 D H 1 / 2 model and the experimental measurements for this pure tin solidification experiment in the AFRODITE setup.

International audience; The present paper focuses on solute segregation occurring in directional solidification processes with sharp solid/liquid interface, like silicon crystal growth. A major difficulty for the simulation of such processes is their inherently multi-scale nature: the impurity segregation problem is controlled at the solute boundary layer scale (micrometers) while the thermal problem is ruled at the crucible scale (meters). The thickness of the solute boundary layer is controlled by the convection regime and requires a specific refinement of the mesh of numerical models. In order to improve numerical simulations, wall functions describing solute boundary layers for convecto-diffusive regimes are derived from a scaling analysis. The aim of these wall functions is to obtain segregation profiles from purely thermo-hydrodynamic simulations , which do not require solute boundary layer refinement at the solid/liquid interface. Regarding industrial applications, various stirring techniques can be used to enhance segregation, leading to fully turbulent flows in the melt. In this context, the scaling analysis is further improved by taking into account the turbulent solute transport. The solute boundary layers predicted by the analytical model are compared to those obtained by transient segregation simulations in a canonical 2D lid driven cavity configuration for validation purposes. Convective regimes ranging from laminar to fully turbulent are considered. Growth rate and molecular diffusivity influences are also investigated. Then, a procedure to predict concentration fields in the solid phase from a hydrodynamic simulation of the solidification process is proposed. This procedure is based on the analytical wall functions and on solute mass conservation. It only uses wall shear-stress profiles at the solidification front as input data. The 2D analytical concentration fields are directly compared to the results of the complete simulation of segregation in the lid driven cavity configuration. Finally, an additional output from the analytical model is also presented. We put in light the correlation between different species convecto-diffusive behaviour; we use it to propose an estimation method for the segregation parameters of various chemical species knowing segregation parameters of one specific species.

The main objective of this work is to study experimentally the primary instability of non-Newtonian film flows down an inclined plane. We focus on low-concentration shear-thinning aqueous solutions obeying the Carreau law. The experimental study essentially consists of measuring wavelengths in marginal conditions, which yields the primary stability threshold for a given slope. The experimental results for neutral curves presented in the $(Re,f_{c})$ and $(Re,k)$ planes (where $f_{c}$ is the driving frequency, $k$ is the wavenumber and $Re$ is the Reynolds number) are in good agreement with the numerical results obtained by a resolution of the generalized Orr–Sommerfeld equation. The long-wave asymptotic extension of our results is consistent with former theoretical predictions of the critical Reynolds number. This is the first experimental evidence of the destabilizing effect of the shear-thinning behaviour in comparison with the Newtonian case: the critical Reynolds number is smaller, and the ratio between the critical wave celerity and the flow velocity at the free surface is larger.

International audience; We combine an effective diffusivity model with a numerical approach initially proposed by Meyer (1981) to simulate transient heat and mass transfer phenomena in a directionally solidifying Sn-Bi rod. This particularly efficient 1D numerical model is light enough to be used within the frame of optimization methods at reasonable numerical cost. This approach is tested against reference in situ measurements obtained under microgravity conditions during the Mephisto program. We simulate the final homogenization transient of several experimental runs with different pulling velocities. The solid/liquid interface temperature evolution with time is extracted from the simulations and compared with that obtained by Seebeck in line measurements. After optimization of the model the observed discrepancy between the simulated and measured data is less than 1.5%. This validates both the proposed very efficient 1D numerical approach and the consistency of the set of thermophysical parameters values for dilute Sn-Bi alloys.

International audience; Acoustic streaming designates the ability to drive quasisteady flows by acoustic propagation in dissipative fluids and results from an acoustohydrodynamics coupling. It is a noninvasive way of putting a fluid into motion using the volumetric acoustic force and can be used for different applications such as mixing purposes. We present an experimental investigation of a kind of square flow driven by acoustic streaming, with the use of beam reflections, in a water tank. Time-resolved experiments using particle image velocimetry have been performed to investigate the velocity field in the reference plane of the experiments for six powers: 0.5, 1, 2, 4, 6, and 8 W. The evolution of the flow regime from almost steady to strongly unsteady states is characterized using different tools: the plot of time-averaged and instantaneous velocity fields, the calculation of presence density maps for vortex positions and for the maximal velocity and vorticity crest lines, and the use of spatiotemporal maps of the waving observed on the jets created by acoustic streaming. A transition is observed between two regimes at moderate and high acoustic forcing.