Search NASA⌕ Search

SEARCH · Search NASA

Results for “Melting point”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experimental characterization of a triply periodic minimal surface PCM-to-air thermal energy storage device

Here, the expansion of energy production has intensified research into thermal energy storage (TES) to manage variability and improve system efficiency. Heat exchanger design is critical to achieving high power densities in TES systems. This study presents a high-surface-area phase change material (PCM)-to-air heat exchanger fabricated via resin-based stereolithography with a novel gyroid-based geometry tailored for enhanced performance. Material properties of the commercial PCM and resin were characterized using analytical techniques. A controlled air loop was used to evaluate heat transfer and pressure drop at various flow rates and inlet temperatures. Results show a strong dependence on the inlet temperature difference (ΔT) relative to the PCM melting point. During charging at the highest flow rate, increasing the inlet air temperature from ΔT = 5 °C to 20 °C above the melting point increased the average heat transfer rate by 187%. During discharging, decreasing the inlet temperature by the same ΔT below the melting point led to a 232% increase. Notably, the high-surface-area design enabled nearly symmetric charging and discharging behavior, a novel result for PCM-based TES systems which are often restricted by natural convection and other effects. The overall heat transfer coefficient was calculated and compared to values from standard design correlations. The maximum thermal effectiveness reached 95% at ΔT = 20 °C and a moderate flow rate of 34 m 3 /h. Peak coefficients of performance (COP) of 7 during charging and 6.4 during discharging were observed at ΔT = 20 °C and a low flow rate of 20 m 3 /h. These results demonstrate the viability of additively manufactured geometries for advanced TES applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Predicting the influence of particle size on the glass transition temperature and viscosity of secondary organic material

Abstract Atmospheric aerosols can assume liquid, amorphous semi-solid or glassy, and crystalline phase states. Particle phase state plays a critical role in understanding and predicting aerosol impacts on human health, visibility, cloud formation, and climate. Melting point depression increases with decreasing particle diameter and is predicted by the Gibbs–Thompson relationship. This work reviews existing data on the melting point depression to constrain a simple parameterization of the process. The parameter $$\xi $$ ξ describes the degree to which particle size lowers the melting point and is found to vary between 300 and 1800 K nm for a wide range of particle compositions. The parameterization is used together with existing frameworks for modeling the temperature and RH dependence of viscosity to predict the influence of particle size on the glass transition temperature and viscosity of secondary organic aerosol formed from the oxidation of $$\alpha $$ α -pinene. Literature data are broadly consistent with the predictions. The model predicts a sharp decrease in viscosity for particles less than 100 nm in diameter. It is computationally efficient and suitable for inclusion in models to evaluate the potential influence of the phase change on atmospheric processes. New experimental data of the size-dependence of particle viscosity for atmospheric aerosol mimics are needed to thoroughly validate the predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗
DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Triaxial Constraint and Tensile Strength Enhancement in Brazed Joints

A brazed joint consists of a low-melting point and thin interlayer sandwiched between the high-melting-point base materials, in which the interlayer strength is typically lower than that of the base material. When this butt-joined composite is loaded uniaxially in the direction perpendicular to the plane of the brazing layer, the tensile strength is found to be much higher than that of the braze. This seems to violate the iso-stress condition in such a butt-joint serial configuration. Furthermore, the stress triaxiality has been usually ascribed, but without a quantitative rationalization, as being responsible for this tensile strength enhancement. Here a complete finite element simulation has been conducted to study the dependence of triaxiality and strength enhancement on geometric and material parameters. Two asymptotic limit solutions (based on Bridgman and Xia–Shih solutions, respectively) have been identified to understand the simulation results. The critical role of void evolution has been revealed when making a quantitative comparison to available experiments. In addition, ductility of the brazed joint, which has not been fully addressed in literature, is investigated by the Gurson–Tvergaard–Needleman model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of Heaters on Molten Salt Reactors Dynamics

