Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactions and forces”

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,045 records · Page 58

The effect of exhaust plume/afterbody interaction on installed scramjet performance

An upwind implicit Navier-Stokes code has been used to study hypersonic exhaust plume/afterbody flow fields. It is found that afterbody forces varied linearly with the nozzle exit pressure for moderately underexpanded jets, and that exhaust gases with low isentropic exponents (gamma) contribute up to 25 percent more force than high-gamma exhaust gases. Highly underexpanded jets are shown to create a strong plume shock, and the interaction of this shock with the afterbody produces a complicated pattern of crossflow separation.

Edwards, T. A.↗

SPARC-IGAC Symposium on Climate-Chemistry Interactions. Climate Feedback by Water Vapor in the Tropical Upper Troposphere

The strong greenhouse forcing by atmospheric water vapor is expected to play an important role in shaping the direction of any future changes in climate. We present calculations that provide a new perspective on the sensitivity of upper tropospheric water vapor to changes in surface temperature. Equilibrium states of our atmospheric model show unambiguously that as the surface warms, changes in the vertical distribution and temperature of detraining air parcels from tropical convection lead to higher water vapor mixing ratios in the upper troposphere. However, the increase in mixing ratio is not as large as the increase in saturation mixing ratio due to warmer environmental temperatures, so that the relative humidity decreases. Our analysis suggests that models that maintain a fixed relative humidity are likely overestimating the magnitude of the water vapor feedback.

Dessler, A. E.↗

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Substituent and Heteroatom Effects on π–π Interactions: Evidence That Parallel-Displaced π-Stacking is Not Driven by Quadrupolar Electrostatics

Stacking interactions are a recurring motif in supramolecular chemistry and biochemistry, where a persistent theme is a preference for parallel-displaced aromatic rings rather than face-to-face π-stacking. This is usually explained in terms of quadrupole–quadrupole interactions between the arene moieties but that interpretation is inconsistent with accurate calculations, which reveal that the quadrupolar picture is qualitatively wrong. At typical π-stacking distances, quadrupolar electrostatics may differ in sign from an exact calculation based on charge densities of the interacting arenes. We apply symmetry-adapted perturbation theory to dimers composed of substituted benzene and various aromatic heterocycles, which display a wide range of electrostatic interactions, and we investigate the interplay of Pauli repulsion, dispersion, and electrostatics as it pertains to parallel-displaced π-stacking. Profiles of energy components along cofacial slip-stacking coordinates support a prominent role for the “van der Waals model” (dispersion in competition with Pauli repulsion), even for polar monomers where electrostatic interactions are significant. While electrostatic interactions are necessary to explain the optimal face-to-face π-stacking distance and to account for the relative orientation of one polar arene with respect to another, we find no evidence to support continued invocation of quadrupolar electrostatics as a basis for π-stacking. Our results suggest that a driving force for offset-stacking exists even in the absence of electrostatic interactions. Consequently, tuning electrostatics via functionalization does not guarantee that slip-stacking can be avoided. This has implications for rational design of soft materials and other supramolecular architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scientific Objectives, Measurement Needs, and Challenges Motivating the PARAGON Aerosol Initiative

Aerosols are involved in a complex set of processes that operate across many spatial and temporal scales. Understanding these processes, and ensuring their accurate representation in models of transport, radiation transfer, and climate, requires knowledge of aerosol physical, chemical, and optical properties and the distributions of these properties in space and time. To derive aerosol climate forcing, aerosol optical and microphysical properties and their spatial and temporal distributions, and aerosol interactions with clouds, need to be understood. Such data are also required in conjunction with size-resolved chemical composition in order to evaluate chemical transport models and to distinguish natural and anthropogenic forcing. Other basic parameters needed for modeling the radiative influences of aerosols are surface reflectivity and three-dimensional cloud fields. This large suite of parameters mandates an integrated observing and modeling system of commensurate scope. The Progressive Aerosol Retrieval and Assimilation Global Observing Network (PARAGON) concept, designed to meet this requirement, is motivated by the need to understand climate system sensitivity to changes in atmospheric constituents, to reduce climate model uncertainties, and to analyze diverse collections of data pertaining to aerosols. This paper highlights several challenges resulting from the complexity of the problem. Approaches for dealing with them are offered in the set of companion papers.

optical depth↗

Ponderomotive recoil for electromagnetic waves

