Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DYNAMIC PROPERTY”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Molecular dynamics modeling of PPTA crystallite mechanical properties in the presence of defects

Here, the mechanical properties of PPTA crystallites, the fundamental building blocks of aramid polymer fibers such as Kevlar® and Twaron®, are studied here using molecular dynamics simulations. The ReaxFF interatomic potential is employed to study crystallite failure via covalent and hydrogen bond rupture in constant strain-rate tensile loading simulations. Emphasis is placed on analyzing how chain-end defects in the crystallite influence its mechanical response and fracture strength. Chain-end defects are found to affect the behavior of nearby chains in a region of the PPTA crystallite that is small relative to the typical crystallite size in manufactured aramid fibers. The central Csingle bondN bond along the backbone chain is identified as the weakest in the PPTA polymer chain backbone in dynamic strain-to-failure simulations of the crystallite. It is found that clustering of chain-ends leads to reduced crystallite strength and crystallite failure via hydrogen bond rupture and chain sliding, whereas randomly scattered defects impact the strength less and failure is by covalent bond rupture and chain scission. The axial crystallite modulus increases with increasing chain length and is independent of chain-end defect locations. On the basis of these findings, a theoretical model is proposed to predict the axial modulus as a function of chain length.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Application of the new IAPWS Guideline on the fast and accurate calculation of steam and water properties with the Spline-Based Table Look-Up Method (SBTL) in RELAP-7

The numerical simulation of thermalhydraulic processes in nuclear power plants requires very accurate and extremely fast algorithms for calculating the thermophysical properties of water and steam. In order to provide such algorithms, the International Association for the Properties of Water and Steam (IAPWS) has adopted the new “IAPWS Guideline on the Fast Calculation of Steam and Water Properties with the Spline-Based Table Look-Up Method (SBTL)”. In this article, the SBTL method is applied to property functions of specific volume and specific internal energy (v,e) based on the scientific formulation IAPWS-95 and the latest IAPWS formulations for transport properties. From the newly generated SBTL functions, thermodynamic and transport properties as well as their derivatives and inverse functions are calculable in the fluid range of state for pressures up to 100 MPa and for temperatures up to 1273 K, including the metastable liquid and the metastable vapor regions. The SBTL functions reproduce the underlying formulations with an accuracy of 10–100 ppm and significantly reduced computing times. The SBTL method has been implemented into the nuclear reactor system safety analysis code RELAP-7 [2] to consider the real fluid behavior of water and steam in a novel 7-equation two-phase flow model.

42 ENGINEERING↗
DOE OSTI2023

Hugoniot properties of porous stainless steel: Insights from molecular dynamics simulations

Here, we present a systematic study of Hugoniot properties of porous 316L stainless steel using both a simple interpolation scheme and direct shock simulations in order to probe pore collapse kinetics as well as final thermodynamic states. Both methods indicate that equilibrated Hugoniot properties depend on pore density only and not on the pore distribution or size. We then create a simple porous equation of state model that is shown to be accurate for a range of validation data. This allows us to extend our simulations to make direct comparison to experimental data that have generally significantly larger system sizes and durations. In addition, our direct shock simulations indicate that the relaxation time after hotspot formation is system size dependent and can reach nanosecond timescales for the largest pores investigated in our study, thereby possibly having a measurable effect on fast dynamic loading experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Here, using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

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

Splashing Droplets

Current data on droplet breakup is scarce for the sizes and velocities typical of practical applications such as in spray combustion processes and coating processes. While much more representative of practical applications, the small spatial scales and rapid time-scales prevent detailed measurement of the internal fluid dynamics and liquid property gradients produced by impinging upon surfaces. Realized through the extended spatial and temporal scales afforded by a microgravity environment, an improved understanding of drop breakup dynamics is sought to understand and ultimately control the impingement dynamics of droplets upon surfaces in practical situations. The primary objective of this research will be to mark the onset of different 'splashing modes' and to determine their temperature, pressure and angle dependence for impinging droplets representative of practical fluids. In addition, we are modeling the evolution of droplets that do not initially splash but rather undergo a 'fingering' evolution observed on the spreading fluid front and the transformation of these fingers into splashed products. An example of our experimental data is presented below. These images are of Isopar V impacting a mirror-polished surface. They were acquired using a high-speed camera at 1000 frames per second. They show the spreading of a single droplet after impact and ensuing finger instabilities. Normal gravity experimental data such as this will guide low gravity measurements in the 2.2 second drop tower and KC-135 aircraft as available. Presently we are in the process of comparing the experimental data of droplet shape evolution to numerical models, which can also capture the internal fluid dynamics and liquid property gradients such as produced by impingement upon a heated surface. To-date isothermal numerical data has been modeled using direct numerical simulations of representative splashing droplets. The data obtained so far indicates that the present model describes well the droplet wall interactions to a point in time just before splash. Additional information is included in the original extended abstract.

