Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical studies”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Enhanced heat transfer through filler-polymer interface by surface-coupling agent in heat-dissipation material: A non-equilibrium molecular dynamics study

Developing a composite material of polymers and micrometer-sized fillers with higher heat conductance is crucial to realize modular packaging of electronic components at higher densities. Enhancement mechanisms of the heat conductance of the polymer-filler interfaces by adding the surface-coupling agent in such a polymer composite material are investigated through the non-equilibrium molecular dynamics (MD) simulation. A simulation system is composed of α-alumina as the filler, bisphenol-A epoxy molecules as the polymers, and model molecules for the surface-coupling agent. The inter-atomic potential between the α-alumina and surface-coupling molecule, which is essential in the present MD simulation, is constructed to reproduce the calculated energies with the electronic density-functional theory. Through the non-equilibrium MD simulation runs, we find that the thermal resistance at the interface decreases significantly by increasing either number or lengths of the surface-coupling molecules and that the effective thermal conductivity of the system approaches to the theoretical value corresponding to zero thermal-resistance at the interface. Detailed analyses about the atomic configurations and local temperatures around the interface are performed to identify heat-transfer routes through the interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Developing of Quaternary Pumped Storage Hydropower for Dynamic Studies

Quaternary pumped storage hydropower (Q-PSH) technology, as one of the new advanced-PSH technology, has been developed by taking advantage of Conventional-PSH (C-PSH) and Adjustable Speed-PSH (AS-PSH). By combining adjustable-speed pump unit and conventional hydropower turbine unit in the quaternary configuration, Q-PSH has the more competitive capability of providing fast power support in the future high renewable penetrated power system. Acting as energy storage (ES), Q-PSH provides promising power supply to deal with the uncertainty and variability from renewable energy generation. This paper focuses on the dynamic modeling of Q-PSH technology employing full-converter machine and the impact of Q-PSH on the frequency response in a system. The detailed model of Q-PSH is developed and implemented in the IEEE 14-bus system based on GE Positive Sequence Load Flow (PSLF) platform, which captures the dynamic of multiple operation modes, especially hydraulic short-circuit (HSC) operation mode. Several cases are set up to reveal the advantages of Q-PSH technology when power electronic based renewable energy generation is deployed in the system. Sensitivity studies of the controller in pump governor show the impact of parameters in pump response performance. The comparison case illustrates the impact of frequency response provided by the Q-PSH in the system.

dynamic modeling↗
DOE OSTI2016

Rotation mechanism of methylammonium molecules in organometal halide perovskite in cubic phase: An ab initio molecular dynamics study

Rotation of methylammonium (CH 3 NH 3 or MA) molecules is believed to govern the excellent transport properties of photocarriers in the MA lead iodide (MAPbI 3 ) perovskite. Of particular interest is its cubic phase, which exists in industrially important films at room temperature. In order to investigate the rotational behaviors of the MA molecules, we have performed ab initio molecular dynamics simulations of cubic-MAPbI 3 at room temperature. There are two types of rotational motions of MA molecules in a crystalline PbI 3 cage: reorientation of a whole molecule and intramolecular rotation around the C–N bond within MA molecules. Using a cubic symmetry-assisted analysis (CSAA), we found that the prominent orientation of the C–N bond is the crystalline ❬110❭directions, rather than the ❬100❭ and ❬111❭ directions. Rapid rotation around the C–N bond is also observed, which easily occurs when the rotational axis is parallel to the ❬110❭ directions according to the CSAA. To explain the atomistic mechanisms underlying these CSAA results, we have focused on the relation between H–I hydrogen bonds and the orientation of an MA molecule. Here, the hydrogen bonds were defined by population analysis, and it has been found that, while H atoms in the CH 3 group (HC) hardly interacts with I atoms, those in the NH3 group (HN) form at least one hydrogen bond with I atoms and their interatomic distances are in a wide range, 2.2–3.7 Å. Based on these findings, we have given a possible explanation to why the ❬110❭ directions are preferred. Namely, the atomic arrangement and interatomic distance between MA and surrounding I atoms are most suitable for the formation of hydrogen bonds. In addition to films, these results are potentially applicable to the rotational behaviors in bulk MAPbI 3 as well, considering that the atomistic structure and time constants regarding the rotation of MA molecules statistically agree with bulk experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Transient analysis techniques in performing impact and crash dynamic studies

