Search NASA⌕ Search

SEARCH · Search NASA

Results for “Body 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

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

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials

In this study, we present the equation of state of infinite neutron matter as obtained from highly realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne v 18 and the Argonne v' 6 nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ( ρ 0 =0.16 fm ⁻3 ) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and nd doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium β decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Mass-Energy Compensation Effect of 3$\alpha$ Hamiltonian

The 3α phenomenological model describes the structure of the carbon-12 nucleus as a cluster of three alpha particles. This model includes a pairwise α–α interaction and a three-body force. To fit the three-body potential, the 12 C data are used, while ensuring that the pair potential reproduces the α–α scattering data. Alternatively, the mass-energy compensation (MEC) effect can be used to simulate the effect of the three-body potential by adjusting the mass of the α particle within the effective-mass approach. We demonstrate the MEC effect for the 3α ground state by numerically solving the differential Faddeev equation, in which the α–α interaction is described by the Ali-Bodmer potential. The effective masses of α particles are evaluated for the ground and excited 0 + and bound 2 + states. Here, we demonstrate a coupling between the ground and first excited 0 + states, indicated by an anti-crossing of these energy levels in the energy–mass coordinates. A correspondence between the effective mass and a three-body potential is demonstrated. We discuss the results of the 0$^{+}_{2}$ calculations for various models of the α–α interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

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

Elastic proton scattering off nonzero spin nuclei

In recent years, we constructed a microscopic optical potential (OP) for elastic nucleon-nucleus ( NA ) scattering using modern approaches based on chiral theories for the nucleon-nucleon ( NN ) interaction. The OP was derived at first order of the spectator expansion in Watson multiple scattering theory and its final expression was a folding integral between the NN t matrix and the nuclear density of the target. Two- and three-body forces are consistently included both in the target and in the projectile description. The purpose of this work is to apply our microscopic OP to nuclei characterized by a ground state of spin-parity quantum numbers J π ≠ 0 + . We extended our formalism to include the spin of the target nucleus. The full amplitudes of the NN reaction matrix are retained in the calculations starting from two- and three-body chiral forces. The microscopic OP can be applied in the energy range 100 ≤ E ≤ 350 MeV. Finally, we show a remarkable agreement with experimental data for the available observables and, simultaneously, provide reliable estimates for the theoretical uncertainties. This work paves the way toward a full microscopic approach to inelastic NA scattering, showing that the derivation of optical potentials between states with J π ≠ 0 + is completely under control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Tensor force role in β decays analyzed within the Gogny-interaction shell model

The half-life of the famous C 14 β decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. In this study, we revisit and analyze the role of the tensor force for the β decay of C 14 as well as of neighboring isotopes. We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a p – s d model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum L in order to analyze the effects of the tensor force and cross-shell mixing. The inclusion of the tensor force significantly improves the shell-model calculations of the β -decay properties of carbon isotopes. In particular, the anomalously slow β decay of C 14 can be explained by the isospin T = 0 part of the tensor force, which changes the components of N 14 with the orbital angular momentum L = 0 , 1 , and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby β decays are improved. Decomposition of wave function into L components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of β -decay matrix elements. Cross-shell mixing also has a visible impact on the β -decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stemflow Hydrodynamics

Stemflow hydrodynamics is the study of water movement along the exterior surface area of plants. Its primary goal is to describe water velocity and water depth along the stem surface area. Its significance in enriching the rhizosphere with water and nutrients is not in dispute. Yet, the hydrodynamics of stemflow have been entirely overlooked. This review seeks to fill this knowledge gap by drawing from thin film theories to seek outcomes at the tree scale. The depth‐averaged conservation equations of water and solute mass are derived at a point. These equations are then supplemented with the conservation of momentum that is required to describe water velocities or relations between water velocities and water depth. Relevant forces pertinent to momentum conservation are covered and include body forces (gravitational effects), surface forces (wall friction), line forces (surface tension), and inertial effects. The inclusion of surface tension opens new vistas into the richness and complexity of stemflow hydrodynamics. Flow instabilities such as fingering, pinching of water columns into droplets, accumulation of water within fissures due to surface tension and their sudden release are prime examples that link observed spatial patterns of stemflow fronts and morphological characteristics of the bark. Aggregating these effects at the tree‐ and storm‐ scales are featured using published experiments. The review discusses outstanding challenges pertaining to stemflow hydrodynamics, the use of dynamic similarity and 3D printing to enable the interplay between field studies and controlled laboratory experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Light nuclei with semilocal momentum-space regularized chiral interactions up to third order

