Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interaction force”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Compact two-electron wave function for bond dissociation and Van der Waals interactions: A natural amplitude assessment

Electron correlations in molecules can be divided in short range dynamical correlations, long range Van der Waals type interactions, and near degeneracy static correlations. In this work, we analyze for a one-dimensional model of a two-electron system how these three types of correlations can be incorporated in a simple wave function of restricted functional form consisting of an orbital product multiplied by a single correlation function f (r{sub 12}) depending on the interelectronic distance r{sub 12}. Since the three types of correlations mentioned lead to different signatures in terms of the natural orbital (NO) amplitudes in two-electron systems, we make an analysis of the wave function in terms of the NO amplitudes for a model system of a diatomic molecule. In our numerical implementation, we fully optimize the orbitals and the correlation function on a spatial grid without restrictions on their functional form. Due to this particular form of the wave function, we can prove that none of the amplitudes vanishes and moreover that it displays a distinct sign pattern and a series of avoided crossings as a function of the bond distance in agreement with the exact solution. This shows that the wave function ansatz correctly incorporates the long range Van der Waals interactions. We further show that the approximate wave function gives an excellent binding curve and is able to describe static correlations. We show that in order to do this the correlation function f (r{sub 12}) needs to diverge for large r{sub 12} at large internuclear distances while for shorter bond distances it increases as a function of r{sub 12} to a maximum value after which it decays exponentially. We further give a physical interpretation of this behavior.

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

A new smoothing function to introduce long-range electrostatic effects in QM/MM calculations

A new method to account for long range electrostatic contributions is proposed and implemented for quantum mechanics/molecular mechanics long range electrostatic correction (QM/MM-LREC) calculations. This method involves the use of the minimum image convention under periodic boundary conditions and a new smoothing function for energies and forces at the cutoff boundary for the Coulomb interactions. Compared to conventional QM/MM calculations without long-range electrostatic corrections, the new method effectively includes effects on the MM environment in the primary image from its replicas in the neighborhood. QM/MM-LREC offers three useful features including the avoidance of calculations in reciprocal space (k-space), with the concomitant avoidance of having to reproduce (analytically or approximately) the QM charge density in k-space, and the straightforward availability of analytical Hessians. The new method is tested and compared with results from smooth particle mesh Ewald (PME) for three systems including a box of neat water, a double proton transfer reaction, and the geometry optimization of the critical point structures for the rate limiting step of the DNA dealkylase AlkB. As with other smoothing or shifting functions, relatively large cutoffs are necessary to achieve comparable accuracy with PME. For the double-proton transfer reaction, the use of a 22 Å cutoff shows a close reaction energy profile and geometries of stationary structures with QM/MM-LREC compared to conventional QM/MM with no truncation. Geometry optimization of stationary structures for the hydrogen abstraction step by AlkB shows some differences between QM/MM-LREC and the conventional QM/MM. These differences underscore the necessity of the inclusion of the long-range electrostatic contribution.

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

Implementation of a generalized actuator line model for wind turbine parameterization in the Weather Research and Forecasting model

A generalized actuator line (GAL) wind turbine parameterization is implemented within the Weather Research and Forecasting model to enable high-fidelity large-eddy simulations of wind turbine interactions with boundary layer flows under realistic atmospheric forcing conditions. Numerical simulations using the GAL parameterization are evaluated against both an already implemented generalized actuator disk (GAD) wind turbine parameterization and two field campaigns that measured the inflow and near-wake regions of a single turbine. The representation of wake wind speed, variance, and vorticity distributions is examined by comparing fine-resolution GAL and GAD simulations and GAD simulations at both fine and coarse-resolutions. The higher-resolution simulations show slightly larger and more persistent velocity deficits in the wake and substantially increased variance and vorticity when compared to the coarse-resolution GAD. The GAL generates distinct tip and root vortices that maintain coherence as helical tubes for approximately one rotor diameter downstream. Coarse-resolution simulations using the GAD produce similar aggregated wake characteristics to both fine-scale GAD and GAL simulations at a fraction of the computational cost. Furthermore, the GAL parameterization provides the capability to resolve near wake physics, including vorticity shedding and wake expansion.

17 WIND ENERGY↗
DOE OSTI2021

First measurement of the asymmetry and the Gerasimov-Drell-Hearn integrand from the 3 $\vec{H}$e ($\vec{γ},p$) 2 H reaction at an incident photon energy of 29 MeV