Recently, there has been renewed interest in the Molten Salt Reactor (MSR) concept. This interest is mainly due to its advantages over Light Water Reactors (LWRs). Among these advantages are the flexibility in the fuel choice and the possibility to burn actinides [1]. Currently, two molten research reactors are under development. A 2 MW molten salt experimental reactor has been constructed and is planned for operation in China. A 1 MW Molten Salt Research Reactor is under construction in the United States [2]. Thus, Research and Development (R&D) programs are needed for the MSR technology. Reactor dynamics studies are crucial safety evaluation studies. The development of reactor dynamics tools for MSRs started during the Molten Salt Reactor Program (MSRP) at ORNL, focusing on the Molten Salt Reactor Experiment (MSRE) [3]. Although various dynamic simulation tools are available, R&D efforts to develop specialized tools for MSRs are still ongoing [4]. The primary motivation for this interest is the goal of commercializing MSRs. Developing a dynamics tool for MSRs is crucial for analyzing their safety and supporting their demonstration efforts. A crucial safety issue in MSRs is ensuring that the fuel salt in the primary loop remains molten. This concern becomes particularly important when the reactor operates at lower power levels. Most chloride and fluoride salts used as fuel have a melting point of approximately 450 °C [5]. Therefore, maintaining fuel salt temperatures above this melting point is essential to prevent it from freezing. Freezing of the fuel salt can lead to volume expansion, potentially causing damage to reactor components [6]. Moreover, it can create local blockages within fuel salt channels, reducing cooling efficiency and resulting in hot spots [7]. Various methods can prevent fuel salt from freezing at low power levels. These include reducing the fuel mass flow rate in the primary loop, decreasing heat exchange between the primary and secondary loops in the reactor, and using electric fuel salt heaters. Reducing the mass flow rate of fuel salt in the primary loop can increase the fuel salt transit time in the core, enhancing heat production and increasing temperature. However, it can potentially lead to challenges. These include increasing the potential for fuel salt deposition on the reactor channel walls [8], affecting the overall heat transfer processes. Thus, alternative methods should be used [9]. This research investigates the effects of fuel salt heaters on the dynamics and stability during MSR operation. The investigation includes studying the reactor stability with the fuel salt heaters on and off. In addition, the reactor response to transients is compared with and without the operation of fuel salt heaters. These transients include reactivity insertion, primary and secondary pump failure, and heat sink temperature increase. The results of the reactivity insertion transient are presented in the summary.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Superheating gold beyond the predicted entropy catastrophe threshold

In their landmark study, Fecht and Johnson unveiled a phenomenon that they termed the ‘entropy catastrophe’, a critical point where the entropy of superheated crystals equates to that of their liquid counterparts. This point marks the uppermost stability boundary for solids at temperatures typically around three times their melting point. Despite the theoretical prediction of this ultimate stability threshold, its practical exploration has been prevented by numerous intermediate destabilizing events, colloquially known as a hierarchy of catastrophes, which occur at far lower temperatures. Here we experimentally test this limit under ultrafast heating conditions, directly tracking the lattice temperature by using high-resolution inelastic X-ray scattering. Our gold samples are heated to temperatures over 14 times their melting point while retaining their crystalline structure, far surpassing the predicted threshold and suggesting a substantially higher or potentially no limit for superheating. We point to the inability of our samples to expand on these very short timescales as an important difference from previous estimates. These observations provide insights into the dynamics of melting under extreme conditions.

Laser-produced plasmas↗
DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2022

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. Here, we used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

58 GEOSCIENCES↗
DOE OSTI2020

Ceramic encapsulated metal phase change material for high temperature thermal energy storage

