Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Tensor Force Manifestations in Ab Initio Study of the {sup 2}H(d,{gamma}){sup 4}He, {sup 2}H(d,p){sup 3}H, and {sup 2}H(d,n){sup 3}He Reactions

The {sup 2}H(d,p){sup 3}H, {sup 2}H(d,n){sup 3}He, and {sup 2}H(d,{gamma}){sup 4}He reactions are studied at low energies in a multichannel ab initio model that takes into account the distortions of the nuclei. The internal wave functions of these nuclei are given by the stochastic variational method with the AV8{sup '} realistic interaction and a phenomenological three-body force included to reproduce the two-body thresholds. The obtained astrophysical S factors are all in very good agreement with the experiment. The most important channels for both transfer and radiative capture are identified by comparing to calculations with an effective central force. They are all found to dominate thanks to the tensor force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Quantum mechanical force field for hydrogen fluoride with explicit electronic polarization

The explicit polarization (X-Pol) theory is a fragment-based quantum chemical method that explicitly models the internal electronic polarization and intermolecular interactions of a chemical system. X-Pol theory provides a framework to construct a quantum mechanical force field, which we have extended to liquid hydrogen fluoride (HF) in this work. The parameterization, called XPHF, is built upon the same formalism introduced for the XP3P model of liquid water, which is based on the polarized molecular orbital (PMO) semiempirical quantum chemistry method and the dipole-preserving polarization consistent point charge model. We introduce a fluorine parameter set for PMO, and find good agreement for various gas-phase results of small HF clusters compared to experiments and ab initio calculations at the M06-2X/MG3S level of theory. In addition, the XPHF model shows reasonable agreement with experiments for a variety of structural and thermodynamic properties in the liquid state, including radial distribution functions, interaction energies, diffusion coefficients, and densities at various state points.

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

DEM simulation of experimental dense granular packing

In this study we present numerical analysis performed on the experimental results of sphere packings of mono-sized hard sphere whose packing fraction spans across a wide range of 0.59<{Phi}<0.72. Using X-ray Computed Tomography (XCT), we have full access to the 3D structure of the granular packings. Numerical analysis performed on thr data provides the first experimental proofs of how densification affects local order parameters. Furthermore by combining Discrete Element Method (DEM) and the experimental results from XCT, we investigate how the intergranular forces change with the onset of crystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deuteron-α scattering: Separable versus nonseparable Faddeev approach

Deuteron-induced reactions are widely used to probe nuclear structure and astrophysical information. Those (d,p) reactions may be viewed as three-body reactions and described with Faddeev techniques. Faddeev-AGS equations in momentum space have a long tradition of utilizing separable interactions in order to arrive at sets of coupled integral equations in one variable. However, it needs to be demonstrated that observables calculated based on separable interactions agree exactly with those based on nonseparable forces. Momentum space AGS equations are solved with separable and nonseparable forces as coupled integral equations. Deuteron-α scattering is calculated via momentum space AGS equations using the CD-Bonn neutron-proton force and a Woods-Saxon type neutron(proton)- 4 He force, for which the Pauli-forbidden S-wave bound state is projected out. Elastic as well as breakup observables are calculated and compared to results in which the interactions in the two-body sub-systems are represented by separable interactions derived in the Ernst-Shakin-Thaler (EST) framework. We find that the calculations based on the separable representation of the interactions and the original interactions give results that are in excellent agreement. Specifically, integrated cross sections and angular distributions for elastic scattering agree within ≈ 1 %, which is well below typical experimental errors. In addition, the fivefold differential cross sections corresponding to breakup of the deuteron agree extremely well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Diffusive-to-ballistic transition in grain boundary motion studied by atomistic simulations

An adapted simulation method is used to systematically study grain boundary motion at velocities and driving forces across more than five orders of magnitude. This analysis reveals that grain boundary migration can occur in two modes, depending upon the temperature ( T ) and applied driving force ( P ). At low P and T , grain boundary motion is diffusional, exhibiting the kinetics of a thermally activated system controlled by grain boundary self-diffusion. At high P and T , grain boundary migration exhibits the characteristic kinetic scaling behavior of a ballistic process. A rather broad transition range in both P and T lies between the regimes of diffusive and ballistic grain boundary motion, and is charted here in detail. In conclusion, the recognition and delineation of these two distinct modes of grain boundary migration also leads to the suggestion that many prior atomistic simulations might have probed a different kinetic regime of grain boundary motion (ballistic) as compared to that revealed in most experimental studies (diffusional).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Structural and functional impact of troponin C-mediated Ca 2+ sensitization on myofilament lattice spacing and cross-bridge mechanics in mouse cardiac muscle

Acto-myosin cross-bridge kinetics are important for beat-to-beat regulation of cardiac contractility; however, physiological and pathophysiological mechanisms for regulation of contractile kinetics are incompletely understood. Here we explored whether thin filament-mediated Ca 2+ sensitization influences cross-bridge kinetics in permeabilized, osmotically compressed cardiac muscle preparations. We used a murine model of hypertrophic cardiomyopathy (HCM) harboring a cardiac troponin C (cTnC) Ca 2+ -sensitizing mutation, Ala8Val in the regulatory N-domain. We also treated wild-type murine muscle with bepridil, a cTnC-targeting Ca 2+ sensitizer. Our findings suggest that both methods of increasing myofilament Ca 2+ sensitivity increase cross-bridge cycling rate measured by the rate of tension redevelopment (k TR ); force per cross-bridge was also enhanced as measured by sinusoidal stiffness and I 1,1 -1/I 1,0 ratio from X-ray diffraction. Computational modeling suggests that Ca 2+ sensitization through this cTnC mutation or bepridil accelerates k TR primarily by promoting faster cross-bridge detachment. To elucidate if myofilament structural rearrangements are associated with changes in k(TR), we used small angle X-ray diffraction to simultaneously measure myofilament lattice spacing and isometric force during steady-state Ca 2+ activations. Within in vivo lattice dimensions, lattice spacing and steady-state isometric force increased significantly at submaximal activation. We conclude that the cTnC N-domain controls force by modulating both the number and rate of cycling cross-bridges, and that the both methods of Ca 2+ sensitization may act through stabilization of cTnC's D-helix. Furthermore, we propose that the transient expansion of the myofilament lattice during Ca 2+ activation may be an additional factor that could increase the rate of cross-bridge cycling in cardiac muscle. These findings may have implications for the pathophysiology of HCM.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Excitation of Soliton-Type Waves in Crystals of the A{sub 3}B Stoichiometry

Using the method of molecular dynamics and taking Ni{sub 3}Al and Pt{sub 3}Al as examples, crystals of the A{sub 3}B composition are considered for the possibility of excitation of soliton-type waves in them. The potentials obtained by the embedded-atom method were used to describe interatomic interactions. It is shown that the harmonic external stimulus can excite waves of the soliton type in a Pt{sub 3}Al crystal, but not in Ni{sub 3}Al. Such compression–expansion waves are generated because of excitation of discrete breathers with soft nonlinearity that cannot exist in a Ni{sub 3}Al crystal near the affected region. The detected waves are capable of propagating to thousands of nanometers along the Pt{sub 3}Al crystal without losses of integrity and speed. The shape of the obtained wave corresponds to the kink solution of the sine-Gordon equation. The aggregate amount of energy transferred by a wave is determined by the number of rows of atoms involved in fluctuations; this may involve dozens and hundreds of electron volts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Dynamics in the quantum/classical limit based on selective use of the quantum potential

A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed “quantum,” defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.

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

Electrochemical Synthesis of Monodisperse Platinum-Cobalt Nanocrystals

The synthesis of platinum-cobalt nanocrystals (NCs) using colloidal solvothermal techniques is well understood. However, for monodisperse NCs to form, high temperatures and environmentally detrimental solvents are needed. We report a room temperature, aqueous method of platinum-cobalt NC synthesis using electrochemical reduction as the driving force for nucleation and growth. It is found that colloidal NCs will form in both the presence and absence of surfactant. Additionally, we report a monodisperse electrochemical deposition of NCs utilising a transparent conducting oxide electrode. The methods developed here will allow for a synthetic method to produce nanocatalysts with minimal environmental impact and should be readily applicable to other NC systems, including single- and multi-component alloys.

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

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Rapid Freeform Sheet Metal Forming: Technology Development and System Verification

The objective of this project is to develop a transformational RApid Freeform sheet metal Forming Technology (RAFFT) in an industrial environment, which has the potential to increase manufacturing energy efficiency up to ten times, at a fraction of the cost of conventional technologies. The RAFFT technology is a flexible and energy-efficient process that eliminates the need for having geometry-specific forming dies. The innovation lies in the idea of using the energy resource at the local deformation area which provides greater formability, process control, and process flexibility relative to traditional methods. Double-Sided Incremental Forming (DSIF), the core technology in RAFFT, is a new concept for sheet metal forming. A blank sheet is clamped around its periphery and gradually deformed into a complex 3D freeform part by two strategically aligned stylus-type tools that follow a pre-described toolpath. The two tools, one on each side of the blank, can form a part with sharp features for both concave and convex shapes. Since deformation happens locally, the forming force at any instant is significantly decreased when compared to traditional methods. The key advantages of DSIF are its high process flexibility, high energy-efficiency, low capital investment, and the elimination of the need for massive amounts of die casting and machining. Additionally, the enhanced formability and process flexibility of DSIF can open up design spaces and result in greater weight savings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: Using density functional theory to accelerate wavefunction methods

The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modeling the quantum nature of the nuclei, the cost is still prohibitive when combined with advanced electronic structure methods. Here we present how multiple time step integrators can be combined with ring-polymer contraction techniques (effectively, multiple time stepping in imaginary time) to reduce virtually to zero the overhead of modelling nuclear quantum effects, while describing inter-atomic forces at high levels of electronic structure theory. This is demonstrated for a combination of MP2 and semi-local DFT applied to the Zundel cation. The approach can be seamlessly combined with other methods to reduce the computational cost of path integral calculations, such as high-order factorizations of the Boltzmann operator or generalized Langevin equation thermostats.

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

Spectral neighbor analysis method for automated generation of quantum-accurate interatomic potentials

We present a new interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projected onto a basis of hyperspherical harmonics in four dimensions. The bispectrum components are the same bond-orientational order parameters employed by the GAP potential [1]. The SNAP potential, unlike GAP, assumes a linear relationship between atom energy and bispectrum components. The linear SNAP coefficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum components. The calculation of the bispectrum components and the SNAP potential are implemented in the LAMMPS parallel molecular dynamics code. We demonstrate that a previously unnoticed symmetry property can be exploited to reduce the computational cost of the force calculations by more than one order of magnitude. We present results for a SNAP potential for tantalum, showing that it accurately reproduces a range of commonly calculated properties of both the crystalline solid and the liquid phases. In addition, unlike simpler existing potentials, SNAP correctly predicts the energy barrier for screw dislocation migration in BCC tantalum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗
DOE OSTI2019

Future climate emulations using quantile regressions on large ensembles

The study of climate change and its impacts depends on generating projections of future temperature and other climate variables. For detailed studies, these projections usually require some combination of numerical simulation and observations, given that simulations of even the current climate do not perfectly reproduce local conditions. We present a methodology for generating future climate projections that takes advantage of the emergence of climate model ensembles, whose large amounts of data allow for detailed modeling of the probability distribution of temperature or other climate variables. The procedure gives us estimated changes in model distributions that are then applied to observations to yield projections that preserve the spatiotemporal dependence in the observations. We use quantile regression to estimate a discrete set of quantiles of daily temperature as a function of seasonality and long-term change, with smooth spline functions of season, long-term trends, and their interactions used as basis functions for the quantile regression. A particular innovation is that more extreme quantiles are modeled as exceedances above less extreme quantiles in a nested fashion, so that the complexity of the model for exceedances decreases the further out into the tail of the distribution one goes. We apply this method to two large ensembles of model runs using the same forcing scenario, both based on versions of the Community Earth System Model (CESM), run at different resolutions. The approach generates observation-based future simulations with no processing or modeling of the observed climate needed other than a simple linear rescaling. The resulting quantile maps illuminate substantial differences between the climate model ensembles, including differences in warming in the Pacific Northwest that are particularly large in the lower quantiles during winter. We show how the availability of two ensembles allows the efficacy of the method to be tested with a “perfect model” approach, in which we estimate transformations using the lower-resolution ensemble and then apply the estimated transformations to single runs from the high-resolution ensemble. Finally, we describe and implement a simple method for adjusting a transformation estimated from a large ensemble of one climate model using only a single run of a second, but hopefully more realistic, climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