skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Effects of the inert phase on solidification near a triple-phase line
The local temperature solution near the triple-phase line of a solidifying front, its melt, and a surrounding inert phase was obtained analytically including all three phases and solidification kinetics. This analytical solution was validated using a three-phase numerical model of the horizontal ribbon growth of silicon and compared to a two-phase analysis that models the effect of the third phase (e.g. the gas) as an applied heat flux. Although the three-phase solutions have additional modes to represent the gas behavior, for many conditions the two-phase and three-phase models predicted consistent behaviors. However, introduction of a non-zero growth angle causes the gas phase heat fluxes to have strong gradients near the triple-phase line. Even with zero growth angle, there are conditions in which the two-phase and three-phase solutions are very different; one predicting infinite heat fluxes while the other predicts finite fluxes. This depended on the ratios of thermal conductivities, and the angle at which the solid-melt interface intersected the free surface. In particular, when the thermal conductivity of the inert phase was comparable to the melt or solid phases there were significant differences.  more » « less
Award ID(s):
2317674
PAR ID:
10528507
Author(s) / Creator(s):
;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Journal of Crystal Growth
Volume:
625
Issue:
C
ISSN:
0022-0248
Page Range / eLocation ID:
127438
Subject(s) / Keyword(s):
computer simulation heat transfer interfaces solidification triple-phase line growth from melt
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Three-phase (solid, melt, and gas) and two-phase (solid and melt) models of horizontal ribbon growth were compared to identify the significance of different gas effects. The boundary conditions at the melt–gas and solid–gas interfaces for two-phase simulations were obtained from decoupled simulations of the gas phase. The results showed that the gas shear stress strongly changes the flow and temperature fields and the position of the triple-phase line. Also, the gas pressure distribution determined the vertical position of the triple-phase line. In the absence of growth angle effects, the results of the two-phase model with specified convective heat transfer coefficient, shear stress, and pressure as boundary conditions along the gas phase interface closely matched that of the three-phase model. Even with non-zero growth angle effects, the two-phase model with all the boundary conditions agreed well with three-phase simulation results, despite increased deviations at higher pull speeds. Finally, the results indicated that gas-induced velocities are significant compared to the Marangoni and buoyancy velocities, which could lead to flow instabilities and the variations in solid shape as observed in HRG experiments. 
    more » « less
  2. Three-phase (solid, melt, and gas) and two-phase (solid and melt) models of horizontal ribbon growth were compared to identify the significance of different gas effects. The boundary conditions at the melt-gas and solid-gas interfaces for two-phase simulations were obtained from decoupled simulations of the gas phase. The results showed that the gas shear stress strongly changes the flow and temperature fields and the position of the triple-phase line. Also, the gas pressure distribution determined the vertical position of the triple-phase line. In the absence of growth angle effects, the results of the two-phase simulation with convective heat transfer coefficient, shear stress, and pressure specified closely matched that of the three-phase model. Even with non-zero growth angle effects the two-phase model with all the boundary conditions applied, agreed well with three-phase simulation results despite increased deviations at higher pull speeds. Finally, the results indicated that gas-induced velocities are significant compared to the Marangoni and buoyancy velocities, which could lead to flow instabilities and the variations in solid shape as observed in HRG experiments. 
    more » « less
  3. Studies in atomic-scale modeling of surface phase equilibria often focus on temperatures near zero Kelvin due to the challenges in calculating the free energy of surfaces at finite temperatures. The Bayesian-inference-based nested sampling (NS) algorithm allows for modeling phase equilibria at arbitrary temperatures by directly and efficiently calculating the partition function, whose relationship with free energy is well known. This work extends NS to calculate adsorbate phase diagrams, incorporating all relevant configurational contributions to the free energy. We apply NS to the adsorption of Lennard-Jones (LJ) gas particles on low-index and vicinal LJ solid surfaces and construct the canonical partition function from these recorded energies to calculate ensemble averages of thermodynamic properties, such as the constant-volume heat capacity and order parameters that characterize the structure of adsorbate phases. Key results include determining the nature of phase transitions of adsorbed LJ particles on flat and stepped LJ surfaces, which typically feature an enthalpy-driven condensation at higher temperatures and an entropy-driven reordering process at lower temperatures, and the effect of surface geometry on the presence of triple points in the phase diagrams. Overall, we demonstrate the ability and potential of NS for surface modeling. 
    more » « less
  4. none (Ed.)
    SUMMARY The geophysical detection of magma bodies and the estimation of the dimensions, physical properties and the volume fraction of each phase composing the magma is required to improve the forecasting of volcanic hazards and to understand transcrustal magmatism. We develop an analytical model to calculate P waves velocity in a three-phase magma consisting of crystals and gas bubbles suspended in a viscous melt. We apply our model to calculate the speed of sound as a function of the temperature in three magmas with different chemical compositions, representative of the diversity that is encountered in arc magmatism. The model employs the coupled phase theory that explicitly accounts for the exchanges of momentum and heat between the phases. We show that the speed of sound varies nonlinearly with the frequency of an acoustic perturbation between two theoretical bounds. The dispersion of the sound in a magma results from the exchange of heat between the melt and the dispersed phases that affects the magnitude of their thermal expansions. The lower bound of the sound speed occurs at low frequencies for which all the constituents can be considered in thermal equilibrium, whereas the upper bound occurs at high frequencies for which the exchange of heat between the phases may be neglected. The presence of gas in a magma produces a sharp decrease in the velocity of compressional waves and generates conditions in which the dispersion of the sound is significant at the frequencies usually considered in geophysics. Finally, we compare the estimates of our model with the ones from published relationships. Differences are largest at higher frequencies and are <10 per cent for typical magma. 
    more » « less
  5. The evaporation of droplets on surfaces is a ubiquitous phenomenon essential in nature and industrial applications ranging from thermal management of electronics to self-assembly-based fabrication. In this study, water droplet evaporation on a thin quartz substrate is analyzed using an unsteady two-step arbitrary Lagrangian-Eulerian (ALE) moving mesh model, wherein the evaporation process is simulated during the constant contact radius (CCR) and contact angle (CCA) modes. The numerical model considers mass transfer in the gas domain, flow in the liquid and gas domains, and heat transfer in the solid, liquid, and gas domains. Besides, the model also accounts for interfacial force balance, including thermocapillary stresses, to obtain the instantaneous droplet shape. Experiments involving droplet evaporation on unheated quartz substrates agree with model predictions of contact radius, contact angle, and droplet volume. Model results indicating temperature and velocity distribution across an evaporating water droplet show that the lowest temperatures are at the liquid-gas interface, and a single vortex exists for the predominant duration of the droplet's lifetime. The temperature of the unheated substrate is also significantly reduced due to evaporative cooling. The interfacial evaporation flux distribution, which depends on heat transfer across the droplet and advection in the surrounding medium, shows the highest values near the three-phase contact line. In addition, the model also predicts evaporation dynamics when the substrate is heated and exposed to different advection conditions. Generally, higher evaporation rates result from higher substrate heating and advection rates. However, substrate heating and advection in the surrounding gas have minimal effects on the relative durations of CCR and CCA modes for a given receding contact angle. Specifically, in this case, a 40× increase in substrate heating rate or 7.5× increase in gas velocity can only change these relative durations by 3%. This study also highlights the importance of surface wettability, which affects evaporation dynamics for all the conditions explored by the numerical model. 
    more » « less