The first measurement of the 3 $\vec{H}$e ($\vec{γ},p$) 2 H process was performed at the High Intensity γ-ray Source ( HIγS ) facility at Triangle Universities Nuclear Laboratory using a circularly polarized, monoenergetic γ -ray beam and a longitudinally polarized 3 He target. The spin-dependent asymmetry and the contribution from the two-body photodisintegration to the 3 He Gerasimov-Drell-Hearn integrand are extracted and compared with state-of-the-art three-nucleon system calculations at the incident photon energy of 29 MeV. The data are in general agreement with the various theoretical predictions based on the Siegert theorem or on explicit inclusion of meson-exchange currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory (χ EFT) that are provided by bound-state observables in the A = 3 and A = 4 sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the χ EFT expansion at next-to-next-to-leading order. A consistent solution for the 3 H binding energy, the 4 He binding energy and radius, and the 3 H β-decay rate can only be obtained if χ EFT truncation errors are included in the analysis. Here, the β-decay rate is the only one of these that yields a nondegenerate constraint on the 3NF low-energy constants, which makes it crucial for the parameter estimation. We use eigenvector continuation for fast and accurate emulation of no-core shell model calculations of the few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the parameters that quantify the truncation error. We find a χ EFT expansion parameter of Q = 0.33 ± 0.06 for these observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

A predictive theory for elastic scattering and recoil of protons from 4 He

Low-energy cross sections for elastic scattering and recoil of protons from 4 He nuclei (also known as α particles) are calculated directly by solving the Schrodinger equation for five nucleons interacting through accurate two- and three-nucleon forces derived within the framework of chiral effective field theory. Precise knowledge of these processes at various proton backscattering/recoil angles and energies is needed for the ion-beam analysis of numerous materials, from the surface layers of solids, to thin films, to fusion-reactor materials. Indeed, the same elastic scattering process, in two different kinematic configurations, can be used to probe the concentrations and depth profiles of either hydrogen or helium. Furthermore, we compare our results to available experimental data and show that direct calculations with modern nuclear potentials can help to resolve remaining inconsistencies among data sets and can be used to predict these cross sections when measurements are not available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fire as a fundamental ecological process: Research advances and frontiers

Fire is a powerful ecological and evolutionary force that regulates organismal traits, population sizes, species interactions, community composition, carbon and nutrient cycling and ecosystem function. It also presents a rapidly growing societal challenge, due to both increasingly destructive wildfires and fire exclusion in fire–dependent ecosystems. As an ecological process, fire integrates complex feedbacks among biological, social and geophysical processes, requiring coordination across several fields and scales of study. Here, we describe the diversity of ways in which fire operates as a fundamental ecological and evolutionary process on Earth. We explore research priorities in six categories of fire ecology: (a) characteristics of fire regimes, (b) changing fire regimes, (c) fire effects on above–ground ecology, (d) fire effects on below–ground ecology, (e) fire behaviour and (f) fire ecology modelling. We identify three emergent themes: the need to study fire across temporal scales, to assess the mechanisms underlying a variety of ecological feedbacks involving fire and to improve representation of fire in a range of modelling contexts. As fire regimes and our relationships with fire continue to change, prioritizing these research areas will facilitate understanding of the ecological causes and consequences of future fires and rethinking fire management alternatives.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Studies of Dynamic, Radiative Macroscopic Magnetized HED Plasmas with Closed B-Field Lines

The purpose of this research has been to study the physics of macroscopic magnetized high-energy-density laboratory plasmas (HEDLPs) created through the compression of a high-beta compact toroid (CT) plasma having closed magnetic field lines. The high-beta CT chosen for this work is a field-reversed configuration (FRC). The basic approach is to investigate CT plasmas as they are compressed to a HED state by the electromagnetic implosion of a surrounding metallic shell or solid liner (Figure 1). The shell provides an axisymmetric, electrically-conducting boundary around the plasma and its supporting magnetic field and is imploded by means of the magnetic pressure force arising from axial current flow in the liner interacting with its associated azimuthal magnetic field. Compression of the CT will bring the plasma to fusion temperatures at higher densities and magnetic fields (multi-MegaGauss [MG]) than have previously been present in conventional magnetic fusion approaches. The resulting energy densities will be ~1 Mbar or greater and thus will place the plasma in a parameter space intermediate to MFE and IFE. This work has been a collaboration between the Air Force Research Laboratory, Los Alamos National Laboratory, and NumerEx, LLC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

