Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DYNAMICS”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Application of Dynamical System Scaling to Bubble Dynamics

Dynamical System Scaling (DSS) provides a useful method for analyzing, categorizing and scaling time-dependent processes. A key feature of DSS is the temporal displacement rate, D, which relates the natural process time to the reference clock time. Because of its property of being invariant under a two-parameter affine transformation, it provides an underlying basis for process scaling. Application of DSS to bubble dynamics serves to elucidate the physical significance of the temporal displacement rate and its role in scaling a variety of bubble dynamics processes. This paper shows that the temporal displacement rate consists of the sum of dimensionless groups that govern the bubble dynamics. Preserving the temporal displacement rate for a prototype and a scaled model results in similitude of the time-dependent normalized bubble size, growth rates and interfacial accelerations for different fluid conditions. For vapor bubble growth in a superheated liquid, it is shown that inertial controlled bubble growth occurs when D=0 and thermally controlled bubble growth occurs when D=1 .

dynamical system scaling, SMR, bubble dynamics, DS↗
DOE OSTI2017

Atomic-level heterogeneity and defect dynamics in concentrated solid-solution alloys

Performance enhancement of structural materials in extreme radiation environments has been activelyinvestigated for many decades. Traditional alloys, such as steel, brass and aluminum alloys, normallycontain one or two principal element(s) with a low concentration of other elements. While these existin either a mixture of metallic phases (multiple phases) or in a solid solution (single phase), limited orlocalized chemical disorder is a common characteristic of the main matrix. Fundamentally different fromtraditional alloys, recently developed single-phase concentrated solid-solutio n alloys (CSAs) contain mul-tiple elemental species in equiatomic or high concentrations with different elements randomly arrangedon a crystalline lattice. Due to the lack of ordered elemental arrangement in these CSAs, they exhibit sig-nificant chemical disorder and unique site-to-site lattice distortion. While it is well recognized intraditional alloys that minor additions lead to enhanced radiation resistance, it remains unclear inCSAs how atomic-level heterogeneity affects defect formation, damage accumulation, and microstruc-tural evolution. These knowledge gaps have acted as roadblocks to the development of future-generation energy technology. CSAs with a simple crystal structure, but complex chemical disorder,are unique systems that allow us, through replacing principal alloying elements and modifying concen-trations, to study how compositional complexity influences defect dynamics, and to bridge the knowl-edge gaps through understanding intricate electronic- and atomic-level interactions, mass and energytransfer processes, and radiation resistance performance. Recent advances in defect dynamics and irradi-ation performance of CSAs are reviewed, intrinsic chemical effects on radiation performance arediscussed, and direction for future studies is suggested.

Zhang, Yanwen↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Dynamic Process Intensification of Binary Distillation via Periodic Operation

This paper applies the concept of dynamic intensification (defined as changes to the dynamics, operation strategy, and/or control of a process that lead to a substantially more efficient processing path) to binary distillation columns. The resulting strategy consists of manufacturing a target product as a blend of two auxiliary products, both having lower energy demands than a reference value, which corresponds to producing the target product(s) in a column operating at steady state. Here, a discussion of the appropriate control structures and switching strategies between the two auxiliary products is provided. An extensive case study concerning the separation of a methanol–1-propanol mixture was carried out, demonstrating that energy savings in the order of 1.4% are possible with no disruption in product quality or production rate.

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

Self-Consistent Theory for Structural Relaxation, Dynamic Bond Exchange Times, and the Glass Transition in Polymeric Vitrimers

