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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Derivation of a realistic forcing term to reproduce the turbulent characteristics of round jets on the centerline

Turbulence forcing techniques are often required in the numerical simulation of statistically stationary turbulent flows. However, the existing forcing techniques are not based on physics, but rather arbitrary numerical methods that sustain the turbulent kinetic energy. In this work, a forcing technique is devised to reproduce the centerline turbulent characteristics of round jets in a triply periodic box. It is derived from the Navier-Stokes equations by applying aReynolds decompositionwith themean velocity of the axisymmetric jet. The result is an anisotropic linear forcing term, which is intended to be used in a three-dimensional box to create turbulence. Four direct numerical simulations with different Re λ have been performed with these forcing terms. We report the budget of the terms in the kinetic energy equation is very close to the experimental measurement on the centerline. The anisotropy, kinetic energy k, and dissipation rate ε of the simulations are also comparable to experimental values. Finally, the kinetic energy spectrum in the axial direction, $\phi(κ_1)$, is presented. With appropriate normalizations, the spectrum agrees well with the round jet spectrum on its centerline.

42 ENGINEERING↗
DOE OSTI2022

Calculation of Nuclear Reactor Cooling Tower Performance With Limited Data Streams

Monitoring of cooling tower performance in a nuclear reactor facility is necessary to ensure safe operation; however, instrumentation for measuring performance characteristics can be difficult to install and may malfunction or break down over long duration experiments. This paper describes employing a thermodynamic approach to quantify cooling tower performance, the Merkel model, which requires only five parameters, namely, inlet water temperature, outlet water temperature, liquid mass flowrate, gas mass flowrate, and wet bulb temperature. Using this model, a general method to determine cooling tower operation for a nuclear reactor was developed in situations when neither the outlet water temperature nor gas mass flowrate are available, the former being a critical piece of information to bound the Merkel integral. Furthermore, when multiple cooling tower cells are used in parallel (as would be in the case of large-scale cooling operations), only the average outlet temperature of the cooling system is used as feedback for fan speed control, increasing the difficulty of obtaining the outlet water temperature for each cell. To address these shortcomings, this paper describes a method to obtain individual cell outlet water temperatures for mechanical forced-air cooling towers via parametric analysis and optimization. In this method, the outlet water temperature for an individual cooling tower cell is acquired as a function of the liquid-to-gas ratio (L/G). Leveraging the tight tolerance on the average outlet water temperature, an error function is generated to describe the deviation of the parameterized L/G to the highly controlled average outlet temperature. The method was able to determine the gas flowrate at rated conditions to be within 3.9% from that obtained from the manufacturer’s specification, while the average error for the four individual cooling cell outlet water temperatures were 1.6 °C, -0.5 °C, -1.0 °C, and 0.3 °C.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Curved-line search algorithm for ab initio atomic structure relaxation

Ab initio atomic relaxations often take large numbers of steps and long times to converge, especially when the initial atomic configurations are far from the local minimum or there are curved and narrow valleys in the multidimensional potentials. An atomic relaxation method based on on-the-flight force learning and a corresponding curved-line search algorithm is presented to accelerate this process. In this work, results demonstrate the superior performance of this method for metal and magnetic clusters when compared with the conventional conjugate-gradient method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Advancing a Model-Validated Statistical Method for Decomposing the Key Oceanic Drivers of Regional Climate: Focus on Northern and Tropical African Climate Variability in the Community Earth System Model (CESM)

This study advances the practicality and stability of the traditional multivariate statistical method, generalized equilibrium feedback assessment (GEFA), for decomposing the key oceanic drivers of regional atmospheric variability, especially when available data records are short. An advanced stepwise GEFA methodology is introduced, in which unimportant forcings within the forcing matrix are eliminated through stepwise selection. Method validation of stepwise GEFA is performed using the CESM, with a focused application to northern and tropical Africa (NTA). First, a statistical assessment of the atmospheric response to each primary oceanic forcing is carried out by applying stepwise GEFA to a fully coupled control run. Then, a dynamical assessment of the atmospheric response to individual oceanic forcings is performed through ensemble experiments by imposing sea surface temperature anomalies over focal ocean basins. Finally, to quantify the reliability of stepwise GEFA, the statistical assessment is evaluated against the dynamical assessment in terms of four metrics: the percentage of grid cells with consistent response sign, the spatial correlation of atmospheric response patterns, the area-averaged seasonal cycle of response magnitude, and consistency in associated mechanisms between assessments. In CESM, tropical modes, namely El Niño–Southern Oscillation and the tropical Indian Ocean Basin, tropical Indian Ocean dipole, and tropical Atlantic Niño modes, are the dominant oceanic controls of NTA climate. In complementary studies, stepwise GEFA is validated in terms of isolating terrestrial forcings on the atmosphere, and observed oceanic and terrestrial drivers of NTA climate are extracted to establish an observational benchmark for subsequent coupled model evaluation and development of process-based weights for regional climate projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Energetics of the nucleation and glide of disconnection modes in symmetric tilt grain boundaries