Because of the emphasis being placed on crashworthiness as a design requirement, increasing demands are being made by various organizations to analyze a wide range of complex structures that must perform safely when subjected to severe impact loads, such as those generated in a crash event. The ultimate goal of crashworthiness design and analysis is to produce vehicles with the ability to reduce the dynamic forces experienced by the occupants to specified levels, while maintaining a survivable envelope around them during a specified crash event. DYCAST is a nonlinear structural dynamic finite element computer code that started from the plans systems of a finite element program for static nonlinear structural analysis. The essential features of DYCAST are outlined.

Pifko, A. B.↗
DOE OSTI2023

Dynamic study of direct CO 2 capture from indoor air using poly(ethylenimine)-impregnated fiber sorbents

Supported amine adsorbents are promising materials for direct air capture (DAC) of CO 2 due to their high CO 2 capacity and relatively low energy requirement for regeneration. For a DAC process, it is essential to properly define operating parameters to achieve high sorbent productivity (amount of CO 2 captured per unit quantity of sorbent material over unit time). It is furthermore essential to understand the kinetic behavior of the process under the influence of various operating conditions such as the inlet air velocity, sorbent composition, and humidity to select an effective range of operating conditions to maximize sorbent productivity. Here, the dynamic behavior of a DAC process is probed using a fixed fiber sorbent contactor containing poly(ethylenimine) (PEI)-impregnated composite silica/cellulose acetate (CA) fibers. Throughout this study, experiments are conducted using both simulated air (398 ppm CO 2 balanced by N 2 ) and real indoor air (~400–500 ppm CO 2 ). The experimental behavior of the fibers using simulated air and indoor air is compared, and the influence of the inlet air velocity on the breakthrough behavior is assessed. By changing operating conditions, the impact on the fiber sorbent productivity (mmol CO 2 g fiber –1 h –1 ) is quantified to identify conditions that could favor high rates of CO 2 removal. The kinetics of steam-assisted CO 2 desorption are studied, identifying achievable desorption times. Productivities of 1.2 mmol CO 2 g fiber –1 h –1 are obtained using an inlet air velocity of 1.1 m s –1 . Performance trends show that further increasing the inlet air velocity will likely lead to even higher productivities.

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

Experimentally studied dynamic dose interplay does not meaningfully affect target dose in VMAT SBRT lung treatments