Here, we formulate a statistical mechanical theory for how dynamic bond exchange influences the activated hopping-driven relaxation of Kuhn segments in dynamically cross-linked networks or vitrimers over a wide range of temperatures and cross-link densities. The key new methodological aspect is to address in a self-consistent manner the dynamic consequences of bond exchange on the Kuhn segmental alpha relaxation, and vice versa. The predicted temperature dependence of the segmental alpha time of vitrimers at high temperatures remains the same as that of permanent networks, but at lower temperatures, a significant acceleration of relaxation occurs due to bond exchanges. From a mechanistic perspective, the vitrimer local cage barrier is very weakly affected by bond exchange, while the collective elastic barrier contribution decreases significantly in the deeply supercooled regime. The vitrimer glass transition temperature is predicted to grow linearly with the square root of the cross-link density, as previously found for permanent networks. Material-specific chemical effects such as cross-linker size relative to that of the normal Kuhn segment or special attraction of a cross-linker with polymers are crudely considered based on model calculations. The theory is quantitatively applied to recent experiments on dry ethylene vitrimers. Good agreements are found including that the bond exchange time follows an Arrhenius law at high enough temperatures but upward non-Arrhenius deviations emerge in the deeply supercooled regime, and a collapsed master curve of the segment alpha time exists over a wide range of cross-link densities and temperatures. Possible extensions to treat dynamic heterogeneity effects and penetrant transport in vitrimers are briefly discussed.

chemical calculations↗
DOE OSTI2022

MultiSector Dynamics: Advancing the Science of Complex Adaptive Human-Earth Systems

The field of MultiSector Dynamics (MSD) explores the dynamics and co-evolutionary pathways of human and Earth systems with a focus on critical goods, services, and amenities delivered to people through interdependent sectors. This commentary lays out core definitions and concepts, identifies MSD science questions in the context of the current state of knowledge, and describes ongoing activities to expand capacities for open science, leverage revolutions in data and computing, and grow and diversify the MSD workforce. Central to our vision is the ambition of advancing the next generation of complex adaptive human-Earth systems science to better address interconnected risks, increase resilience, and improve sustainability. This will require convergent research and the integration of ideas and methods from multiple disciplines. Understanding the tradeoffs, synergies, and complexities that exist in coupled human-Earth systems is particularly important in the context of energy transitions and increased future shocks.

Reed, Patrick↗
DOE OSTI2020

Biomolecular complex viewed by dynamic nuclear polarization solid-state NMR spectroscopy

Solid-state nuclear magnetic resonance (ssNMR) is an indispensable tool for elucidating the structure and dynamics of insoluble and non-crystalline biomolecules. The recent advances in the sensitivity-enhancing technique magic-angle spinning dynamic nuclear polarization (MAS-DNP) have substantially expanded the territory of ssNMR investigations and enabled the detection of polymer interfaces in a cellular environment. This article highlights the emerging MAS-DNP approaches and their applications to the analysis of biomolecular composites and intact cells to determine the folding pathway and ligand binding of proteins, the structural polymorphism of low-populated biopolymers, as well as the physical interactions between carbohydrates, proteins, and lignin. Furthermore, these structural features provide an atomic-level understanding of many cellular processes, promoting the development of better biomaterials and inhibitors. It is anticipated that the capabilities of MAS-DNP in biomolecular and biomaterial research will be further enlarged by the rapid development of instrumentation and methodology.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Slow dynamics of a protein backbone in molecular dynamics simulation revealed by time-structure based independent component analysis

We recently proposed the method of time-structure based independent component analysis (tICA) to examine the slow dynamics involved in conformational fluctuations of a protein as estimated by molecular dynamics (MD) simulation [Y. Naritomi and S. Fuchigami, J. Chem. Phys. 134, 065101 (2011)]. Our previous study focused on domain motions of the protein and examined its dynamics by using rigid-body domain analysis and tICA. However, the protein changes its conformation not only through domain motions but also by various types of motions involving its backbone and side chains. Some of these motions might occur on a slow time scale: we hypothesize that if so, we could effectively detect and characterize them using tICA. In the present study, we investigated slow dynamics of the protein backbone using MD simulation and tICA. The selected target protein was lysine-, arginine-, ornithine-binding protein (LAO), which comprises two domains and undergoes large domain motions. MD simulation of LAO in explicit water was performed for 1 μs, and the obtained trajectory of C{sub α} atoms in the backbone was analyzed by tICA. This analysis successfully provided us with slow modes for LAO that represented either domain motions or local movements of the backbone. Further analysis elucidated the atomic details of the suggested local motions and confirmed that these motions truly occurred on the expected slow time scale.

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

