In this paper, a numerical modeling of the impact, spreading, and eventually rebound of a viscoelastic droplet is reported. The numerical model is based on the volume of fluid (VOF) method coupled with the FENE-CR constitutive equations, and accounts for both the surface tension and the substrate wettability. The FENE-CR constitutive equations are used to model the polymer solution, while taking advantage of its rheological characterization. The comparison is performed between droplets of Newtonian solvent and a monodisperse polymer solution. The droplet impact on both hydrophilic and superhydrophobic substrate is analyzed through a detailed analysis of the spreading diameter evolution. It is found that while the droplet kinematic phase seems independent of the substrate and fluids properties, the recoiling phase is highly related to all of them. In addition, the model infers a critical polymer concentration above which the droplet rebound from a superhydrophobic substrate is suppressed. The simulation is of particular interest to ink-jet processing, and demonstrates the capability of the model to handle complex non-Newtonian droplet dynamics.

Accurate evaluation of thermo-fluid dynamic characteristics in tank are critically important for designing liquid hydrogen tank of small-scale hydrogen liquefier to minimize heat in-leak into the liquid and ullage. Due to the huge cost, most future liquid hydrogen storage tank designs will have to rely on predictive computational models for storage tank pressurization and heat-leak minimization. Thus, reliable on predictive numerical models to aid design process of tank design have been required. Therefore, in this study, in order to improve the storage efficiency of the small-scale hydrogen liquefier, a three-dimensional CFD model that can predict the boil-off rate and the thermo-fluid characteristics due to heat penetration has been developed. The prediction performance and accuracy of the CFD model is validated on the basis of the comparisons between its results and previous experimental data and a good agreement is obtained. To evaluate insulation performance of polyurethane foam according to thickness, the pressure change and thermo-fluid characteristics in a partially liquid hydrogen tank, subject to fixed ambient temperature and wind velocity, has been investigated numerically. Three different insulation thickness were considered to investigate not only interfacial characteristics between liquid and ullage but thermo-fluid dynamics in hydrogen storage tank for various heat-leak situation. The results indicate that as the insulation thickness becomes thinner and more heat penetration exists, natural convection by buoyancy develops stronger, so strong vortices are created near the interfacial zone by upward flow driven by buoyancy. In particular, as the amount of heat in-leak increased, wider and stronger heat accumulation was found in the top dished head area of the ullage area. In addition, it was confirmed that the CFD model developed in this study could well describe the phenomenon in which the temperature of the tank wall and the insulation mat increased at a faster rate due to faster heat transfer from the vapor and interface. In the future, the numerical model developed in this study will be used for optimizing insulation system of storage tank for small-scale hydrogen which is cost-effective and highly efficient.

In this paper, we investigate the accuracy and robustness of three classes of methods for solving two-phase incompressible flows on a staggered grid. Here, the unsteady two-phase flow equations are simulated by finite volumes and penalty methods using implicit and monolithic approaches (such as the augmented Lagrangian and the fully coupled methods), where all velocity components and pressure variables are solved simultaneously (as opposed to segregated methods). The interface tracking is performed with a Volume-of-Fluid (VOF) method, using the Piecewise Linear Interface Construction (PLIC) technique. The home code Fugu is used for implementing the various methods. Our target application is the simulation of two-phase flows at high density and viscosity ratios, which are known to be challenging to simulate. The resulting strategies of monolithic approaches will be proven to be considerably better suited for these two-phase cases, they also allow to use larger time step than segregated methods.