VanderWal, Randall L.↗
DOE OSTI2019

Future physics opportunities for high-density QCD at the LHC with heavy-ion and proton beams

The future opportunities for high-density QCD studies with ion and proton beams at the LHC are presented. Four major scientific goals are identified: the characterisation of the macroscopic long wavelength Quark–Gluon Plasma (QGP) properties with unprecedented precision, the investigation of the microscopic parton dynamics underlying QGP properties, the development of a unified picture of particle production and QCD dynamics from small (pp) to large (nucleus–nucleus) systems, the exploration of parton densities in nuclei in a broad ($x, Q^2$) kinematic range and the search for the possible onset of parton saturation. In order to address these scientific goals, high-luminosity Pb–Pb and p–Pb programmes are considered as priorities for Runs 3 and 4, complemented by high-multiplicity studies in pp collisions and a short run with oxygen ions. High-luminosity runs with intermediate-mass nuclei, for example Ar or Kr, are considered as an appealing case for extending the heavy-ion programme at the LHC beyond Run 4. The potential of the High-Energy LHC to probe QCD matter with newly-available observables, at twice larger center-of-mass energies than the LHC, is investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Dependence of Dynamic Modeling Accuracy on Sensor Measurements, Mass Properties, and Aircraft Geometry

The NASA Generic Transport Model (GTM) nonlinear simulation was used to investigate the effects of errors in sensor measurements, mass properties, and aircraft geometry on the accuracy of identified parameters in mathematical models describing the flight dynamics and determined from flight data. Measurements from a typical flight condition and system identification maneuver were systematically and progressively deteriorated by introducing noise, resolution errors, and bias errors. The data were then used to estimate nondimensional stability and control derivatives within a Monte Carlo simulation. Based on these results, recommendations are provided for maximum allowable errors in sensor measurements, mass properties, and aircraft geometry to achieve desired levels of dynamic modeling accuracy. Results using additional flight conditions and parameter estimation methods, as well as a nonlinear flight simulation of the General Dynamics F-16 aircraft, were compared with these recommendations

Grauer, Jared A.↗
NASA NTRS2018

Molecular Dynamics of ULTEM 9085 for 3D Manufacturing: Spectra, Thermodynamic Properties, and Shear Viscosity.

We present results of a molecular dynamic analysis of welding at the polymer-polymer interface. The analysis is performed for polyetherimide/ polycarbonate polymer blends. The work is motivated by the applications to 3D manufacturing in space. In the first part of the report, we discuss bulk and spectral characteristics of the amorphous polymer blends. The vibrational and infra-red spectra obtained using auto-correlation functions calculations in molecular dynamics are compared with the experimental spectra. The mechanical and thermal properties of the samples including heat capacity, bulk modulus, and thermal expansion coefficients are estimated and compared with experimental values. In the second part of the report, we discuss the result of molecular dynamical modeling of shear viscosity in a fully atomistic model of amorphous polymer blends with flat interface. The key result of the research is the demonstration of shear thinning behavior of the shear viscosity as a function of shear rate which is in good agreement with experimental data.

Luchinsky, Dmitry G.↗
DOE OSTI2016

Effect of temperature and geometric parameters on elastic properties of tungsten nanowire: A molecular dynamics study

Tungsten is a promising material and has potential use as battery anode. Tungsten nanowires are gaining attention from researchers all over the world for this wide field of application. In this paper, we investigated effect of temperature and geometric parameters (diameter and aspect ratio) on elastic properties of Tungsten nanowire. Aspect ratios (length to diameter ratio) considered are 8:1, 10:1, and 12:1 while diameter of the nanowire is varied from 1-4 nm. For 2 nm diameter sample (aspect ratio 10:1), temperature is varied (10 K ~ 1500 K) to observe elastic behavior of Tungsten nanowire under uniaxial tensile loading. EAM potential is used for molecular dynamic simulation. We applied constant strain rate of 10{sup 9} s{sup −1} to deform the nanowire. Elastic behavior is expressed through stress vs. strain plot. We also investigated the fracture mechanism of tungsten nanowire and radial distribution function. Investigation suggests peculiar behavior of Tungsten nanowire in nano-scale with double peaks in stress vs. strain diagram. Necking before final fracture suggests that actual elastic behavior of the material is successfully captured through atomistic modeling.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Effect of magnetic fields on the methyl rotation in a paramagnetic cobalt( ii ) complex. Quasielastic neutron scattering studies

