Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulation”

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 181 records · Page 10

ICME for NASA Aerospace Applications: Batteries for Electric Aviation

NASA’s approach to computational materials modeling is detailed in the NASA Vision 2040 Roadmap for Multiscale Modeling and Simulation of Materials and Systems. This report is in the spirit of national initiatives such as the Material Genome Initiative (MGI), Integrated Computational Materials Engineering (ICME), and others. We utilize a combination of fundamental modeling, computational high-throughput screening, and data science methods, e.g., machine learning, are used to find innovative solutions to NASA or national technology challenges. Applications of interest are wide ranging from advanced alloys to batteries to coatings, among others. In this talk, we present three examples for recent work related to NASA applications. First, doping advanced sulfur battery cathodes with selenium boosts electrical conductivity important for electric aircraft applications. First principles calculations will be discussed that result in compositional design maps for these materials. Second, development of icephobic coatings is important to mitigate safety hazards associated with icing for aircraft. Molecular dynamics simulations are reported for ice-surface interfaces to understand adhesion mechanisms and help screen optimal ice-phobic coatings. Third, shape memory alloys have numerous applications as actuators, superelastic materials, etc. for aerospace. We report machine learning models that predict martensitic transition temperatures across a broad swath of compositional space.

John Lawson↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale repeating unit cell model of a single spherulite containing four disparate length scales was developed to predict the thermoelastic behavior of semicrystalline thermoplastic materials for composite aerospace applications. The continuum level scales were fully coupled and modeled using the generalized method of cells and the high fidelity generalized method of cells micromechanics theories. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents in the multiscale continuum models. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, showing good agreement with the available experimental data from the open literature. Moreover, it is shown that predicted properties are fairly insensitive to the fidelity of the micromechanics model used at the highest continuum scale or the assumed shape of the spherulite.

thermoplastics↗

Interactions of anesthetics with the membrane-water interface

Although the potency of conventional anesthetics correlates with lipophilicity, an affinity to water also is essential. It was recently found that compounds with very low affinities to water do not produce anesthesia regardless of their lipophilicity. This finding implies that clinical anesthesia might arise because of interactions at molecular sites near the interface of neuronal membranes with the aqueous environment and, therefore, might require increased concentrations of anesthetic molecules at membrane interfaces. As an initial test of this hypothesis, we calculated in molecular dynamics simulations the free energy profiles for the transfer of anesthetic 1,1,2-trifluoroethane and nonanesthetic perfluoroethane across water-membrane and water-hexane interfaces. Consistent with the hypothesis, it was found that trifluoroethane, but not perfluoroethane, exhibits a free energy minimum and, therefore, increased concentrations at both interfaces. The transfer of trifluoroethane from water to the nonpolar hexane or interior of the membrane is accompanied by a considerable, solvent-induced shift in the conformational equilibrium around the C-C bond.

NASA Center ARC↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics↗

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗

Interactions of anesthetics with the water-hexane interface. A molecular dynamics study

The free energy profiles characterizing the transfer of nine solutes across the liquid-vapor interfaces of water and hexane and across the water-hexane interface were calculated from molecular dynamics simulations. Among the solutes were n-butane and three of its halogenated derivatives, as well as three halogenated cyclobutanes. The two remaining molecules, dichlorodifluoromethane and 1,2-dichloroperfluoroethane, belong to series of halo-substituted methanes and ethanes, described in previous studies (J. Chem. Phys. 1996, 104, 3760; Chem. Phys. 1996, 204, 337). Each series of molecules contains structurally similar compounds that differ greatly in anesthetic potency. The accuracy of the simulations was tested by comparing the calculated and the experimental free energies of solvation of all nine compounds in water and in hexane. In addition. the calculated and the measured surface excess concentrations of n-butane at the water liquid-vapor interface were compared. In all cases, good agreement with experimental results was found. At the water-hexane interface, the free energy profiles for polar molecules exhibited significant interfacial minima, whereas the profiles for nonpolar molecules did not. The existence of these minima was interpreted in terms of a balance between the free energy contribution arising from solute-solvent interactions and the work to form a cavity that accommodates the solute. These two contributions change monotonically, but oppositely, across the interface. The interfacial solubilities of the solutes, obtained from the free energy profiles, correlate very well with their anesthetic potencies. This is the case even when the Meyer-Overton hypothesis, which predicts a correlation between anesthetic potency and solubility in oil, fails.

NASA Center ARC↗

Filtering Water by Use of Ultrasonically Vibrated Nanotubes