Grain boundaries (GBs) evolve by the nucleation and glide of disconnections, which are dislocations with a step character. In this work, motivated by recent success in predicting GB properties such as the shear coupling factor and mobility from the intrinsic properties of disconnections, we develop a systematic method to calculate the energy barriers for the nucleation and glide of individual disconnection modes under arbitrary driving forces and a quasi-2D setting. This method combines tools from bicrystallography to enumerate disconnection modes and the Nudged elastic band (NEB) method to calculate their energetics, yielding minimum energy paths and atomistic mechanisms for the nucleation and glide of each disconnection mode. We apply the method to accurately predict shear coupling factors of $[001]$ symmetric tilt grain boundaries in Cu. Particular attention is paid to the boundaries where the dislocation-based disconnection nucleation model produces incorrect nucleation barriers. We demonstrate that the method can accurately compute energy barriers and predict shear-coupling factors in the low-temperature regime. For certain disconnection modes in which the assumptions underlying our method do not hold, we report upper bounds on the energy barriers for disconnection nucleation and glide. In addition, the NEB trajectories reveal interesting phenomena such as the dissociation of a higher energy mode into lower energy modes, and in some cases, shear coupling being mediated by partial disconnections, wherein the GB structure temporarily changes to a metastable state before reverting back to its original structure. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Continuous-Flow Centrifugal Solid/Liquid Separation for the Recovery of Rare-Earth Elements Containing Particles from Phosphoric Acid Sludge

Phosphoric acid sludge contains acid (~54% P 2 O 5 ) and solid precipitates including rare earth elements (REEs) at concentrations of ~2200 ppm. Low-cost recovery of valuable P 2 O 5 and simultaneous solid separation could be an economically feasible approach to recovering REEs while increasing the production of phosphoric acid. The sludge, however, is a complicated stream that cannot be separated by traditional technologies because of high viscosity and a large solid content (30–40%). Guided by a force balance model, an efficient solid/liquid separation method is demonstrated, using a continuous-flow centrifugal contactor. Additionally, the shear regime is bypassed by introducing the sludge directly into the rotor where a centrifugal force is exerted on the fluid, inducing phase separation. Solid particles are trapped in the contactor. High liquid recovery is demonstrated and the effects of process parameters on solid capture are investigated. Three contactors in series yield 94% recovery of solids containing 1500–2895 ppm of REEs.

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

Force Matching Approaches to Extend Density Functional Theory to Large Time and Length Scales

We present methods for the creation of semi-empirical quantum approaches and reactive force fields through force matching to quantum simulation data for materials under reactive conditions. Our methodologies overcome the extreme computational cost of standard Kohn-Sham Density Functional Theory (DFT) by mapping DFT computed simulation data onto functional forms with linear dependence on their parameters. This allows for quick parameterization of our models by avoiding the non-linear fitting bottlenecks associated with most molecular dynamics model development. We illustrate our approach with two different systems: (i) determination of density functional tight binding models for aqueous glycine dimerization, and (ii) determination of the Chebyshev Interactional Model for Efficient Simulation (ChIMES) reactive force field for metallic liquid carbon. In each case, we observe that our approach is easy to parametrize and yields a model that is orders of magnitude faster than DFT while largely retaining its accuracy. Overall, our methods have potential use for studying complex long time and length scale chemical reactivity at extreme conditions, where there is a significant need for computationally efficient atomistic simulations methods to aid in the interpretation and design of experiments.

Lindsey, R.↗
DOE OSTI2018

Few-nucleon and many-nucleon systems with semilocal coordinate-space regularized chiral nucleon-nucleon forces

We employ a variety of ab initio methods, including Faddeev-Yakubovsky equations, no-core configuration interaction approach, coupled-cluster theory, and in-medium similarity renormalization group, to perform a comprehensive analysis of the nucleon-deuteron elastic and breakup reactions and selected properties of light and medium-mass nuclei up to 48 C a using the recently constructed semilocal coordinate-space regularized chiral nucleon-nucleon potentials. We compare the results with those based on selected phenomenological and chiral EFT two-nucleon potentials, discuss the convergence pattern of the chiral expansion, and estimate the achievable theoretical accuracy at various chiral orders using an approach to quantify truncation errors of the chiral expansion without relying on cutoff variation. We also address the robustness of this method and explore alternative ways to estimate the theoretical uncertainty from the truncation of the chiral expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Mode-synthesizing atomic force microscopy and mode-synthesizing sensing