Purpose: The effects of respiratory motion on the tumor dose can be divided into the gradient and interplay effects. While the interplay effect is likely to average out over a large number of fractions, it may play a role in hypofractionated [stereotactic body radiation therapy (SBRT)] treatments. This subject has been extensively studied for intensity modulated radiation therapy but less so for volumetric modulated arc therapy (VMAT), particularly in application to hypofractionated regimens. Also, no experimental study has provided full four-dimensional (4D) dose reconstruction in this scenario. The authors demonstrate how a recently described motion perturbation method, with full 4D dose reconstruction, is applied to describe the gradient and interplay effects during VMAT lung SBRT treatments.Methods: VMAT dose delivered to a moving target in a patient can be reconstructed by applying perturbations to the treatment planning system-calculated static 3D dose. Ten SBRT patients treated with 6 MV VMAT beams in five fractions were selected. The target motion (motion kernel) was approximated by 3D rigid body translation, with the tumor centroids defined on the ten phases of the 4DCT. The motion was assumed to be periodic, with the period T being an average from the empirical 4DCT respiratory trace. The real observed tumor motion (total displacement ≤8 mm) was evaluated first. Then, the motion range was artificially increased to 2 or 3 cm. Finally, T was increased to 60 s. While not realistic, making T comparable to the delivery time elucidates if the interplay effect can be observed. For a single fraction, the authors quantified the interplay effect as the maximum difference in the target dosimetric indices, most importantly the near-minimum dose (D{sub 99%}), between all possible starting phases. For the three- and five-fractions, statistical simulations were performed when substantial interplay was found.Results: For the motion amplitudes and periods obtained from the 4DCT, the interplay effect is negligible (<0.2%). It is also small (0.9% average, 2.2% maximum) when the target excursion increased to 2–3 cm. Only with large motion and increased period (60 s) was a significant interplay effect observed, with D{sub 99%} ranging from 16% low to 17% high. The interplay effect was statistically significantly lower for the three- and five-fraction statistical simulations. Overall, the gradient effect dominates the clinical situation.Conclusions: A novel method was used to reconstruct the volumetric dose to a moving tumor during lung SBRT VMAT deliveries. With the studied planning and treatment technique for realistic motion periods, regardless of the amplitude, the interplay has nearly no impact on the near-minimum dose. The interplay effect was observed, for study purposes only, with the period comparable to the VMAT delivery time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Temperature dependent structural behavior of Mg{sub 70}Zn{sub 30}: A molecular dynamics study

Mg{sub 70}Zn{sub 30} is simplest of metallic glasses having good mechanical properties, biocompatibility with human body and is biodegradable, so it is widely studied material in its class. We present in this paper first principle calculations showing the change in the structure of the system with temperature. The computations are performed using PWSCF tool of Quantum ESPRESSO package, which uses density functional theory for calculations of electronic ground state properties and dynamical structures. With increasing temperature the number of peaks in RDFs reduce in number leaving only a single peak,the system thus makes a transition from non-crystalline solid to liquid state within the range of 300–1000K. First coordination number was found to decrease with rising temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mammalian complex III heme dynamics studied with pump-probe spectroscopy and red light illuminations

The electronic or molecular mechanisms which initiate photobiomodulation (PBM) in cells are not yet fully understood. Porcine complex III (C-III) of the electron transport chain was characterized with transient absorption spectroscopy (TAS). Then we applied our recently developed continuous wave laser coupled TAS procedure (CW-TAS) to investigate the effect of red light irradiances on the heme dynamics of C-III in its c 1 reduced state. The time constants associated with the axial ligand photodissociation were found to be 3.3 ± 0.3 ps for the oxidized state and 4.9 ± 0.4 ps for the c 1 reduced state. The analysis of the CW-TAS procedure yielded no significant changes in the C-III heme dynamics. We rule out the possibility of 635 nm CW light at 4.7 mW/cm 2 inducing a PBM effect on the heme dynamic of C-III, specifically with the photodissociation of its axial ligand.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

ReaxFF-based molecular dynamics study of bio-derived polycyclic alkanes as potential alternative jet fuels