Thermal energy storage (TES) is a broad-based technology for reducing CO₂ emissions and advancing concentrating solar, fossil, and nuclear power through improvements in efficiency and economics. Phase change materials (PCMs) are of interest as TES media because of their ability to store large amounts of heat in relatively small volumes. The volume expansion during a phase change, typically between a solid and liquid, can cause breakage of protective coatings. This paper reports on the fabrication of a ceramic encapsulated metal (CEM) high temperature TES technology using a rotary calcining furnace and a fluidized bed chemical vapor deposition coating technique. Aluminum beads were chosen as the PCM because Al has a high melting point (660 °C), low cost, high heat of fusion, and an ability to form a thin, strong alumina layer capable of supporting the Al melt for subsequent processing. Quite remarkably, this study shows that 1 mm diameter Al can be fluidized up to at least 1500 °C in an appropriate atmosphere while maintaining a spheroid geometry. This allowed for producing a first of a kind CEM whereby Al particles were encapsulated in pyro-carbon (PyC) and high purity, dense chemical vapor deposited SiC. The CEM with a PyC only coating was exposed to thermal cycling to test the performance with a differential scanning calorimeter; the melting point and latent heat were measured to be 648.4 ± 2.8 °C and 293.3 ± 6.2 J/g respectively. It was demonstrated that the CEM design is possible to produce, laying the foundation for manufacturing of high temperature, tunable, TES media.

25 ENERGY STORAGE↗
DOE OSTI2023

Hot Isostatic Pressing Control of Tungsten-Based Composites

Metal-oxide composites are commonly used in high temperature environments for their thermal stability and high melting points. Commonly employed with refractory oxides or carbides such as ZrC and HfC, these materials may be improved with the use of a low density, high melting point ceramic such as CeO2. In this work, the consolidation of W-CeO2 metal matrix composites in the high CeO2 concentration regime is explored. The CeO2 concentrations of 50, 33, and 25 wt.%, the CeO2 particle size from nanometer to micrometer, and various hot isostatic pressing temperatures are investigated. Decreasing the CeO2 concentration is observed to increase the composite density and increase the Vickers hardness. The CeO2 oxidation state is observed to be a combination of Ce3+ and Ce4+, which is hypothesized to contribute to the porosity of the composites. The hardness of the metal-oxide composite can be improved more than 2.5 times compared to pure W processed by the same route. This work offers processing guidelines for further consolation of oxide-doped W composites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mg(BH 4 ) 2 -Based Hybrid Metal–Organic Borohydride System Exhibiting Enhanced Chemical Stability in Melt

The formation of a chemically stable and thermally cyclable Mg(BH 4 ) 2 melt is proven to be possible through the utilization of an organic borohydride salt additive. While extensive exploration of additives for lowering the melting point and modifying the chemical stability of Mg(BH 4 ) 2 has been reported, this is the first study to use the organic borohydride salt, tetramethylammonium borohydride (TMAB), to modify the melting behavior of Mg(BH 4 ) 2 . Here, examination of a 5:1 molar mixture of Mg(BH 4 ) 2 and TMAB revealed a reversible melt between 180 and 195 °C, which was reproduced for five thermal cycles. The mixture melt exhibited an enhanced chemical stability compared to melts of the individual Mg(BH 4 ) 2 and TMAB species. It was observed that between room temperature and 250 °C (over 50 °C above the melting point), the mixture releases <0.1 wt % mass, consisting primarily of H 2 . The mixture also exhibits greatly reduced evolution of volatile boron containing compounds compared to either Mg(BH 4 ) 2 or TMAB. The use of TMAB to chemically stabilize a Mg(BH 4 ) 2 -rich melt demonstrated in this work represents an exciting pathway to modification of Mg(BH 4 ) 2 relevant to both hydrogen-storage and magnesium battery fields.

08 HYDROGEN↗
DOE OSTI2024

IC w22_phadiagractox Highlight: Composition-dependent melting temperature of a compound [Poster]