Structural, dynamic, and vibrational properties during heat transfer in Si/Ge superlattices: A Car-Parrinello molecular dynamics study

The structural, dynamic, and vibrational properties during heat transfer process in Si/Ge superlattices are studied by analyzing the trajectories generated by the ab initio Car-Parrinello molecular dynamics simulation. The radial distribution functions and mean square displacements are calculated and further discussions are made to explain and probe the structural changes relating to the heat transfer phenomenon. Furthermore, the vibrational density of states of the two layers (Si/Ge) are computed and plotted to analyze the contributions of phonons with different frequencies to the heat conduction. Coherent heat conduction of the low frequency phonons is found and their contributions to facilitate heat transfer are confirmed. The Car-Parrinello molecular dynamics simulation outputs in the work show reasonable thermophysical results of the thermal energy transport process and shed light on the potential applications of treating the heat transfer in the superlattices of semiconductor materials from a quantum mechanical molecular dynamics simulation perspective.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Ab initio molecular dynamics simulations reveal localization and time evolution dynamics of an excess electron in heterogeneous CO{sub 2}–H{sub 2}O systems

In view of the important implications of excess electrons (EEs) interacting with CO{sub 2}–H{sub 2}O clusters in many fields, using ab initio molecular dynamics simulation technique, we reveal the structures and dynamics of an EE associated with its localization and subsequent time evolution in heterogeneous CO{sub 2}–H{sub 2}O mixed media. Our results indicate that although hydration can increase the electron-binding ability of a CO{sub 2} molecule, it only plays an assisting role. Instead, it is the bending vibrations that play the major role in localizing the EE. Due to enhanced attraction of CO{sub 2}, an EE can stably reside in the empty, low-lying π{sup *} orbital of a CO{sub 2} molecule via a localization process arising from its initial binding state. The localization is completed within a few tens of femtoseconds. After EE trapping, the ∠OCO angle of the core CO{sub 2}{sup −} oscillates in the range of 127°∼142°, with an oscillation period of about 48 fs. The corresponding vertical detachment energy of the EE is about 4.0 eV, which indicates extreme stability of such a CO{sub 2}-bound solvated EE in [CO{sub 2}(H{sub 2}O){sub n}]{sup −} systems. Interestingly, hydration occurs not only on the O atoms of the core CO{sub 2}{sup −} through formation of O⋯H–O H–bond(s), but also on the C atom, through formation of a C⋯H–O H–bond. In the latter binding mode, the EE cloud exhibits considerable penetration to the solvent water molecules, and its IR characteristic peak is relatively red-shifted compared with the former. Hydration on the C site can increase the EE distribution at the C atom and thus reduce the C⋯H distance in the C⋯H–O H–bonds, and vice versa. The number of water molecules associated with the CO{sub 2}{sup −} anion in the first hydration shell is about 4∼7. No dimer-core (C{sub 2}O{sub 4}{sup −}) and core-switching were observed in the double CO{sub 2} aqueous media. This work provides molecular dynamics insights into the localization and time evolution dynamics of an EE in heterogeneous CO{sub 2}–H{sub 2}O media.

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

Kinetically constrained ring-polymer molecular dynamics for non-adiabatic chemical reactions