This work investigates the initial stages of the pyrolysis of HtH-1 (C 18 H 32 ; 2,2,7,7,8a,8b-hexamethyl-dodecahydrobiphenylene) and HtH-2 (C 18 H 34 ; 1,1',3,3,3',3'-hexamethyl-1,1'-bi(cyclohexane)), which are bio-derived polycyclic alkanes and potential jet fuels, using ReaxFF force field based molecular dynamics (MD) simulations. Global Arrhenius parameters, such as activation energies and pre-exponential factors, are calculated and used to analyze the overall decomposition kinetics of the fuels. HtH-1 decomposes faster than HtH-2 at the same temperature and density conditions, and they have a faster decomposition rate compared to some existing jet-fuels, such as JP-10. A systematic reaction analysis framework developed in this work is applied to determine a temperature-dependent decomposition mechanism. At lower temperature, the central C-C bond connecting the two cyclohexane rings is dominantly broken in both HtH-1 and HtH-2. However, C-CH 3 bond breaking becomes dominant with increasing temperature due to the large increase in entropy during this reaction. Major products from HtH-1 are C 5 H 8 and C 4 H 8 , and those from HtH-2 are C 4 H 8 and C 2 H 4 . The major products predict that HtH-1 has a higher sooting tendency than HtH-2, which is consistent with measurements. The impact of HtH-2 on the pyrolysis of HtH-1 is also investigated in their binary mixtures. HtH-1 and HtH-2 decompose by unimolecular reactions, and they rarely interact with each other during the pyrolysis of the mixtures. Furthermore, this work demonstrates that ReaxFF can be used to investigate pyrolysis and combustion chemistry of existing or future fuels and to contribute to the development of their chemical kinetic models without any a priori input and chemical intuition.

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

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗
DOE OSTI2021

High-salinity brine desalination with amine-based temperature swing solvent extraction: A molecular dynamics study

Solvent extraction desalination (SED) is a promising technique for efficiently treating high salinity brines. However, only a few studies in the literature have investigated the influence of the molecular characteristics of the organic solvents on the emergent performance of the separation process. Thus, in this work we investigate the influence of the molecular structure of four different secondary amines on the desalination process of a concentrated aqueous solution of NaCl. Molecular dynamics simulations and free energy calculations are used to obtain the thermodynamic and structural properties of the pure solvent bulk phase, as well as a homogeneous solvent-water bulk mixture. Furthermore, solvent-water and solvent-brine interfacial systems are used to explicitly analyze the water partitioning behavior in the relevant temperature range for SED. The water solubilization in these secondary amines is found to be an exothermic process. Additionally, although the hydrogen-bonding interaction between the water and solvent molecules are beneficial, this interaction is not necessarily the dominant factor for ranking the water solubility performance of the different solvents. Finally, the excess volume of mixing of the bulk solvent-water mixture, as well as the electrostatic surface features of the solvents are also found to provide important clues for screening solvent behavior.

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

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A full-dimensional quantum dynamical study of H{sub 2}+H{sub 2} collisions: Coupled-states versus close-coupling formulation

Collision-induced energy transfer involving H{sub 2} molecules plays an important role in many areas of physics. Kinetic models often require a complete set of state-to-state rate coefficients for H{sub 2}+H{sub 2} collisions in order to interpret results from spectroscopic observations or to make quantitative predictions. Recent progress in full-dimensional quantum dynamics using the numerically exact close-coupling (CC) formulation has provided good agreement with existing experimental data for low-lying states of H{sub 2} and increased the number of state-to-state cross sections that may be reliably determined over a broad range of energies. Nevertheless, there exist many possible initial states (e.g., states with high rotational excitation) that still remain elusive from a computational standpoint even at relatively low collision energies. In these cases, the coupled-states (CS) approximation offers an alternative full-dimensional formulation. We assess the accuracy of the CS approximation for H{sub 2}+H{sub 2} collisions by comparison with benchmark results obtained using the CC formulation. The results are used to provide insight into the orientation effects of the various internal energy transfer mechanisms. A statistical CS approximation is also investigated and cross sections are reported for transitions which would otherwise be impractical to compute.

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

A molecular dynamics study of laser-excited gold

The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melting process, leading to large discrepancies between the predicted threshold energy density for complete melting, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. [Science 360, 1451–1455 (2018)]. We recreate the experimental electron diffraction data, employing both a constant and temperature-dependent electron–ion equilibration rate. In all cases, we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Na ion dynamics in P2-Na x [Ni 1/3 Ti 2/3 ]O 2 : a combination of quasi-elastic neutron scattering and first-principles molecular dynamics study

