Minh Do-Quang

The splash of a solid sphere impacting on a liquid surface: Numerical simulation of the influence of wetting

The impact of a solid sphere on a liquid surface has challenged researchers for centuries and remains of interest today. Recently, Duez [Nat. Phys. 3, 180 (2007)] published experimental results of ...

Modeling of dynamic wetting far from equilibrium

In this paper we present simulations of dynamic wetting far from equilibrium based on phase field theory. In direct simulations of recent experiments [J. C. Bird, S. Mandre, and H. A. Stone, Phys. ...

Phase-field simulation of weld solidification microstructure in an Al-Cu alloy

Since the mechanical properties and the integrity of the weld metal depend on the solidification behaviour and the resulting microstructural characteristics, understanding weld pool solidification is of importance to engineers and scientists. Thermal and fluid flow conditions affect the weld pool geometry and solidification parameters. During solidification of the weld pool, a columnar grain structure develops in the weld metal. Prediction of the formation of the microstructure during welding may be an important and supporting factor for technology optimization. Nowadays, increasing computing power allows direct simulations of the dendritic and cell morphology of columnar grains in the molten zone for specific temperature conditions. In this study, the solidification microstructures of the weld pool at different locations along the fusion boundary are simulated during gas tungsten arc welding of Al–3wt%Cu alloy using the phase-field model for the directional solidification of dilute binary alloys. A macroscopic heat transfer and fluid flow model was developed to assess the solidification parameters, notably the temperature gradient and solidification growth rate. The effect of the welding speed is investigated. Computer simulations of the solidification conditions and the formation of a cellular morphology during the directional solidification in gas tungsten arc welding are described. Moreover, the simulation results are compared with existing theoretical models and experimental findings.

Efficient preconditioners for large scale binary Cahn–Hilliard models

Abstract In this work we consider preconditioned iterative solution methods for numerical simulations of multiphase flow problems, modelled by the Cahn-Hilliard equation. We focus on diphasic flows and the construction and efficiency of a preconditioner for the algebraic systems arising from finite element discretizations in space and the θ-method in time. The preconditioner utilises to a full extent the algebraic structure of the underlying matrices and exhibits optimal convergence and computational complexity properties. Various numerical experiments, including large scale examples, are presented as well as performance comparisons with other solution methods.

Locomotion by tangential deformation in a polymeric fluid

In several biologically relevant situations, cell locomotion occurs in polymeric fluids with Weissenberg number larger than 1. Here we present results of three-dimensional numerical simulations for the steady locomotion of a self-propelled body in a model polymeric (Giesekus) fluid at low Reynolds number. Locomotion is driven by steady tangential deformation at the surface of the body (the so-called squirming motion). In the case of a spherical squirmer, we show that the swimming velocity is systematically less than that in a Newtonian fluid, with a minimum occurring for Weissenberg numbers of order 1. The rate of work done by the swimmer always goes up compared to that occurring in the Newtonian solvent alone but is always lower than the power necessary to swim in a Newtonian fluid with the same viscosity. The swimming efficiency, defined as the ratio between the rate of work necessary to pull the body at the swimming speed in the same fluid and the rate of work done by swimming, is found to always be increased in a polymeric fluid. Further analysis reveals that polymeric stresses break the Newtonian front-back symmetry in the flow profile around the body. In particular, a strong negative elastic wake is present behind the swimmer, which correlates with strong polymer stretching, and its intensity increases with Weissenberg number and viscosity contrasts. The velocity induced by the squirmer is found to decay in space faster than in a Newtonian flow, with a strong dependence on the polymer relaxation time and viscosity. Our computational results are also extended to prolate spheroidal swimmers and smaller polymer stretching are obtained for slender shapes compared to bluff swimmers. The swimmer with an aspect ratio of two is found to be the most hydrodynamically efficient.

Simulation of free dendritic crystal growth in a gravity environment

In this paper we simulate the evolution and free particle motion of an individual nucleus that grows into a dendritic crystal. The melt flow and the convective heat transfer around the crystal are simulated as they settle due to gravity. There is an intricate coupling between the settling and the evolution of the crystal. The relative flow induced by the settling enhances the growth at the downward facing parts, which in its turn affects the subsequent settling motion.Simulations have been done in two dimensions using a semi-sharp phase-field model. The flow was constrained to a rigid body motion by using Lagrange multipliers inside the solidified part. The model was formulated using two different meshes. One is a fixed background mesh, which covers the whole domain. The other is an adaptive mesh, where the node points are also translated and rotated with the movement of the solid particle. In the latter, the dendritic growth is simulated by the semi-sharp phase-field method.

Numerical simulation of the coupling problems of a solid sphere impacting on a liquid free surface

This paper presents a model, using a phase-field method, that is able to simulate the motion of a solid sphere impacting on a liquid surface, including the effects of capillary and hydrodynamic forces. The basic phenomena that were the subject of our research effort are the small scale mechanism such as the wetting property of the solid surface which control the large scale phenomena of the interaction. The coupled problem during the impact will be formulated by the inclusion of the surface energies of the solid surface in the formulation, which gives a reliable prediction of the motion of solid objects in/on/out of a liquid surface and the hydrodynamic behaviours at small scales when the inertia of fluid is less important than its surface tension. Numerical results at different surface wettabilities and impact conditions will be presented and compared with the experiments of Duez el al. [C. Duez, C. Ybert, C. Clanet, L. Bocquet, Nat. Phys. 3 (2007) 180-183] and Lee and Kim [D. Lee, H. Kim, Langmuir 24 (1) (2008) 142].

Element-by-Element schur complement approximations for general nonsymmetric matrices of two-by-two block form

We consider element-by-element Schur complement approximations for indefinite and general nonsymmetric matrices of two-by-two block form, as arising in finite element discretized systems of PDEs The paper provides some analysis of the so-obtained approximation and attempts to quantify the quality of the underlying two-by-two matrix splitting in a way similar to that used for symmetric positive definite matrices The quality of the approximation is illustrated numerically.

Early Onset of Nucleate Boiling on Gas-covered Biphilic Surfaces

For phase-change cooling schemes for electronics, quick activation of nucleate boiling helps safeguard the electronics components from thermal shocks associated with undesired surface superheating at boiling incipience, which is of great importance to the long-term system stability and reliability. Previous experimental studies show that bubble nucleation can occur surprisingly early on mixed-wettability surfaces. In this paper, we report unambiguous evidence that such unusual bubble generation at extremely low temperatures—even below the boiling point—is induced by a significant presence of incondensable gas retained by the hydrophobic surface, which exhibits exceptional stability even surviving extensive boiling deaeration. By means of high-speed imaging, it is revealed that the consequently gassy boiling leads to unique bubble behaviour that stands in sharp contrast with that of pure vapour bubbles. Such findings agree qualitatively well with numerical simulations based on a diffuse-interface method. Moreover, the simulations further demonstrate strong thermocapillary flows accompanying growing bubbles with considerable gas contents, which is associated with heat transfer enhancement on the biphilic surface in the low-superheat region.

Thermohydrodynamics of boiling in binary compressible fluids

We numerically study the thermohydrodynamics of boiling for a CO(2) + ethanol mixture on lyophilic and lyophobic surfaces in both closed and open systems, based on a diffuse interface model for a two-component system. The corresponding wetting boundary conditions for an isothermal system are proposed and verified in this paper. New phenomena due to the addition of another component, mainly the preferential evaporation of the more volatile component, are observed. In the open system and the closed system, the physical process shows very different characteristics. In the open system, except for the movement of the contact line, the qualitative features are rather similar for lyophobic and lyophilic surfaces. In the closed system, the vortices that are observed on a lyophobic surface are not seen on a lyophilic surface. More sophisticated wetting boundary conditions for nonisothermal, two-component systems might need to be further developed, taking into account the variations of density, temperature, and surface tension near the wall, while numerical results show that the boundary conditions proposed here also work well even in boiling, where the temperature is nonuniform.