We extend ring-polymer molecular dynamics (RPMD) to allow for the direct simulation of general, electronically non-adiabatic chemical processes. The kinetically constrained (KC) RPMD method uses the imaginary-time path-integral representation in the set of nuclear coordinates and electronic states to provide continuous equations of motion that describe the quantized, electronically non-adiabatic dynamics of the system. KC-RPMD preserves the favorable properties of the usual RPMD formulation in the position representation, including rigorous detailed balance, time-reversal symmetry, and invariance of reaction rate calculations to the choice of dividing surface. However, the new method overcomes significant shortcomings of position-representation RPMD by enabling the description of non-adiabatic transitions between states associated with general, many-electron wavefunctions and by accurately describing deep-tunneling processes across asymmetric barriers. We demonstrate that KC-RPMD yields excellent numerical results for a range of model systems, including a simple avoided-crossing reaction and condensed-phase electron-transfer reactions across multiple regimes for the electronic coupling and thermodynamic driving force.

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

Ring-polymer molecular dynamics: Rate coefficient calculations for energetically symmetric (near thermoneutral) insertion reactions (X + H{sub 2}) → HX + H(X = C({sup 1}D), S({sup 1}D))

Following our previous study of prototypical insertion reactions of energetically asymmetric type with the RPMD (Ring-Polymer Molecular Dynamics) method [Y. Li, Y. Suleimanov, and H. Guo, J. Phys. Chem. Lett. 5, 700 (2014)], we extend it to two other prototypical insertion reactions with much less exothermicity (near thermoneutral), namely, X + H{sub 2} → HX + H where X = C({sup 1}D), S({sup 1}D), in order to assess the accuracy of this method for calculating thermal rate coefficients for this class of reactions. For both chemical reactions, RPMD displays remarkable accuracy and agreement with the previous quantum dynamic results that make it encouraging for the future application of the RPMD to other barrier-less, complex-forming reactions involving polyatomic reactants with any exothermicity.

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

Ion transport and structural dynamics in homologous ammonium and phosphonium-based room temperature ionic liquids

Charge transport and structural dynamics in a homologous pair of ammonium and phosphonium based room temperature ionic liquids (ILs) have been characterized over a wide temperature range using broadband dielectric spectroscopy and quasi-elastic light scattering spectroscopy. We have found that the ionic conductivity of the phosphonium based IL is significantly enhanced relative to the ammonium homolog, and this increase is primarily a result of a lower glass transition temperature and higher ion mobility. Additionally, these ILs exhibit pronounced secondary relaxations which are strongly influenced by the atomic identity of the cation charge center. While the secondary relaxation in the phosphonium IL has the expected Arrhenius temperature dependence characteristic of local beta relaxations, the corresponding relaxation process in the ammonium IL was found to exhibit a mildly non-Arrhenius temperature dependence in the measured temperature range—indicative of molecular cooperativity. These differences in both local and long-range molecular dynamics are a direct reflection of the subtly different inter-ionic interactions and mesoscale structures found in these homologous ILs.

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

Stochastic dynamics of penetrable rods in one dimension: Entangled dynamics and transport properties

The dynamical properties of a system of soft rods governed by stochastic hard collisions (SHCs) have been determined over a varying range of softness using molecular dynamics simulations in one dimension and analytic theory. The SHC model allows for interpenetration of the system’s constituent particles in the simulations, generating overlapping clustering behavior analogous to the spatial structures observed in systems governed by deterministic bounded potentials. Through variation of an assigned softness parameter δ, the limiting ranges of intermolecular softness are bridged, connecting the limiting ensemble behavior from hard to ideal (completely soft). Various dynamical and structural observables are measured from simulation and compared to developed theoretical values. The spatial properties are found to be well predicted by theories developed for the deterministic penetrable-sphere model with a transformation from energetic to probabilistic arguments. While the overlapping spatial structures are complex, the dynamical properties can be adequately approximated through a theory built on impulsive interactions with Enskog corrections. Our theory suggests that as the softness of interaction is varied toward the ideal limit, correlated collision processes are less important to the energy transfer mechanism, and Markovian processes dominate the evolution of the configuration space ensemble. For interaction softness close to hard limit, collision processes are highly correlated and overlapping spatial configurations give rise to entanglement of single-particle trajectories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Microwave soft x-ray microscopy for nanoscale magnetization dynamics in the 5–10 GHz frequency range