In this work, the P2-type layered material Na 2/3 [Ni 1/3 Ti 2/3 ]O 2 was studied as a promising bi-functional electrode material for sodium-ion batteries. To assess the electrochemical performance of this material, we investigated the diffusion mechanism as well as ionic and electronic conductivity with a combination of experimental and computational techniques. The quasi-elastic neutron scattering (QENS) experiments and first-principles molecular dynamics (FPMD) simulations were performed to identify the diffusion mechanism. The QENS data showed that Na ion diffusion can be well described by the Singwi–Sjölander jump diffusion model, where the obtained mean jump length matched the distances between the neighboring edge-share and face-share Na sites. FPMD predicted diffusivity values similar to those from QENS. The computed composition dependence of ionic and electronic conductivity of Na x [Ni 1/3 Ti 2/3 ]O 2 suggested that electronic conductivity changes significantly when x deviates from 2/3 as the redox couple of Ni and Ti is activated, while the change of ionic conductivity with x is relatively small.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A molecular dynamics study of collisional heat transfer to nanoclusters in the gas phase

Nanoclusters in the gas phase grow by condensation and coagulation, which heat up nanoclusters, creating thermal non-equilibrium with the surrounding gas. The extent of non-equilibrium for nanoclusters is dependent upon the rate of heating relative to the rate of cooling brought about by nanocluster-gas molecule collisions, hence the latter process plays an important role in nanocluster energy evolution during growth. We apply molecular dynamics simulations to investigate heat transfer between metal nanoclusters and gas molecules relevant to high temperature aerosol systems. In analysis, we first define a thermal reemission coefficient (TRC), which is a correction factor for the portion of energy possessed by reemitted gas molecules after collision. Here, the TRC differs from the thermal accommodation coefficient (TAC), which is a correction factor for the maximum possible energy transfer rate from nanoclusters to gas molecules based on the diffusive collision model. We find the TRC is insensitive to nanocluster size but is affected by the nanocluster to gas temperature ratio and the nanocluster to gas atomic mass ratio. Comparison to the traditionally-defined TAC shows that the TAC is insensitive to temperature, a phenomenon that has been widely indicated in literature but not quantitatively explained. We show how the temperature insensitivity of the TAC arises because of the manner in which TRCs vary with both nanocluster and gas temperature, i.e. the manner in which the reemitted gas molecule energy from a surface varies with the temperature of that surface relative to that of the gas. In addition, while calculations ultimately reveal temperature insensitive TACs and justify their continued use in modeling gas-phase heat transfer, we argue that the TRC is more appropriate to calculate a priori in analyzing simulations; it is physically more realistic to first introduce a correction factor solely for the energy of reemitted gas molecules after a collision (as TRCs) instead of a correction factor applying to both the reemission energy and the initial energy (as TACs), as the latter is exactly calculable and independent of nanocluster properties. In line with prior work, for monoatomic gases we find the TAC decreases as the nanocluster to gas atomic mass increases, approaching 0.08 for gold nanoclusters in helium.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Experimental nonlinear dynamical studies in cesium magneto-optical trap using time-series analysis

A magneto-optical trap of neutral atoms is essentially a dissipative quantum system. The fast thermal atoms continuously dissipate their energy to the environment via spontaneous emissions during the cooling. The atoms are, therefore, strongly coupled with the vacuum reservoir and the laser field. The vacuum fluctuations as well as the field fluctuations are imparted to the atoms as random photon recoils. Consequently, the external and internal dynamics of atoms becomes stochastic. In this paper, we have investigated the stochastic dynamics of the atoms in a magneto-optical trap during the loading process. The time series analysis of the fluorescence signal shows that the dynamics of the atoms evolves, like all dissipative systems, from deterministic to the chaotic regime. The subsequent disappearance and revival of chaos was attributed to chaos synchronization between spatially different atoms in the magneto-optical trap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An ab initio molecular dynamics study of varied compositions of the $\text{LiF-NaF-KF}$ molten salt

With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.

36 MATERIALS SCIENCE↗