Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical force fields”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Magnetothermal Convection in Nonconducting Diamagnetic and Paramagnetic Fluids

Nonuniform magnetic fields exert a magnetic body force on electrically nonconducting classical fluids. These include paramagnetic fluids such as gaseous and liquid oxygen and diamagnetic fluids such as helium. Recent experiments show that this force can overwhelm the force of gravity even at the surface of the earth; it can levitate liquids and gases, quench candle flames, block gas flows, and suppress heat transport. Thermal gradients render the magnetic force nonuniform through the temperature-dependent magnetic susceptibility. These thermal gradients can therefore drive magnetic convection analogous to buoyancy-driven convection. This magnetothermal convection can overwhelm convection driven by gravitational buoyancy in terrestrial experiments. The objectives of the proposed ground-based theoretical study are (a) to supply the magnetothermohydrodynamic theory necessary to understand these recent experiments and (b) to explore the consequences of nonuniform magnetic fields in microgravity. Even the linear theory for the onset of magnetothermal convection is lacking in the literature. We intend to supply the linear and nonlinear theory based on the thermohydrodynamic equations supplemented by the magnetic body force. We intend to investigate the effect of magnetic fields on gas blockage and heat transport in microgravity. Since magnetic fields provide a means of creating arbitrary, controllable body force distributions, we intend to investigate the possibility of using magnetic fields to position and control fluids in microgravity. We also intend to investigate the possibility of creating stationary terrestrial microgravity environments by using the magnetic force to effectively cancel gravity. These investigations may aid in the design of space-based heat-transfer, combustion, and human-life-support equipment.

Edwards, Boyd F.↗
NASA NTRS1978

Expansion and broadening of coronal loop transients - A theoretical explanation

Consequences are examined of the assumption that an observed coronal loop transient is a twisted rope of magnetic-field lines expanding and broadening in the background coronal plasma and magnetic field. It is shown that the expansion can be accounted for by the azimuthal component of the field; the observed broadening of the loop as it moves outward can be accounted for by the longitudinal component of the field. In order to have a net outward force and at the same time avoid a classical pinch (sausage) instability, the two components of the field must satisfy a certain inequality. It is predicted that, as the loop rises, the width (h) of its top portion should vary proportionally with distance (R) from the sun's center. This is in good agreement with measurements that show h is proportional to the 0.8 power of R. The prediction that the radius of curvature of the top portion of the loop should be proportional to R differs from the measured variation. The difference could be accounted for by a drag due to the background coronal field that flattens the loop's top.

Mouschovias, T. CH.↗
DOE OSTI2025

Anomalous resistivity and electron heating by lower hybrid drift waves inside reconnecting current sheets

Inside an electron diffusion region of laboratory reconnection experiments, the quasi-electrostatic lower hybrid drift wave (ES-LHDW) is observed when a significant guide field component is present. Through direct measurement of the anomalous drag term and quasilinear analysis, it is shown that ES-LHDW can account for approximately 20% of the mean reconnection electric field in a case with moderate guide field. This value exceeds the contribution from classical resistivity, which is around 10%. The effects of the Lorentz force term, often neglected for electrostatic waves, are crucial for the observed correlation between electric field and density fluctuations. Anomalous electron heating by the perturbed current and resistivity (2.6 MW/m 3 ) also surpasses the classical Ohmic heating, which is about 2.0 MW/m 3 . For the case with a high guide field, significantly higher local electron temperatures were observed during periods of strong ES-LHDW activity. A statistical analysis further supports electron heating by LHDW, showing a larger increase in electron temperature with a high guide field. Finally, data from the Magnetospheric Multiscale mission provide evidence of Landau damping of ES-LHDW, suggesting that ES-LHDW may contribute to the generation of nonthermal electrons along the direction parallel to the magnetic field.

Magnetic reconnection↗
NASA NTRS1998

Computational Investigation of Helical Traveling Wave Tube Transverse RF Field Forces

In a previous study using a fully three-dimensional (3D) helical slow-wave circuit cold- test model it was found, contrary to classical helical circuit analyses, that transverse FF electric fields have significant amplitudes compared with the longitudinal component. The RF fields obtained using this helical cold-test model have been scaled to correspond to those of an actual TWT. At the output of the tube, RF field forces reach 61%, 26% and 132% for radial, azimuthal and longitudinal components, respectively, compared to radial space charge forces indicating the importance of considering them in the design of electron beam focusing.

Kory, Carol L.↗
DOE OSTI2023

Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence

With the advent of faster computer processors and especially graphics processing units (GPUs) over the last few decades, the use of data-intensive machine learning (ML) and artificial intelligence (AI) has increased greatly, and the study of crystal nucleation has been one of the beneficiaries. In this study, we outline how ML and AI have been applied to address four outstanding difficulties of crystal nucleation: how to discover better reaction coordinates (RCs) for describing accurately non-classical nucleation situations; the development of more accurate force fields for describing the nucleation of multiple polymorphs or phases for a single system; more robust identification methods for determining crystal phases and structures; and as a method to yield improved course-grained models for studying nucleation.

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

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

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

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2022

Hydration of divalent lanthanides, Sm 2+ and Eu 2+ : A molecular dynamics study with polarizable AMOEBA force field

The chemistry of divalent lanthanides, Ln 2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this study, we address the hydration of two of the classical Ln 2+ cations, Sm 2+ and Eu 2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm 2+ and Eu 2+ , and 8.8 for both Sm 3+ and Eu 3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln 2+ compared with Ln 3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln 2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln 2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu 2+ and Sm 2+ in aqueous chloride media.

AMOEBA force field parameterization↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

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

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1990

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

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

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
DOE OSTI2023

Electric-Field-Induced Assembly of an Ionic Liquid–Water Interphase Enables Efficient Heavy Metal Electrosorption

Controlling ion desolvation, transport, and charge transfer at the electrode-electrolyte interface (EEI) is critical to enabling the rational design of efficient and selective separation of targeted heavy metals and decontamination of industrial waste water. One challenge is to sufficiently resolve and interrogate the intermediate transformation steps between solvated metal cations and bare cations for electroreduction at the EEI and establish pathways to modulate these steps to achieve efficient energy transfer for targeted reactive separations. Herein, we obtained a predictive understanding of the effect of adding a hydrophobic ionic liquid, 1-ethyl-3-methylimidazolium chloride (EMIMCl), to an aqueous electrolyte on modulating the desolvation and electrosorption of Pb 2+ cations using a combination of experimental and theoretical techniques. Here, we reveal formation of a compact interphase layer consisting of EMIMCl-Pb complex clusters under an applied electric field using operando electrochemical Raman spectroscopy, atomic force microscopy, and electrochemical impedance spectroscopy measurements combined with classical molecular dynamics simulations. The application of a larger negative potential is shown to result in formation of a well oriented layer with the positive imidazolium ring of the EMIMCl lying on the electrode and the long hydrophobic alkyl chain extending into the bulk electrolyte. This oriented layer is demonstrated to facilitate desolvation of incoming solvated Pb 2+ cations and decrease the charge transfer resistance for Pb electrodeposition, which has important implications for selective removal of Pb from contaminated mixtures. Overall, our studies open up new opportunities to modulate ion desolvation using hydrophobic ionic liquids in aqueous electrolytes for efficient heavy metal separation.

36 MATERIALS SCIENCE↗