LAMMPS strong scaling performance optimization on Blue Gene/Q

LAMMPS "Large-scale Atomic/Molecular Massively Parallel Simulator" is an open-source molecular dynamics package from Sandia National Laboratories. Significant performance improvements in strong-scaling and time-to-solution for this application on IBM's Blue Gene/Q have been achieved through computational optimizations of the OpenMP versions of the short-range Lennard-Jones term of the CHARMM force field and the long-range Coulombic interaction implemented with the PPPM (particle-particle-particle mesh) algorithm, enhanced by runtime parameter settings controlling thread utilization. Additionally, MPI communication performance improvements were made to the PPPM calculation by re-engineering the parallel 3D FFT to use MPICH collectives instead of point-to-point. Performance testing was done using an 8.4-million atom simulation scaling up to 16 racks on the Mira system at Argonne Leadership Computing Facility (ALCF). Speedups resulting from this effort were in some cases over 2x.

Coffman, Paul↗
DOE OSTI2016

Improving the Representation of Soluble Iron in Climate Models

Mineral dust produced in the arid and semi-arid regions of the world is the dominant source of iron (Fe) in atmospheric aerosol inputs to the open ocean. The bioavailable Fe fraction of atmospheric dust is thought to regulate and occasionally limit the primary productivity in large oceanic regions, which influences the CO2 uptake from the atmosphere affecting the Earth’s climate. Because Fe bioavailability cannot be directly measured, it is assumed that the dissolved Fe or highly reactive Fe in the dust is bioavailable. The fraction of soluble Fe in dust is mainly controlled by: (1) the mineral composition of the soils and the emitted dust from the source areas; (2) the atmospheric processing that converts the Fe in Fe-bearing minerals into highly soluble forms of Fe. The project has mainly focused on constraining the mineral composition of dust aerosols (1), a previously neglected, yet a key issue to constrain the deposition of soluble iron. Deriving aerosol mineral composition requires global knowledge of the soil mineral content, which is available from poorly constrained global atlases. In addition, the mineral content of the emitted aerosol differs from that of the parent soil. Measurements of soil mineral fractions are based upon wet sedimentation (or ’wet sieving’) techniques that disturb the soil sample, breaking aggregates that are found in the original, undispersed soil that is subject to wind erosion. Wet sieving alters the soil size distribution, replacing aggregates that are potentially mobilized as aerosols with a collection of smaller particles. A major challenge is to derive the size-distributed mineral fractions of the emitted dust based upon their fractions measured from wet-sieved soils. Finally, representations of dust mineral composition need to account for mixtures of minerals. Examination of individual particles shows that iron, an element that is central to many climate processes, is often found as trace impurities of iron oxide attached to aggregates of other minerals. This is another challenge that has been tackled by the project. The project has produced a major step forward on our understanding of the key processes needed to predict the mineral composition of dust aerosols by connecting theory, modeling and observations. The project has produced novel semi-empirical and theoretical methods to estimate the emitted size distribution and mineral composition of dust aerosols. These methods account for soil aggregates that are potentially emitted from the original undisturbed soil but are destroyed during wet sieving. The methods construct the emitted size distribution of individual minerals building upon brittle fragmentation theory, reconstructions of wet-sieved soil mineral size distributions, and/or characteristic mineral size distributions estimated from observations at times of high concentration. Based on an unprecedented evaluation with a new global compilation of observations produced with the project support, we showed that the new methods remedy some key deficiencies compared to the previous state-of-the-art. This includes the correct representation of Fe-bearing phyllosilicates at silt sizes, where they are abundant according to observations. In addition, the quartz fraction of silt particles is in better agreement with measured values. In addition, we represent an additional class of iron oxide aerosol that is a small impurity embedded within other minerals, allowing it to travel farther than in its pure crystalline state. We assume that these impurities are least frequent in soils rich in iron oxides (as a result of the assumed effect of weathering that creates pure iron oxide crystals). The mineral composition of dust is also important to other interaction with climate - through shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, and the heterogeneous formation of sulfates and nitrates - and to its impacts upon human health. Despite the importance of mineral composition, models have typically assumed that soil dust aerosols have globally uniform composition. The results of this project will allow an improved estimation of the dust effects upon climate and health.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Development and Rigorous Validation of Modules for Predicting Organic Aerosol CCN Activity from Molecular Properties