Devices that could be characterized as acoustically driven molecular sieves have been proposed for filtering water to remove all biological contaminants and all molecules larger than water molecules. Originally intended for purifying wastewater for reuse aboard spacecraft, these devices could also be attractive for use on Earth in numerous settings in which there are requirements to obtain potable, medical-grade, or otherwise pure water from contaminated water supplies. These devices could also serve as efficient means of removing some or all water from chemical products . for example, they might be useful as adjuncts or substitutes for stills in the removal of water from alcohols and alcoholic beverages. These devices may be constructed using various materials, such as ceramics, metallics, or polymers, depending on end-use requirements. A representative device of this type (see figure) would include a polymeric disk, about 1 mm in diameter and between 1 and 40 microns thick, within which would be embedded single-wall carbon nanotubes aligned along the thickness axis. The polymeric disk would be part of a unitary polymeric ring assembly. An acoustic transducer in the form of a piezoelectric-film-and-electrode subassembly - typically 9 microns thick and made of poly(vinylidene fluoride) coated with copper 150 nm thick -. would be affixed to the outside of the outer polymeric ring by means of an electrically nonconductive epoxy. The nanotubes would be chosen to have diameters between about 8 and about 13.5 A because water molecules could fit into the nanotubes, but larger molecules could not. Water to be purified would be placed in contact with one face (typically, the upper face) of the filter disk. The surface tension of water is low enough that water molecules should enter and travel along the nanotubes, and computational simulations of molecular dynamics and experimental measurements have shown that the water molecules inside the nanotubes in this size range can be expected to become aligned into helical columns that exhibit properties of both hexagonal ice crystals and liquid water

Gavalas, Lillian Susan↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

molecular dynamics↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

ionic liquids↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI] and [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-Nbutylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-Npropylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of x(sub Li) we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing x(sub Li), the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of x(sub Li) is approximately 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1 - 0.3 mS/cm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

lithium batteries↗

Nonequilibrium shock wave structure, 1. Kinetics of nitric oxide formation and decomposition, 2

Theoretical and numerical aspects of molecular flow simulation focussed on the development of collision models for diatomic and triatomic gases undergoing rotational excitation and energy exchange. A laser absorption technique is proposed and analyzed as sufficiently sensitive for measuring rotational populations in a nonequilibrium shock wave. Also reported is a survey of available literature on NO chemical rate data in preparation for experiments on NO decompensation kinetics in atmospheric pollution.

Kruger, C. H.↗

Heat transfer and forces on concave surfaces in free molecule flow.

A Monte Carlo modeling technique is described for mathematically simulating free molecular flows over a concave spherical surface and a concave cylindrical surface of finite length. The half-angle of the surfaces may vary from 0 to 90 degrees, and the incident flow may have an arbitrary speed ratio and an arbitrary angle of attack. Partial diffuse reflection and imperfect energy accommodation for molecules colliding with the surfaces are also considered. Results of heat transfer, drag and lift coefficients are presented for a variety of flow conditions. The present Monte Carlo results are shown to be in very good agreement with certain available theoretical solutions.

Fan, C.↗

A molecular model for ice nucleation and growth, attachment 2

The central field empirical pair potential model is applied to studying the effects of kinks, ledges, and vacancies on the absorption of water molecules from the vapor. Molecular dynamics simulations indicate that cluster and/or surface modes play a primary role in the absorption process, the flexibility of the hydrogen bond serves to decrease the energy required for structural interconversion, and the rapid distribution of added energy in a hydrogen bonded system lead to aggregate stability which greatly exceeds that predicted by static energy calculations.

Hale, B. N.↗

Algorithmic commonalities in the parallel environment

The ultimate aim of this project was to analyze procedures from substantially different application areas to discover what is either common or peculiar in the process of conversion to the Massively Parallel Processor (MPP). Three areas were identified: molecular dynamic simulation, production systems (rule systems), and various graphics and vision algorithms. To date, only selected graphics procedures have been investigated. They are the most readily available, and produce the most visible results. These include simple polygon patch rendering, raycasting against a constructive solid geometric model, and stochastic or fractal based textured surface algorithms. Only the simplest of conversion strategies, mapping a major loop to the array, has been investigated so far. It is not entirely satisfactory.

Mcanulty, Michael A.↗

Molecular dynamics computations of two dimensional supersonic rarefied gas flow past blunt bodies

This paper presents results of molecular dynamics computations of supersonic flow past a circular cylinder and past a flat plate perpendicular to a supersonic stream. The results are for Mach numbers of approximately 5 and 10, for several Knudsen numbers and several ratios of surface to free stream temperatures. A special feature of the computations is the use of relatively small numbers of particles in the molecular dynamics simulation, and an examination of the adequacy of using small numbers of particles to obtain physically useful results.

Greber, Isaac↗

Dynamics of charged particles in a Paul radio-frequency quadrupole trap

A molecular-dynamics simulation of hundreds of ions confined in a Paul trap has been performed. The simulation includes the trapped particles' micromotion and interparticle Coulomb interactions. A random walk in velocity was implemented to bring the secular motion to a given temperature which was numerically measured. When the coupling Gamma is large the ions from concentric shells which undergo a quadrupole oscillation at the RF frequency, while the ions within a shell form a 2D hexagonal lattice. Ion clouds at 5 mK show no RF heating for q(z) less than about 0.6, whereas rapid heating is seen for qz = 0.8.

Prestage, J. D.↗