An official website of the United States government
Abstract We present a phase-field crystal model for solidification that accounts for thermal transport and a temperature-dependent lattice parameter. Elasticity effects are characterized through the continuous elastic field computed from the microscopic density field. We showcase the model capabilities via selected numerical investigations which focus on the prototypical growth of two-dimensional crystals from the melt, resulting in faceted shapes and dendrites. This work sets the grounds for a comprehensive mesoscale model of solidification including thermal expansion.
more »
« less
A non-isothermal phase-field crystal model with lattice expansion: analysis and benchmarks
Abstract We introduce a non-isothermal phase-field crystal model including heat flux and thermal expansion of the crystal lattice. The fundamental thermodynamic relation between internal energy and entropy, as well as entropy production, is derived analytically and further verified by numerical benchmark simulations. Furthermore, we examine how the different model parameters control density and temperature evolution during dendritic solidification through extensive parameter studies. Finally, we extend our framework to the modeling of open systems considering external mass and heat fluxes. This work sets the ground for a comprehensive mesoscale model of non-isothermal solidification including thermal expansion within an entropy-producing framework, and provides a benchmark for further meso- to macroscopic modeling of solidification.
more »
« less
- Award ID(s):
- 2309547
- PAR ID:
- 10642159
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Modelling and Simulation in Materials Science and Engineering
- Volume:
- 33
- Issue:
- 2
- ISSN:
- 0965-0393
- Page Range / eLocation ID:
- 025007
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract A bulk gadolinium (Gd) single crystal exhibits virtually zero remnant magnetization, a common trait among soft uniaxial ferromagnets. This characteristic is reflected in our magnetometry data showing virtually hysteresis free isothermal magnetization loops with large saturation magnetization. The absence of hysteresis allows to model the measured easy axis magnetization as a function of temperature and applied magnetic field, rather than a relation, which permits the application of Maxwell relations from equilibrium thermodynamics. Demagnetization effects broaden the isothermal first-order transition from negative to positive magnetization. By analyzing magnetization data within the coexistence regime, we deduce the isothermal entropy change and the field-induced heat capacity change. Comparing the numerically inferred heat capacity with relaxation calorimetric data confirms the applicability of the Maxwell relation. Analysis of the entropy in the mixed phase region suggests the presence of hitherto unresolved nanoscale magnetic structures in the demagnetized state of Gd. To support this prediction, Monte Carlo simulations of a 3D Ising model with dipolar interactions are performed. Analyzing the cluster size statistics and magnetization from the model provides strong qualitative support of our analytic approach.more » « less
-
null (Ed.)Thermoelastic behavior of transition metal boride Os2B3 was studied under quasi-hydrostatic and isothermal conditions in a Paris-Edinburgh cell to 5.4 GPa and 1273 K. In-situ Energy Dispersive X-ray diffraction was used to determine interplanar spacings of the hexagonal crystal structure and thus the volume and axial compression. P-V-T data were fitted to a 3rd Order Birch-Murnaghan equation of state with a temperature modification to determine thermal elastic constants. The bulk modulus was shown to be K0 = 402 ± 21 GPa when the first pressure derivative was held to K0’ = 4.0 from the room temperature P-V curve. Under a quadratic fit α=α_0+α_1 T-α_2 T^(-2), the thermal expansion coefficients were determined to be α_0=1.862×10^(-5) K-1, α_1=0.841×10^(-9) K-2, and α_2=-0.525 K. Density functional theory (DFT) with the quasi-harmonic approximation (QHA) were further employed to study Os2B3, including its P-V-T curves, phonon spectra, bulk modulus, specific heat, thermal expansion, and the Grüneisen parameter. A good agreement between the first-principle theory and experimental observations was achieved, highlighting the success of the Armiento-Mattsson 2005 generalized gradient approximation functional employed in this study and QHA for describing thermodynamic properties of Os2B3.more » « less
-
Refractory high-entropy alloys (RHEAs) are strong candidates for use in high-temperature engineering applications. As such, the thermodynamic properties as a function of temperature for a variety of RHEA systems need to be studied. In the present work, thermodynamic quantities such as entropy, enthalpy, heat capacity at constant volume, and linear thermal expansion are calculated for three quaternary and three quinary single-phase, BCC RHEAs: AlMoNbV, NbTaTiV, NbTaTiZr, AlNbTaTiV, HfNbTaTiZr, and MoNbTaVW. First-principle calculations based on density functional theory are used for the calculations, and special quasirandom structures (SQSs) are used to represent the random solid solution nature of the RHEAs. A code for the finite temperature thermodynamic properties using the Debye-Grüneisen model is written and employed. For the first time, the finite temperature thermodynamic properties of all 24 atomic configuration permutations of a quaternary RHEA are calculated. At most, 1.7% difference is found between the resulting properties as a function of atomic configuration, indicating that the atomic configuration of the SQS has little effect on the calculated thermodynamic properties. The behavior of thermodynamic properties among the RHEAs studied is discussed based on valence electron concentration and atomic size. Among the quaternary RHEAs studied, namely AlMoNbV, NbTaTiZr, and NbTaTiV, it is found that the presence of Zr contributes to higher entropy. Additionally, at lower temperatures, Zr contributes to higher heat capacity and thermal expansion compared to the alloys without Zr, possibly due to its valence electron concentration. At higher temperatures, Al contributes to higher heat capacity and thermal expansion, possibly due its ductility. Among the quinary systems, the presence of Mo, W, and/or V causes the RHEA to have a lower thermal expansion than the other systems studied. Finally, when comparing the systems with the NbTaTi core, the addition of Al increases thermal expansion, while the removal of Zr lowers the thermal expansion.more » « less
-
Compressed air energy storage (CAES) is a low-cost, long-duration storage option under research development. Several studies suggest that near-isothermal compression may be achieved by injecting water droplets into the air during the process to increase the overall efficiency. However, little is known about the thermal-fluid mechanisms and the controlling nondimensional parameters of the expansion process, which has previously been assumed to mirror the compression process. Furthermore, the isothermal round-trip efficiency and the impact of spray-based CAES have not been investigated. This study uses a validated 1-D model for compression and expansion with spray injection to complete a parametric analysis to analyze the thermal-fluid time-dependent physics and resulting roundtrip isothermal efficiency of a CAES system. Comparing the results for compression and expansion simulations, compression is found to have a higher isothermal efficiency than expansion for the same set up. The polytropic index for both compression and expansion tends to decrease and approach the ideal isothermal limit as nondimensional mass loading increases and as nondimensional Crowe number (ratio of thermal response time to domain time) decreases. As such, the highest efficiency designs are those with slow compression speeds and high spray flow rates to achieve high mass loading and those with small droplets to achieve low Crowe numbers—as long as spray work is neglected. If spray work is included, the optimum spray conditions shift to those with lager drop sizes. For example, roundtrip isothermal efficiency peaks around 95 % at a mass loading of 14 and at Crowe numbers <0.1 with a pressure ratio of ten. The results indicate that near-isothermal CAES compression and expansion is possible but that spray work should be included for significant mass loadings (e.g. greater than unity). Further investigation is recommended to consider effects of multi-dimensionality, turbulence, wall-interactions, and droplet dynamics.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1088/1361-651X/ada784