When waves damp or amplify on resonant particles in a plasma, nonresonant particles experience a recoil force that conserves the total momentum between particles and electromagnetic fields. This force is important to understand, as it can completely negate current drive and rotation drive mechanisms that are predicted on the basis of only resonant particles. Here, the existing electrostatic theory of this recoil force is extended to electromagnetic waves. While the result bears close similarity to historical fluid theories of laser–plasma interactions, it now incorporates both resonant and nonresonant particles, allowing momentum conservation to be self-consistently proven. Furthermore, the result is shown to be generally valid for kinetic plasmas, which is verified through single-particle hot-plasma simulations. Here, the new form of the force provides physical insight into the nature of the generalized Minkowski (plasmon) momentum of geometrical optics, which is shown to correspond to the momentum gained by the field and nonresonant particles as the wave is self-consistently ramped up from vanishing amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interannual Variability, Global Teleconnection, and Potential Predictability Associated with the Asian Summer Monsoon

In this Chapter, aspects of global teleconnections associated with the interannual variability of the Asian summer monsoon (ASM) are discussed. The basic differences in the basic dynamics of the South Asian Monsoon and the East Asian monsoon, and their implications on global linkages are discussed. Two teleconnection modes linking ASM variability to summertime precipitation over the continental North America were identified. These modes link regional circulation and precipitation anomalies over East Asia and continental North America, via coupled atmosphere-ocean variations over the North Pacific. The first mode has a large zonally symmetrical component and appears to be associated with subtropical jetstream variability and the second mode with Rossby wave dispersion. Both modes possess strong sea surface temperature (SST) expressions in the North Pacific. Results show that the two teleconnection modes may have its origin in intrinsic modes of sea surface temperature variability in the extratropical oceans, which are forced in part by atmospheric variability and in part by air-sea interaction. The potential predictability of the ASM associated with SST variability in different ocean basins is explored using a new canonical ensemble correlation prediction scheme. It is found that SST anomalies in tropical Pacific, i.e., El Nino, is the most dominant forcing for the ASM, especially over the maritime continent and eastern Australia. SST anomalies in the India Ocean may trump the influence from El Nino in western Australia and western maritime continent. Both El Nino, and North Pacific SSTs contribute to monsoon precipitation anomalies over Japan, southern Korea, northern and central China. By optimizing SST variability signals from the world ocean basins using CEC, the overall predictability of ASM can be substantially improved.

Lau, K. M.↗

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda–Devane–Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and non-polar moieties swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9–C9 interactions lead to the largest enhancement, whereas C9–C6OH and C6OH–C6OH interactions are less favorable for this bilayer system. As a result, ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating water adsorption in metal–organic frameworks with open metal sites using the 12-6-4 Lennard–Jones potential

Metal-organic frameworks (MOFs) with coordinatively unsaturated open metal sites (OMSs) are promising sorbent materials in atmospheric water harvesting (AWH) systems at low relative humidity (RH) due to their strong interactions with water molecules. However, accurate computational modelling of water adsorption in those materials are challenging, as standard force fields (FFs) based on the 12-6 Lennard-Jones (L-J) potential cannot properly describe the water-OMS interactions. In biomolecular simulations, the 12-6-4 L-J potential has been successfully used to model metal ion-water interactions in which the extra 1/r 4 -term is for charge-induced dipole electrostatics. In this work, we adopted this strategy and used the 12-6-4 L-J potential to model water-OMS interactions in MOFs for low RH AWH applications. Notably, without any modifications, the parameters used for aqueous metal ions were able to greatly improve the accuracy in predicting water adsorption isotherms in MOF-74, which highlights the simplicity and transferability of this method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Controls of interaction dynamics of orbital assembly

Building structures and spacecraft in orbit will require technologies for positioning, docking/berthing, and joining orbital structures. A fundamental problem underlying the operation of docking and berthing is that of controlling the contact dynamics of mechanical structures actuated by active mechanisms such as robotic devices. Control systems must be designed to control these active mechanisms so that both the free space motions and contact motions are stable and satisfy specifications on position accuracy and bounds on contact forces. For the large orbital structures of the future, the problem of interactive dynamics and control is fundamentally different in several ways than it was for spacecraft docking in the past. First, future space structures must be treated as flexible structures - the operations of docking, berthing, and assembly will need to respect the vibrations of the structures. Second, the assembly of these structures will require multiple-point contact, rather than the essentially single-point positioning of conventional spacecraft docking. Third, some assembly operations require the subassemblies to be brought and held in contact so that successful joining can be accomplished. A preliminary study of contact stability and compliance control design has resulted in the development of an analytical method and a design method to analyze stability. The analytical method analyzes the problem of stability when an actively-controlled structure contacts a passive structure. This method makes it possible to accurately estimate the stiffness of the passive structures with which the contact motion will become unstable. The analytic results suggest that passivity is neither achievable in practice, nor necessary as a design concept. A contact control system need only be passive up to a certain frequency; beyond that frequency the system can be stabilized with sufficiently small gains. With this concept the Center developed a design methodology for achieving desired compliant contact motions. This design method is based on H-infinity norm optimization, which makes it possible to consider both driving point mechanical impedance and systems robustness to modeling uncertainty. A laboratory facility was set up to verify experimentally the analytical and design theory.

