SEARCH · Engineering Papers
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.
Speech understanding among researchers and managers, current developments in voice technology, and an exchange of information concerning government voice technology efforts are discussed.
Surfactant-induced Marangoni effects strongly alter the stresses exerted along fluid particle interfaces. In low gravity processes, these stresses can dictate the system behavior. The dependence of Marangoni effects on surfactant physical chemistry is not understood, severely impacting our ability to predict and control fluid particle flows. A droplet in an extensional flow allows the controlled study of stretching and deforming interfaces. The deformations of the drop allow both Marangoni stresses, which resist tangential shear, and Marangoni elasticities, which resist surface dilatation, to develop. This flow presents an ideal model system for studying these effects. Prior surfactant-related work in this flow considered a linear dependence of the surface tension on the surface concentration, valid only at dilute surface concentrations, or a non-linear framework at concentrations sufficiently dilute that the linear approximation was valid. The linear framework becomes inadequate for several reasons. The finite dimensions of surfactant molecules must be taken into account with a model that includes surfaces saturation. Nonideal interactions between adsorbed surfactant molecules alter the partitioning of surfactant between the bulk and the interface, the dynamics of surfactant adsorptive/desorptive exchange, and the sensitivity of the surface tension to adsorbed surfactant. For example, cohesion between hydrocarbon chains favors strong adsorption. Cohesion also slows the rate of desorption from interfaces, and decreases the sensitivity of the surface tension to adsorbed surfactant. Strong cohesive interactions result in first order surface phase changes with a plateau in the surface tension vs surface concentration. Within this surface concentration range, the surface tension is decoupled from surface concentration gradients. We are engaged in the study of the role of surfactant physical chemistry in determining the Marangoni stresses on a drop in an extensional flow in a numerical and experimental program. Using surfactants whose dynamics and equilibrium behavior have been characterized in our laboratory, drop deformation will be studied in ground-based experiment. In an accompanying numerical study, predictive drop deformations will be determined based on the isotherm and equation of state determined in our laboratory. This work will improve our abilities to predict and control all fluid particle flows.
OVERFLOW, a widely used CFD simulation code, is combined with a visualization system, pV3, to experiment with an environment for simulation/visualization co-processing on a SGI Origin 2000 computer(O2K) system. The shared memory version of the solver is used with the O2K 'pfa' preprocessor invoked to automatically discover parallelism in the source code. No other explicit parallelism is enabled. In order to study the scaling and performance of the visualization co-processing system, sample runs are made with different processor groups in the range of 1 to 254 processors. The data exchange between the visualization system and the simulation system is rapid enough for user interactivity when the problem size is small. This shared memory version of OVERFLOW, with minimal parallelization, does not scale well to an increasing number of available processors. The visualization task takes about 18 to 30% of the total processing time and does not appear to be a major contributor to the poor scaling. Improper load balancing and inter-processor communication overhead are contributors to this poor performance. Work is in progress which is aimed at obtaining improved parallel performance of the solver and removing the limitations of serial data transfer to pV3 by examining various parallelization/communication strategies, including the use of the explicit message passing.
We present the 2335-2325 cm(exp -1) infrared spectra and band positions, profiles, and strengths (A values) of solid nitrogen and binary mixtures of N2 with other molecules at 12 K. The data demonstrate that the strength of the infrared forbidden N2 fundamental near 2328 cm(exp -1) is moderately enhanced in the presence of NH3, strongly enhanced in the presence of H2O and very strongly enhanced in the presence of CO2, but is not significantly affected by CO, CH4, or O2. The mechanisms for the enhancements in N2-NH3 and N2-H2O mixtures are fundamentally different from those proposed for N2-CO2 mixtures. In the first case, interactions involving hydrogen-bonding are likely the cause. In the latter, a resonant exchange between the N2 stretching fundamental and the O-18=C-12 asymmetric stretch of O-18C-12O-16 is indicated. The implications of these results for several astrophysical issues are briefly discussed.
We present the 2335-2325 cm-1 infrared spectra and band positions, profiles and strengths (A values) of solid nitrogen and binary mixtures of N2 with other molecules at 12 K. The data demonstrate that the strength of the infrared forbidden N2 fundamental near 2328 cm-1 is moderately enhanced in the presence of NH3, strongly enhanced in the presence of H2O and very strongly enhanced (by over a factor of 1000) in the presence of CO2, but is not significantly affected by CO, CH4, or O2. The mechanisms for the enhancements in N2-NH3 and N2-H2O mixtures are fundamentally different from those proposed for N2-CO2 mixtures. In the first case, interactions involving hydrogen-bonding are likely the cause. In the latter, a resonant exchange between the N2 stretching fundamental and the 18O = 12C asymmetric stretch of 18O12C16O is indicated. The implications of these results for several astrophysical issues are briefly discussed.
The process of magnetic reconnection plays an important role during the interaction of the solar wind with the Earth's magnetosphere which leads to the exchange of mass, momentum, and energy between these two highly conducting plasmas.
The first “Changing Microbiomes” conference was successfully held at the Pennsylvania State University from May 31st – June 3rd, 2022. The overall goal of the symposium was to provide space for researchers to connect, exchange knowledge, and partner to advance our understanding of the mechanisms and applications of ecological and evolutionary changes in microbiomes across systems. The symposium spanned terrestrial and aquatic contexts and host-associated and free-living microbial communities. As such, the symposium will enable cross-disciplinary knowledge and idea exchange to promote future interdisciplinary collaborations. In addition, the program allowed formal and informal interactions between and among students, postdocs, PIs, and industry scientists.
For the first time, the (d, 2 He) reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the β + direction from an unstable nucleus. The new technique was made possible by the use of an active-target time-projection chamber and a magnetic spectrometer, and opens a path to addressing a range of scientific challenges, including in astrophysics and neutrino physics. Here, in this Letter, the nucleus studied was 14 O, and the Gamow-Teller transition strength to 14 N was extracted up to an excitation energy of 22 MeV. The data were compared to shell-model and state-of-the-art coupled-cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. In a significant step forward to better understand this quenching, the coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents.
In order to understand how the climate responds to variations in forcing, one necessary component is to understand the full distribution of variability of exchanges of heat and moisture between the atmosphere and ocean. Surface heat and moisture fluxes are critical to the generation and decay of many coupled air-sea phenomena. These mechanisms operate across a number of scales and contain contributions from interactions between the anomalous (i.e. non-mean), often extreme-valued, flux components. Satellite-derived estimates of the surface turbulent and radiative heat fluxes provide an opportunity to assess results from modeling systems. Evaluation of only time mean and variability statistics, however only provides limited traceability to processes controlling what are often regime-dependent errors. This work will present an approach to evaluate the representation of the turbulent fluxes at the air-sea interface in the current and evolving Goddard Earth Observing System (GEOS) model. A temperature and moisture vertical profile-based clustering technique is used to identify robust weather regimes, and subsequently intercompare the turbulent fluxes and near-surface parameters within these regimes in both satellite estimates and GEOS-driven data sets. Both model reanalysis (MERRA) and seasonal-to-interannual coupled GEOS model simulations will be evaluated. Particular emphasis is placed on understanding the distribution of the fluxes including extremes, and the representation of near-surface forcing variables directly related to their estimation. Results from these analyses will help identify the existence and source of regime-dependent biases in the GEOS model ocean surface turbulent fluxes. The use of the temperature and moisture profiles for weather-state clustering will be highlighted for its potential broad application to 3-D output typical of model simulations.
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
The electron bombardment ion thruster has been under development for a number of years and during this time, studies of the plasmas produced by the thrusters and their interactions with spacecraft have been evaluated, based on available data. Due to diagnostic techniques used and facility effects, there is uncertainty as to the reliability of data from these early studies. This paper presents data on the flow of the charge-exchange plasma produced just downstream of the thruster's ion optics. The 'end-effect' of a cylindrical Langmuir probe is used to determine ion density and directed ion velocity. Results are compared with data obtained from a retarding potential analyzer-Faraday cup.
A general computer code was developed to allow calculation of atom-atom and ion-atom transport collision integrals from accurate potential energy curves described by a set of discrete data points for a broad range of scattering conditions. This code is based upon semiclassical approximations that properly account for quantum mechanical behavior such as tunneling effects near a barrier maximum, resonance charge exchange, and nuclear symmetry effects. Transport collision integrals were determined for N-N, O-O, N(+)-N, and O(+)-O interactions from complete sets of accurate potential functions derived from combined experimental and ab initio structure calculations. For the O-O case, this includes results for excited states. The calculated values of the N(+)-N and O(+)-O resonance charge exchange cross section Q(ex) agree well with measurements from beam experiment that are available at high energies where the diffusion cross section Q(d) satisfies Q(d) approximately equal to 2Q(ex).
The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.
The shape of the velocity distribution of water group ions observed by the Giotto ion mass spectrometer on its approach to comet Halley is modeled to derive empirical values for the rates of ionization, energy diffusion, and loss in the midcometosheath. The model includes the effect of rapid pitch angle scattering into a bispherical shell distribution as well as the effect of the magnetization of the plasma on the charge exchange loss rate. It is found that the average rate of ionization of cometary neutrals in this region of the cometosheath appears to be of the order of a factor 3 faster than the `standard' rates approx. 1 x 10(exp -6)/s that are generally assumed to model the observations in most regions of the comet environment. For the region of the coma studied in the present work (approx. 1 - 2 x 10(exp 5) km from the nucleus), the inferred energy diffusion coefficient is D(sub 0) approx. equals 0.0002 to 0.0005 sq km/cu s, which is generally lower than values used in other models. The empirically obtained loss rate appears to be about an order of magnitude greater than can be explained by charge exchange with the `standard' cross section of approx. 2 x 10(exp -15)sq cm. However such cross sections are not well known and for water group ion/water group neutral interactions, rates as high as 8 x 10(exp -15) sq cm have previously been suggested in the literature. Assuming the entire loss rate is due to charge exchange yields a rate of creation of fast neutral atoms of the order of approx. 10(exp -4)/s or higher, depending on the level of velocity diffusion. The fast neutrals may, in turn, be partly responsible for the higher-than-expected ionization rate.
We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). 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 Zr. 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.
Topological surface states (TSS) in the prototypical topological insulator (TI) Bi 2 Se 3 are frequently characterized using optical probes, but electron-hole interactions and their effect on surface localization and optical response of the TSS remain unexplored. Here, we use ab initio calculations to understand excitonic effects in the bulk and surface of Bi 2 Se 3 . We identify multiple series of chiral excitons that exhibit both bulk and TSS character, due to exchange-driven mixing. Further, our results address fundamental questions about the degree to which electron-hole interactions can relax the topological protection of surface states and dipole selection rules for circularly polarized light in TIs by elucidating the complex intermixture of bulk and surface states excited in optical measurements and their coupling to light.
We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). 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 Ta. 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 sites. 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,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.
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.