A continuum model for momentum, heat and species transport in binary solid-liquid phase change systems—I. Model formulation

doi:10.1016/0017-9310(87)90094-9

International Journal of Heat and Mass Transfer

Volume 30, Issue 10, October 1987, Pages 2161-2170

https://doi.org/10.1016/0017-9310(87)90094-9 Get rights and content

Abstract

Semi-empirical laws and microscopic descriptions of transport behavior have been integrated with principles of classical mixture theory to obtain a set of continuum conservation equations for binary, solid-liquid phase change systems. For a restricted, yet frequently encountered, class of phase change systems, the continuum equations have been cast into forms amenable to clear physical interpretation and solution by conventional numerical procedures.

Résumé

Des lois semi-empiriques et des descriptions microscopiques du mécanisme de transfert ont été intégrées à la theorie classique des mélanges pour obtenir un système d'équations de bilan continu pour des systèmes à changement de phase solide-liquide. Pour une classe limitée mais fréquemment rencontrée, de systèmes les équations ont été mise sous des formes permettant une interprétation physique claire et une résolution par des procédures numériques conventionnelles.

Zusammenfassung

Halbempirische Gesetze und mikroskopische Beschreibungen des Transportverhaltens wurden mit den Prinzipien der klassischen Mischungstheorie verbunden. Es ergibt sich ein Satz von Erhaltungsgleichungen für binäre Fest-flüssig-Phasenwechsel-Systeme. Für eine bestimmte, häufig auftretende, Gruppe von Phasenwechselsystemen, wurden die Gleichungen in eine Form gebracht, die eine klare physikalische Interpretation ermöglicht und auβerdem mit herkömmlichen numerischen Verfahren lösbar ist.

Реферат

Для пoлyчeния cиcтeмы фeиoмeиoлoгичecкич ypaвнeний пepeнoca для бинapиыч cиcтeм твepдoe тeлo-жидкocть пpи фaзoвыч измeиeиияч выпoлнeиo oбьeдинeниe пoлyзмпиpичecкич зaкoнoв и микpocкoиичecкич oпиcaинй пpoцeccoв пepeнoca c зaкoиaми тeopии cмeceй. Для oгpaиичeииoгo, oдиaкo, чacтo вcтpeчaюшeгocя клacca cиcтeм пpи фaзoвыч пepeчoдaч пpивoдятcя ypaвнeиия цepeиoca в вндe, yдoбиoм для ич фнзичecкoй иитepпpeтaции и peщeиия c пoмoщью oбычиыч чиcлeииыч мeтoдoв.

References (21)

A.A. Lazaridis
A numerical solution of the multidimensional solidification (or melting) problem
Int. J. Heat Mass Transfer
(1970)
L.S. Yao
Natural convection effects in the continuous casting of a horizontal cylinder
Int. J. Heat Mass Transfer
(1984)
C.J. Ho et al.
Heat transfer during inward melting in a horizontal tube
Int. J. Heat Mass Transfer
(1984)
W.T. Sha et al.
Multidomain multiphase fluid mechanics
Int. J. Heat Mass Transfer
(1978)
D.S. Drumheller
The theoretical treatment of a porous solid using a mixture theory
Int. J. Solids Structures
(1978)
J.W. Nunziato
A multiphase mixture theory for fluid-particle flows
W.D. Bennon et al.
A continuum model for momentum, heat and species transport in binary solid-liquid phase change systems—II. Application to solidification in a rectangular cavity
Int. J. Heat Mass Transfer
(1987)
H.S. Carslaw et al.
Conduction of Heat in Solids
(1959)
M.N. Ozisik
Heat Conduction
(1980)
W.D. Murray et al.
Numerical and machine solutions of transient heat-conduction problems involving melting or freezing : Part I—method of analysis and sample solutions
J. Heat Transfer
(1959)

There are more references available in the full text version of this article.

Cited by (1007)