The indirect interaction of anthropogenic aerosols with clouds (aerosol indirect forcing) represents one of the main uncertainties in climate projections, and reducing this uncertainty requires improved understanding of the aerosol life cycle, including atmospheric aging and the ability to nucleate cloud droplets or crystals. A key parameter that links the chemical properties of aerosol to its nucleating properties, and ultimately to aerosol indirect forcing, is its hygroscopicity (ability to uptake water). One way to compare and systematize the ability of various chemical compounds or particle types to nucleate cloud drops is via the hygroscopicity parameter κ, which can be used directly in global models to simulate aerosol indirect forcing. In this project we further investigated the relationships between the chemical composition of organic aerosol particles, which are known to be important components of atmospheric particulate matter, and their ability to serve as cloud condensation nuclei (CCN). In particular, we isolated the role of viscosity on κ by applying methods similar to those we used previously in our DOE/ASR research. Our motivation for the focus on viscosity was through its link to particle phase state. Particle phase affects water uptake on timescales relevant to the laboratory and also likely relevant to atmospheric processes, with reduced apparent hygroscopicities for semisolid and solid particles, compared with liquid particles. Water uptake itself affects viscosity, an important feedback that we explored in this project. The objectives of the project were as follows: (1) measure the effects of molecular structure, temperature, and relative humidity on the viscosity of organic aerosol formed via controlled laboratory experiments, (2) measure the effects of organic aerosol viscosity on water uptake and CCN activity above and below water saturation, (3) further develop the κ parameterization of organic aerosol hygroscopicity to include viscosity/phase, and (4) experimentally validate the κ parameterization for organic aerosol with chemical composition ranging from relatively simple to complex. It was proposed that particle viscosity exerts a kinetic control on equilibrium water uptake. High viscosity was expected to delay water uptake based on diffusional limitations. This would manifest itself in reduced observed CCN activity or κ. Understanding the interplay between water content and viscosity (composition control on particle viscosity) and between viscosity and equilibrium water content (diffusional control on water uptake) were the overall science objectives of this grant. Targeted experiments were carried out to test the kinetic control hypothesis. These experiments included studies that explicitly measured the equilibration time scale of water uptake, and experiments that tested the role of viscous particle coatings on the observed CCN activity. From the combined work, it is concluded that mass transfer limitations by glassy organic particles are unlikely to affect cloud droplet activation near laboratory temperatures. This implies that viscosity does not need to be included in parameterizations that link particle chemical composition and κ. A parameterization to predict κ through the UNIFAC group contribution method was developed with previous ASR funding. As planned, this computational module was further developed and validated here through targeted laboratory studies. The computational module has been made available as free software under the GNU public license. A database of κ-values to aid model development has been compiled and is available online via an open access model. In conducting this work, significant effort was made to understand how aerosol composition, including the associated aerosol water content, affects viscosity. This included a study detailing the influence of functional group composition on organic aerosol viscosity. We identified a strong link between volatility and viscosity. Functional groups responsible for low vapor pressure also promote high viscosity. Specifically, carboxylic acid and hydroxyl groups were identified as major contributors to high viscosity. Several studies were carried out to characterize the link between composition and the glass-to-liquid phase transition. Compositions tested include models of organic/inorganic particles and models of complex organic aerosol formed through secondary processes. A simple computational model to predict the amorphous phase diagram was developed and tested. Excellent agreement between model and measurements was obtained for organic/inorganic mixtures. The computational model is simple enough to be included in global model simulations. However, constraints of input parameters such as the Gordon-Taylor constant for binary mixtures for a wide range of complex mixtures have not yet been fully evaluated. Combined, the performed work provides new constraints on the amorphous phase state of organic aerosols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗
DOE OSTI2017

Planet-induced Stellar Pulsations in HAT-P-2's Eccentric System

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet’s atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet–star interactions in HAT-P-2's eccentric planetary system gained from the analysis of ∼350 hr of 4.5 μ m observations with the Spitzer Space Telescope . The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b’s 4.5 μ m photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet’s orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

79 ASTRONOMY AND ASTROPHYSICS↗