We present a scanning transmission x-ray microscopy setup combined with a novel microwave synchronization scheme for studying high frequency magnetization dynamics at synchrotron light sources. The sensitivity necessary to detect small changes in the magnetization on short time scales and nanometer spatial dimensions is achieved by combining the excitation mechanism with single photon counting electronics that is locked to the synchrotron operation frequency. Our instrument is capable of creating direct images of dynamical phenomena in the 5-10 GHz range, with high spatial resolution. When used together with circularly polarized x-rays, the above capabilities can be combined to study magnetic phenomena at microwave frequencies, such as ferromagnetic resonance (FMR) and spin waves. We demonstrate the capabilities of our technique by presenting phase resolved images of a ∼6 GHz nanoscale spin wave generated by a spin torque oscillator, as well as the uniform ferromagnetic precession with ∼0.1° amplitude at ∼9 GHz in a micrometer-sized cobalt strip.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Molecular dynamics simulations of the structure and single-particle dynamics of mixtures of divalent salts and ionic liquids

We report a molecular dynamics study of the structure and single-particle dynamics of mixtures of a protic (ethylammonium nitrate) and an aprotic (1-butyl-3-methylimidazolium hexaflurophosphate [BMIM][PF{sub 6}]) room-temperature ionic liquids doped with magnesium and calcium salts with a common anion at 298.15 K and 1 atm. The solvation of these divalent cations in dense ionic environments is analyzed by means of apparent molar volumes of the mixtures, radial distribution functions, and coordination numbers. For the protic mixtures, the effect of salt concentration on the network of hydrogen bonds is also considered. Moreover, single-particle dynamics of the salt cations is studied by means of their velocity autocorrelation functions and vibrational densities of states, explicitly analyzing the influence of salt concentration, and cation charge and mass on these magnitudes. The effect of the valency of the salt cation on these properties is considered comparing the results with those for the corresponding mixtures with lithium salts. We found that the main structural and dynamic features of the local solvation of divalent cations in ionic liquids are similar to those of monovalent salts, with cations being localized in the polar nanoregions of the bulk mixture coordinated in monodentate and bidentate coordination modes by the [NO{sub 3}]{sup −} and [PF{sub 6}]{sup −} anions. However, stronger electrostatic correlations of these polar nanoregions than in mixtures with salts with monovalent cations are found. The vibrational modes of the ionic liquid (IL) are seen to be scarcely affected by the addition of the salt, and the effect of mass and charge on the vibrational densities of states of the dissolved cations is reported. Cation mass is seen to exert a deeper influence than charge on the low-frequency vibrational spectra, giving a red shift of the vibrational modes and a virtual suppression of the higher energy vibrational modes for the heavier Ca{sup 2+} cations. No qualitative difference with monovalent cations was found in what solvation is concerned, which suggests that no enhanced reduction of the mobility of these cations and their complexes in ILs respective to those of monovalent cations is to be expected.

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

Insight into the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} eutectic mixture from classical molecular dynamics: Thermodynamics, structure, and dynamics

We use molecular dynamics simulations to study the thermodynamics, structure, and dynamics of the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} (62:38 mol. %) eutectic mixture. We present a new classical non-polarizable force field for this molten salt mixture, optimized using experimental and first principles molecular dynamics simulations data as reference. This simple force field allows efficient molecular simulations of phenomena at long time scales. We use this optimized force field to describe the behavior of the eutectic mixture in the 900–1100 K temperature range, at pressures between 0 and 5 GPa. After studying the equation of state in these thermodynamic conditions, we present molecular insight into the structure and dynamics of the melt. In particular, we present an analysis of the temperature and pressure dependence of the eutectic mixture’s self-diffusion coefficients, viscosity, and ionic conductivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