A method of analyzing a sample that includes applying a first set of energies at a first set of frequencies to a sample and applying, simultaneously with the applying the first set of energies, a second set of energies at a second set of frequencies, wherein the first set of energies and the second set of energies form a multi-mode coupling. The method further includes detecting an effect of the multi-mode coupling.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

The ABACUS cosmological N -body code

We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct 1/r 2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10 -5 . We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Molecular simulation of fluids with non-identical intermolecular potentials: Thermodynamic properties of 10-5 + 12-6 Mie potential binary mixtures

General methods for combining interactions between particles characterised by non-identical intermolecular potentials are investigated. The combination methods are tested by performing molecular dynamics simulations to determine the pressure, energy, isochoric and isobaric heat capacities, thermal expansion coefficient, isothermal compressibility, Joule-Thomson coefficient, and speed of sound of 10-5 + 12-6 Mie potential binary mixtures. In addition to the two non-identical Mie potentials, mixtures are also studied with non-identical intermolecular parameters. The combination methods are compared with results obtained by simply averaging the Mie exponents. When either the energy or size parameters are non-identical, very significant differences emerge in the thermodynamic properties predicted by the alternative combination methods. The isobaric heat capacity is the thermodynamic property that is most affected by the relative magnitude of the intermolecular potential parameters and the method for combining non-identical potentials. Either the arithmetic or geometric combination of potentials provides a simple and effective way of performing simulations involving mixtures of components characterised by non-identical intermolecular potentials, which is independent of their functional form.

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

Method of mechanical holding of cantilever chip for tip-scan high-speed atomic force microscope

In tip-scan atomic force microscopy (AFM) that scans a cantilever chip in the three dimensions, the chip body is held on the Z-scanner with a holder. However, this holding is not easy for high-speed (HS) AFM because the holder that should have a small mass has to be able to clamp the cantilever chip firmly without deteriorating the Z-scanner’s fast performance, and because repeated exchange of cantilever chips should not damage the Z-scanner. This is one of the reasons that tip-scan HS-AFM has not been established, despite its advantages over sample stage-scan HS-AFM. Here, we present a novel method of cantilever chip holding which meets all conditions required for tip-scan HS-AFM. The superior performance of this novel chip holding mechanism is demonstrated by imaging of the α{sub 3}β{sub 3} subcomplex of F{sub 1}-ATPase in dynamic action at ∼7 frames/s.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Anomalous random correlations of force constants on the lattice dynamical properties of disordered Au-Fe alloys

Here, gold iron (Au-Fe) alloys are of immense interest due to their biocompatibility, anomalous Hall conductivity, and applications in various medical treatments. However, irrespective of the method of preparation, they often exhibit a high level of disorder with properties sensitive to the thermal or magnetic annealing temperatures. We calculate the lattice dynamical properties of Au 1–x Fe x alloys using density functional theory methods where, being multisite properties, reliable interatomic force constant (IFC) calculations in disordered alloys remain a challenge. We follow a twofold approach: an accurate IFC calculation in an environment with nominally zero chemical pair correlations to mimic the homogeneously disordered alloy and a configurational averaging for the desired phonon properties (e.g., dispersion, density of states, and entropy). We find an anomalous change in the IFC's and phonon dispersion (split bands) near x=0.19, which is attributed to the local stiffening of the Au-Au bonds when Au is in the vicinity of Fe. Other results based on mechanical and thermophysical properties reflect a similar anomaly: Phonon entropy, e.g., becomes negative below x=0.19, suggesting a tendency for chemical unmixing, reflecting the onset of a miscibility gap in the phase diagram. Our results match fairly well with reported data wherever available.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High-coercivity magnetism in nanostructures with strong easy-plane anisotropy

Here, we report the fabrication of a rare-earth-free permanent-magnet material Co 3 Si in the form of nanoparticles and investigate its magnetic properties by experiments and density-functional theory (DFT). The DFT calculations show that bulk Co 3 Si has an easy-plane anisotropy with a high K 1 ≈ -64 Merg/cm 3 (-6.4 MJ/m 3 ) and magnetic polarization of 9.2 kG (0.92 T). In spite of having a negative anisotropy that generally leads to negligibly low coercivities in bulk crystals, Co 3 Si nanoparticles exhibit high coercivities (17.4 kOe at 10K and 4.3 kOe at 300 K). This result is a consequence of the unique nanostructure made possible by an effective easy-axis alignment in the cluster-deposition method and explained using micromagnetic analysis as a nanoscale phenomenon involving quantum-mechanical exchange interactions.

36 MATERIALS SCIENCE↗