Su, Renjeng↗

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

Solar-terrestrial research opportunities - A look to the future

Recent and proposed research programs and their impact on the understanding of solar-terrestrial interactions, the processes driving solar activity, and the propagation of the solar wind are outlined. The magnetospheres of the earth and major planets have been mapped and their interactions with the solar wind modeled. MHD activity has been identified as the driving force for stellar surface activity, and Pioneer data are permitting a mapping of the outer areas of the heliosphere. Further studies are required to understand solar magnetic field-earth magnetosphere, ionosphere-atmosphere interactions, the effects of solar heating on earth atmosphere ozone levels, and the effects of solar transient events on the earth atmosphere. The initiation of an International Geosphere-Biosphere Program is now under study and several satellites for solar-terrestrial research are under development as cooperative efforts among the U.S., Japanese, and European agencies.

Intriligator, D. S.↗

NbTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

Impact of solvent forces and broken symmetry on the assembly of designed proteins at a liquid-solid interface

The era of protein design has enabled the creation of hybrid protein-inorganic interfaces, leading to both surface-directed self-assembly of de novo protein architectures and protein-directed formation of inorganic materials. However, the resulting patterns of protein assembly are often unexpected, implying that essential interactions are not accounted for in current design platforms. Here, we use high-speed atomic force microscopy (AFM) analyzed through machine learning to follow the assembly of protein nanorods in aqueous electrolytes on two types of mica exhibiting disparate symmetry elements, which are imprinted on the overlying hydration structure. Using Monte Carlo simulations, we reproduce the observed phases and show that an observed smectic phase, previously thought to be unstable for non-interacting rods in two dimensions, emerges when crystal symmetry introduces a directional bias. The findings demonstrate the importance of incorporating solvent forces as modulated by the hydration structure inherent to interfacial systems when designing protein assemblies at liquid-crystal interfaces. Coupling physics-based simulations that can account for these factors to de novo protein design algorithms can lead to improved design platforms for bio-inspired, hybrid materials.

bioinspired materials↗

A machine learning inversion scheme for determining interaction from scattering

Small angle scattering techniques have now been routinely used to quantitatively determine the potential of mean force in colloidal suspensions. However the numerical accuracy of data interpretation is often compounded by the approximations adopted by liquid state analytical theories. To circumvent this long standing issue, here we outline a machine learning strategy for determining the effective interaction in the condensed phases of matter using scattering. Via a case study of colloidal suspensions, we show that the effective potential can be probabilistically inferred from the scattering spectra without any restriction imposed by model assumptions. Comparisons to existing parametric approaches demonstrate the superior performance of this method in accuracy, efficiency, and applicability. This method can effectively enable quantification of interaction in highly correlated systems using scattering and diffraction experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nucleon-Nucleon Scattering Up to N5LO in Chiral Effective Field Theory

During the past few decades a large effort has been made toward describing the NN interaction in the framework of chiral Effective Field Theory (EFT). The main idea is to exploit the symmetries of QCD to obtain an effective theory for low energy nuclear systems. In 2003, the first accurate charge-dependent NN potential in this scheme was developed and it has been applied to many ab-initio calculations, opening the possibility to study nuclear systems in a systematic and accurate way. It was shown that the fourth order (N3LO) was necessary and sufficient to describe the NN scattering data with a χ 2 /d.o.f on the order of so-called high precision potentials. However the systematics of chiral EFT also allow to relate two- and many-body interactions in a well-defined way. Since many-body forces make their first appearance at higher order, they are substantially smaller than their two-body counterparts, but may never-the-less be crucial for some processes. Thus, there are observables where they can have a big impact and, for example, there are indications that they solve the long standing Ay puzzle of N-d scattering. The last few years, have also seen substantial progress toward higher orders of chiral EFT which was motivated by the fact that only three-body forces of rather high order may solve some outstanding issues in microscopic nuclear structure and reactions. In this chapter we will review the latest contributions of the authors to development of chiral EFT based potentials up to N4LO as well as first calculations conducted for NN scattering at N5LO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