Material migration and damage characteristics of fused silica optical surface treated by CO<inf>2</inf> laser/short – pulse ultraviolet laser
2024, Optics and Laser Technology
Fused silica optics, as key components in high - power laser systems, are prone to be damaged under short - pulse ultraviolet (UV) laser irradiation. CO₂ laser is usually used to repair the laser damage and extend the service life of optics. However, the anti - laser damage characteristics of repaired substrates sometimes can’t meet the expectation. To facilitate better application of CO₂ laser repairing method, this work studied the CO₂ laser irradiating mechanism of fused silica in detail and investigated the damage characteristics of repaired optical surface through experiments. Firstly, finite element level - set (FELS) method is adopted to simulate CO₂ laser irradiating process. The simulation showed that CO₂ laser irradiation would cause molten material migration, vapor material deposition and affect the formation of residual stress. Various tests of morphology, residual stress, and photo - thermal absorption were conducted to study the changes of laser damage characteristics of repaired craters. The experimental results presented that molten materials migration and vapor materials re deposition occurred in CO₂ laser irradiation, which would promote the formation of residual stress, cause stronger laser absorption, or even induce laser damage. Through the experiment, accuracy of simulation was verified further. In this work, the CO₂ laser repairing mechanism of fused silica and laser damage characteristics of repaired optics are clarified through FELS simulation and experiments. Relative results can provide a certain technical support for the performance improvement of repaired fused silica optics.
Effect of withdrawal rate on freckle formation of large-size directional solidified nickel-based superalloy blades
2024, Journal of Materials Research and Technology
The freckle formation and distribution in large-size directional solidified nickel-based superalloy blades at different withdrawal rates were investigated. The microstructure of freckles was analyzed by optical microscope (OM), and scanning electron microscope (SEM). The crystallographic orientation in the freckle regions was characterized by electron backscattered diffraction (EBSD). The solute segregation patterns on the surface of freckle regions were obtained by X-ray fluorescence (XRF). Numerical simulations for the temperature field were performed to analyze the mushy zone evolution during directional solidification. The experimental results indicated that the root region is prone to form freckles. The freckle chains were only distributed on the outward surface of the root region at low and middle withdrawal rates, and the inward side was free of freckles. When a high withdrawal rate was applied, a few freckles occurred on the inward side of the root region. The number and length of freckle chains increased as the withdrawal rate increased. The XRF results indicated that channel segregation formed on the surface of the root region. The segregation channels predicted by the macrosegregation simulation agreed well with the experimental observation. The concave or inclined solidification interface increased the propensity for freckle formation.
Multi-scale and multi-physics simulation of central segregation in an equiaxed dendritic mushy zone during continuous casting of steel
2024, Materialia
Centreline segregation during continuous casting of steel is highly detrimental to the end-product mechanical properties. In this study, a 3D multi-scale and multi-physics model is developed and validated that directly predicts alloy segregation occurring in peritectic steel grades having centreline equiaxed grain morphology. The model consists of four models: (1) a 1D macroscale heat transfer and solidification model, (2) a 3D mesoscale dendritic solidification model, (3) a 3D mesoscale dendritic fluid flow model, and (4) a 3D mesoscale solute transport and redistribution model. The use of a mesoscale domain where the solid and liquid phases are separately discretized allows for consideration of physical phenomena affecting segregation that are classically hidden by the averaging procedures of fully continuum approaches. Thus the new model is able to account for the multiple phenomena occurring during continuous casting while also directly considering the stochastic effects resulting from grains having randomly-placed nuclei as well as interactions between the liquid and solid phases. The results demonstrate that solute partitioning combined with intra-dendritic fluid flow leads eventually to liquid channels enriched with solute that result in centreline segregation. The predicted composition in these discrete liquid channels shows excellent agreement with a Mn centreline segregation profile experimentally-measured via microscopic X-ray fluorescence of a commercially-cast steel slab. Finally, the effects of alloy composition, and soft reduction on the degree of centreline segregation are examined.
A multiphase modeling for investigating temperature history, flow field and solidification microstructure evolution of FeCoNiCrTi coating by laser cladding
2024, Optics and Laser Technology
In this work, a multiphase numerical model was established to investigate the non-isothermal flow and microstructure evolution during laser cladding. The model considered the beam-powder interaction, temperature dependent material properties, the effect of buoyancy and surface tension on convection. The accuracy of the model was verified comparing the simulated coating geometric morphology with the experimental results. Simulation indicated that the temperature and flow of the molten pool were interactive and mutually causal. On the one hand, the temperature field affected the tension gradient and the liquid flow direction. On the other hand, the velocity field significantly affected the heat transfer in the molten pool. Specifically, the maximum fluid velocity was 9.58 mm/s and the Peclet number (Pe_T) value was less than 5 at the initial stage, indicating that conduction dominated in heat transfer. Furthermore, the Pe_T value was greater than 200 when the maximum velocity reached 344 mm/s. The heat transfer in the molten pool was mainly thermal convection. In addition, the distribution diagram of solidification parameters (G and R) was calculated to quantitatively establish the internal correlation between solidification parameters and microstructure combing with experimental observation, and revealed the columnar-to-equiaxed transition (CET). Eventually, the microhardness results showed that HEA coating could improve the Ti6Al4V surface microhardness. Simultaneously, the CET behavior significantly affected the microhardness of the coating.
Anomalously strong effects of plume contraction and material redeposition in nanosecond pulsed laser vaporization
2023, International Journal of Heat and Mass Transfer
Long-term dynamics of vapor plumes induced by irradiation of a copper target in argon background gas by nanosecond laser pulses is studied numerically based on a one-dimensional hybrid computational model. The model includes a thermal model of the irradiated target and a kinetic model of the gaseous plume flow. The latter is implemented in the form of the direct simulation Monte Carlo method. The simulations, performed at moderate laser fluences in the range from 1.5 to 4 Jcm⁻² and background gas pressures between 0.01 and 1 bar, show that the long-term plume dynamics can be divided into three major stages: initial inertial expansion, plume contraction, and subsequent diffusive expansion. The simulations unexpectedly predict extremely strong effects of plume contraction, when the plume size can exhibit a two-fold reduction, and delayed material redeposition back to the irradiated surface, when $\sim 90$ % of vaporized material returns to the irradiated surface. Both plume contraction and delayed redeposition originate from the strong backward flow induced by the internal (secondary) shock wave that propagates from the mixing layer towards the irradiated surface. The major part of vaporized material condenses at the surface with a long time delay with respect to the laser pulse after the plume stopping and before the onset of diffusive expansion. As a result, only a marginal part of vapor, which can be as small as $\sim 5$ %, is retained in the plume by the beginning of the diffusive expansion. The obtained simulation results suggest that the plume contraction and delayed material redeposition are common phenomena for nanosecond laser vaporization and may affect the efficiency and quality of laser surface modification in the ablation regime as well.
Domination and effect of multi-parameters in direct chill casting based on establishment of thermo model by numerical simulation and experiment
2023, Journal of Materials Research and Technology
In this paper, a mathematical model validated by in-situ temperature measurement data of the direct chill (DC) casting process has been established in order to study the domination and effect of multi-parameters including secondary cooling water flow rate (Q), casting speed (V), condition of primary cooling, and spraying water temperature (T_w) on DC casting index of sump depth (d_sump), maximum and minimum temperature gradient (T_g,max and T_g,min) at mushy zone and shape of solidus line. Further, the height of impingement zone ( $H_{i p}$ ) is systematic quantitative analysis by a novel designed test ( $H_{i p} = 1.08865 Q ` + 0.06272 V - 4.75597$ ). The calculated temperature using the numerical simulation exhibits a good match with the actual temperature measurement data. The results of multi-factors on index by the designed orthogonal test indicate that V dominates the d_sump ( $d_{s u m p} = 0.033 + 27.74 V$ ) and T_g,min ( $T_{g, \min} = 2569.15 - 766342.20 V$ ) at mushy zone, Q dominates the T_g,max ( $T_{g, \max} = 9195.12 + 2344770 Q - 27.83 T_{w}$ ) by analysing the range (R), variance (S) and significance (F). Likewise, by discussing the combination of R, S, F and microstructure observation, the priority influencing the shape of solidus line at center part was V. The Q dominates the shape of solidus line at the position between the center and edge part.

View all citing articles on Scopus

View full text

Abstract

Résumé

Zusammenfassung

Реферат

Int. J. Heat Mass Transfer

Int. J. Heat Mass Transfer

Int. J. Heat Mass Transfer

Int. J. Heat Mass Transfer

Int. J. Solids Structures

Int. J. Heat Mass Transfer

Conduction of Heat in Solids

Heat Conduction

Numerical and machine solutions of transient heat-conduction problems involving melting or freezing : Part I—method of analysis and sample solutions

J. Heat Transfer