Molecular dynamics is a fundamental property of metal complexes. These dynamic processes, especially for paramagnetic complexes under external magnetic fields, are in general not well understood. Quasielastic neutron scattering (QENS) in 0–4 T magnetic fields has been used to study the dynamics of Co(acac) 2 (D 2 O) 2 (1- d 4 , acac = acetylacetonate). At 80–100 K, rotation of the methyl groups on the acac ligands is the dominant dynamical process. This rotation is slowed down by the magnetic field increase. Rotation times at 80 K are 5.6(3) × 10 –10 s at 0 T and 2.04(10) × 10 –9 s at 4 T. The QENS studies suggest that methyl groups in these paramagnetic Co(II) molecules do not behave as isolated units, which is consistent with results from earlier magnetic susceptibility studies indicating the presence of intermolecular interactions. DFT calculations show that unpaired electron spin density in 1 is dispersed to the atoms of both acac and H 2 O ligands. Methyl torsions in 1- d 4 have also been observed at 5–100 K in inelastic neutron spectroscopy (INS). In conclusion, the QENS and INS results here help understand the dynamics of the compound in the solid state.

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

Star PolyMOCs with Diverse Structures, Dynamics, and Functions by Three-Component Assembly

We report star polymer metal–organic cage (polyMOC) materials whose structures, mechanical properties, functionalities, and dynamics can all be precisely tailored through a simple three-component assembly strategy. The star polyMOC network is composed of tetra-arm star polymers functionalized with ligands on the chain ends, small molecule ligands, and palladium ions; polyMOCs are formed via metal–ligand coordination and thermal annealing. The ratio of small molecule ligands to polymer-bound ligands determines the connectivity of the MOC junctions and the network structure. The use of large M12L24 MOCs enables great flexibility in tuning this ratio, which provides access to a rich spectrum of material properties including tunable moduli and relaxation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Characterizing the Small-Scale Dynamical, Ice Microphysical, and Residual Aerosol Properties of Mid-Latitude Cold Clouds: A Pilot Study (Final Report)

Our understanding of cirrus clouds has been advancing rapidly due to improvements in measurements made inside clouds, computer cloud models, and laboratory experimentation. This advancement in knowledge has made it clear that cirrus are complicated systems where the intricate interactions among dynamic air motions, solar and terrestrial radiation, ice crystal growth, and water and air chemistry can substantially influence cirrus structure and evolution. Measurements made inside clouds have shown that cirrus crystals are often complex with shapes that are often unclassifiable, especially at small sizes. Laboratory measurements have revealed that the ice crystal surface is complex, even at small sizes. Recent laboratory and satellite studies show that this surface complexity modifies the radiative properties of cirrus clouds in important ways. Moreover, observations and computer cloud modeling work show that the larger and smaller scale properties of cirrus are sensitive to the interplay between small-scale air motions and the aerosol articles upon which cirrus ice crystals form. Taken together, these studies indicate that accurately accounting for the radiative properties of cirrus clouds requires knowing the air motions and aerosol properties that influence the number concentration and crystal size, along with an understanding of the surface properties of vapor grown ice. Although this prior work has been illuminating, there have been few studies that provide detailed shape and surface complexity measurements for the small ice crystals that often compose cirrus clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Molecular Dynamics Simulation of the Structure and Properties of CaO-SiO{sub 2}-CaF{sub 2} Slag Systems

As a fundamental slag system in metallurgical processes, the CaO-SiO{sub 2}-CaF{sub 2} system, its properties, and related microstructures are of critical importance. The effect of the binary basicity and addition of CaF{sub 2} on the structure and viscosity of CaO-SiO{sub 2}-CaF{sub 2} slag were investigated using molecular dynamics simulations and Fourier transform infrared spectrometry. The average bond lengths of Si-O, Ca-O, and O-O were 1.61, 2.31, and 2.61 Å, respectively, and the coordinate number of Si remained constant. O-Si-O maintained an angle of 109.2, which was considered to be unaffected by the basicity changes. With increasing basicity, the network connectivity degrees Q{sup 3} and Q{sup 4} decreased, and the microstructure of the melt was significantly simpler. The effect of CaF{sub 2} on the depolymerization of the network structure was not obvious, and CaF{sub 2} acted mainly as a diluent. The calculations showed that the viscosity decreased upon setting up a ratio of nonbridging oxygen over tetrahedrally coordinated atoms (NBO/T). A good linear relationship was confirmed between the viscosity of the molten slag and the NBO/T.

36 MATERIALS SCIENCE↗