In this work, we present a systematic investigation of few-nucleon systems and light nuclei using the current Low Energy Nuclear Physics International Collaboration interactions comprising semilocal momentum-space regularized two- and three-nucleon forces up to third chiral order (N 2 LO). Following our earlier study utilizing the coordinate-space regularized interactions, the two low-energy constants entering the three-body force are determined from the triton binding energy and the differential cross-section minimum in elastic nucleon-deuteron scattering. Predictions are made for selected observables in elastic nucleon-deuteron scattering and in the deuteron breakup reactions, for properties of the A = 3 and A = 4 nuclei, and for spectra of p-shell nuclei up to A = 16. A comprehensive error analysis is performed including an estimation of correlated truncation uncertainties for nuclear spectra. The obtained predictions are generally found to agree with experimental data within errors. Similarly to the coordinate-space regularized chiral interactions at the same order, a systematic overbinding of heavier nuclei is observed, which sets in for A~10 and increases with A.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence

Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Furthermore, even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98%. These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the falloff rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

van der Waals interaction as the starting point for an effective field theory

We consider the system of three 4 He atoms to assess whether a pure van der Waals potential can be used as a starting point for an effective field theory to describe three-body processes in ultracold atomic systems. Using a long-range van der Waals interaction in combination with short-distance two-body counterterms, we analyze the dependence of two- and three-body observables on the short-distance regulator that is required due to the singular nature of the van der Waals interaction. We benchmark our approach with results obtained with the realistic 4 He- 4 He LM2M2 potential and find good agreement. We furthermore show that in this effective field theory approach no three-body force is required at leading order and that universal van der Waals physics leads to a universal correlation between three-body observables in the absence of an Efimov three-body parameter.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Local position-space two-nucleon potentials from leading to fourth order of chiral effective field theory

We present local, position-space chiral NN potentials through four orders of chiral effective field theory ranging from leading order (LO) to next-to-next-to-next-to-leading order (N 3 LO, fourth order) of the Δ-less version of the theory. The long-range parts of these potentials are fixed by the very accurate πN LECs as determined in the Roy-Steiner equations analysis. At the highest order (N 3 LO), the NN data below 190 MeV laboratory energy are reproduced with the respectable χ 2 /datum of 1.45. A comparison of the N 3 LO potential with the phenomenological Argonne v 18 (AV18) potential reveals substantial agreement between the two potentials in the intermediate range ruled by chiral symmetry, thus, providing a chiral underpinning for the phenomenological AV18 potential. Our chiral NN potentials may serve as a solid basis for systematic ab initio calculations of nuclear structure and reactions that allow for a comprehensive error analysis. In particular, the order by order development of the potentials will make possible a reliable determination of the truncation error at each order. Our new family of local position-space potentials differs from existing potentials of this kind by a weaker tensor force as reflected in relatively low D-state probabilities of the deuteron (P D ≲ 4.0 % for our N 3 LO potentials) and predictions for the triton binding energy above 8.00 MeV (from two-body forces alone). As a consequence, our potentials may lead to different predictions when applied to light and intermediate-mass nuclei in ab initio calculations and, potentially, help solve some of the outstanding problems in microscopic nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shoulder fillet effects in strength distributions of microelectromechanical system components

The failure forces and fracture strengths of polysilicon microelectromechanical system (MEMS) components in the form of stepped tensile bars with shoulder fillets were measured using a sequential failure chain methodology. Approximately 150 specimens for each of four fillet geometries with different stress concentration factors were tested. The resulting failure force and strength distributions of the four geometries were related by a common sidewall flaw population existing within different effective stressed lengths. The failure forces, strengths, and flaw population were well described by a weakest-link based analytical framework. Finite element analysis was used to verify body-force based expressions for the stress concentration factors and to provide insight into the variation of specimen effective length with fillet geometry. Monte Carlo simulations of flaw size and location, based on the strength measurements, were also used to provide insight into fillet shape and size effects. The successful description of the shoulder fillet specimen strengths provides further empirical support for application of the strength and flaw framework in MEMS fabrication and design optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quasielastic $\overrightarrow{^{3}He}$ $(\overrightarrow{e},e′)$ asymmetry in the threshold region with a longitudinally polarized target

We report on the measurements of the double-spin asymmetry from electron- 3 He scattering in the threshold region of two- and three-body breakup of 3 He for 4-momentum transfer values of 0.1 and 0.2 (GeV/c) 2 . The measurement was done at Jefferson Lab using a polarized 3 He target with its spin aligned parallel to the beam direction. This asymmetry is sensitive to a combination of spin-dependent response functions not previously measured and may highlight differences in two body forces, meson-exchange currents and final-state interactions. The results of this measurement serve as a test of our understanding of few-body systems. When compared with calculations from plane wave impulse approximation and Faddeev calculations, we found that the latter, which use modern nuclear potentials and prescriptions for meson-exchange currents, demonstrate an overall good agreement with data.

Few-body nuclei↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