Figure 1 displays, ambient melting points of the Au-Cu system as a function of Cu atomic percent: present model vs. experimental data from four different sources. Figure 2 displays, ambient melting points of the U-O system as a function of O atomic percent: present model vs. both experiment and quantum molecular dynamics (QMD) simulations on the stoichiometric U-O compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Ambient melting behavior of stoichiometric uranium oxides

As UO 2 is easily oxidized during the nuclear fuel cycle it is important to have a detailed understanding of the structures and properties of the oxidation products. Experimental work over the years has revealed many stable uranium oxides including UO 2 , U 4 O 9 (UO 2.25 ), U 3 O 7 (UO 2.33 ), U 2 O 5 (UO 2.5 ), U 3 O 8 (UO 2.67 ), and UO 3 , all with a number of different polymorphs. These oxides are broadly split into two categories, fluorite-based structures with stoichiometries in the range of UO 2 to UO 2.5 and less dense layered-type structures with stoichiometries in the range of UO 2.5 to UO 3 . While UO 2 is well characterized, both experimentally and computationally, there is a paucity of data concerning higher stoichiometry oxides in the literature. In this work we determine the ambient melting points of all the six stoichiometric uranium oxides listed above and compare them to the available experimental and/or theoretical data. We demonstrate that a family of the six ambient melting points map out a solid-liquid transition boundary consistent with the high-temperature portion of the phase diagram of uranium-oxygen system suggested by Babelot et al.

Burakovsky, Leonid↗
DOE OSTI2020

Solid-solid phase transition via the liquid in a Pd 43 Cu 27 Ni 10 P 20 bulk metallic glass under conventional conditions

The microscopic mechanism of solid-solid phase transitions is a long-standing fundamental issue in materials science. Here we directly report the experimental evidence of the existence of an intermediate liquid during the solid-solid phase transition in a bulk Pd 43 Cu 27 Ni 10 P 20 alloy by in situ high temperature high energy X-ray diffraction, combined with the differential scanning calorimetry at slow heating rates down to 5 K/min. We elucidate that this intermediate liquid is attributed to the melting of the low-melting-point crystalline phases, rather than induced by the elastic stress energy along the interface proposed in literatures. This could also be a general feature in multicomponent alloys if a low-melting-point crystalline phase is formed together with other crystalline phases during heating. Furthermore, all results obtained here are scientifically sound and provide new insight into the phase transition theory, which will trigger more studies on the microscopic mechanism of solid-solid phase transition via an intermediate liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Sn microalloying on the nucleation of L1 2 Al 3 Zr precipitates in a dilute aluminum-zirconium alloy

While L1 2 -Al 3 Zr nanoprecipitates provide a balance between strengthening and good electrical conductivity, the precipitation of L1 2 -Al 3 Zr in aluminum requires aggressive heat treatments. An improved age-hardening response was observed during isochronal aging of an Al-0.24Zr (wt%) alloy when microalloyed with Sn. A new mechanism termed Low melting point Element-Assisted Nucleation (LEAN) is proposed to explain the lower temperature nucleation of L1 2 -Al 3 Zr precipitates observed in this alloy based on the addition of a low melting point element, such as Sn. Characterization verified the first-principles density functional theory prediction that Zr and Sn atoms cluster during homogenization owing to the favorable binding energy of Zr-Sn-vacancy triplets. Direct microstructural observations revealed these clusters form Sn nanoprecipitates that assist the nucleation of L1 2 -Al 3 Zr at 200°C, where L1 2 -Al 3 Zr precipitation is not expected due to the low diffusivity of Zr atoms in Al. At higher temperatures (≳350°C), the acceleration of L1 2 -Al 3 Zr precipitation is driven by faster Zr diffusion in Al with Sn microalloying and the nuclei formed via the LEAN mechanism. In conclusion, this combination of mechanisms explains the improvement in age hardening through L1 2 -Al 3 Zr precipitation with Sn microalloying.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts

Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.

36 MATERIALS SCIENCE↗