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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Gravitational Effects on Distortion in Sintering

During sintering a powder compact gains strength through low-temperature interparticle bonding, usually induced by solid-state surface diffusion, followed by further strength contributions from high-temperature densification. In cases where a liquid phase forms, sintering densification is accelerated and shape retention is sustained while open pores remain and contribute capillary forces. Unfortunately, sintering densification requires the compact become thermally softened to a point where creep strain rates reach levels near 10(exp -2)/s when the liquid forms. On the other hand, thermal softening of the powder compact substantially reduces the strength at high temperatures. Therefore, the in situ strength evolution during sintering is a primary focus to separate compact densification (as required for high performance) with minimized distortion (as required for net-shaping). With respect to gravitation effects on distortion during sintering there are two points of substantial weakness - prior to significant interparticle bonding and during final pore closure. This research is focused on understanding the competition among interparticle neck growth, densification, thermal softening, grain boundary wetting, capillary effects associated with liquid wetting and residual porosity, and gravity. Most surprising is the apparent role of gravity, where the deviatoric stress acting on the powder structure induces skeletal formation that reduces distortion. In contrast with theory, microgravity samples exhibit more distortion yet fail to fully densify. Results are presented on the experimental concepts supporting an emerging model of sintering strength evolution that enables understanding of both distortion and densification. The experiments have relied on tungsten heavy alloys, various combinations of dihedral angle, pore size, initial porosity, liquid:solid ratio, and heating rates. On Earth, the dominant factor with respect to distortion is the starting body heterogeneity. Current modeling efforts are seeking some means to uniformly predict the distortion based on a starting pore structure heterogeneity parameter. Densification is largely unaffected by the initial pore structure, but distortion is highly variable, suggesting that nonuniform pore closure might be a significant parameter during sintering. With respect to flight experiments, plans include removal of the solid body forces acting on the solid grains, allowing stabilization of the pore structure and examination of the buoyancy effects with regard to distortion. In microgravity there is the surprising result that compacts do not densify, yet distort a factor that is contrary to all current sintering models. Densification without distortion during liquid phase sintering was achieved by manipulating microstructure and its evolution during sintering. Microstructure parameters such as the solid volume fraction, dihedral angle, initial porosity, and pore size were varied to measure densification and distortion behavior during LPS using W-Ni-Cu alloys. Green compacts were formed using ethylene-bis-stearamide as a pore-forming agent with the amount of polymer controlling the initial porosity. Different initial pore sizes were generated by varying the polymer particle size. Dihedral angle was varied by changing the Ni:Cu ratio in the alloys. Finally, the solid volume fraction was adjusted via the tungsten content. Distortion was quantified using profiles determined with a coordinate measuring machine to calculate a distortion parameter. Sintering results showed that solid volume fraction and dihedral angle are the dominant factors on densification and distortion during liquid phase sintering. Distortion decreases with increasing solid volume fraction and dihedral angle, while initial porosity and pore size have no observable effect on distortion at nearly full densification. Various strategies emerge to improve distortion control in liquid phase sintering.

German, Randall M.↗
DOE OSTI2026

Toward scalable quantum computations of atomic nuclei

We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3 He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2 H and 3 He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Furthermore, based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
DOE OSTI2023

Generalized filtered lifting line theory for arbitrary chord lengths and application to wind turbine blades

Abstract The filtered lifting line theory is an analytical approach used to solve the equations of flow subjected to body forces with a Gaussian distribution, such as used in the actuator line model. In the original formulation, the changes in chord length along the blade were assumed to be small. This assumption can lead to errors in the induced velocities predicted by the theory compared to full solutions of the equations. In this work, we revisit the original derivation and provide a more general formulation that can account for significant changes in chord along the blade. The revised formulation can be applied to wings with significant changes in chord along the span, such as wind turbine blades.

17 WIND ENERGY↗
DOE OSTI2023

Computational studies of impurity migration during induction stirring of molten uranium

Understanding and controlling impurity behavior in actinide metal casting processes are foundational for efficient part production yet remain major challenges for researchers and industry. To help provide insight regarding impurity distribution during actinide metal casting, we have developed computational fluid dynamics (CFD) models for a laboratory-scale system using commercial and open-source code. Multiple simulation frameworks allow for improved confidence in the resulting outputs, while taking advantage of the maturity and convenience offered by commercial platforms and simultaneously maintaining the transparency and flexibility often provided by open-source software. Here, we describe multiple experiment-informed models designed to simulate a specific laboratory system in which uranium melt, containing a known starting concentration of carbon impurity, is electromagnetically stirred in an induction furnace. Here, the goal of the simulations is to predict the motion of uranium carbide microparticles in the velocity field of the melt. Prior to simulating the uranium-carbon system, numerical models were validated using a previously published nonradioactive experimental system. Effects of the size and shape of impurity particles in the models were investigated and agree with experimental findings. Simulation of smaller particles (< 50 µm) shows more homogenous distribution throughout the stirred melt. With increased particle size (100 µm), the body forces, which include the buoyancy force, dominate over the drag force, causing larger particles to move toward the crucible walls and upward in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A generalized analytical energy balance model for evaluating agglomeration from a binary collision of wet particles

Agglomeration of wet particles, i.e., particles coated with a thin liquid layer, is a common phenomenon in many processes like fluidized bed combustion of low rank fuels. The availability of an agglomeration model that can evaluate the outcome of a binary collision between wet particles differing in solid particle properties, liquid layer thicknesses, and initial collision (impact) speeds is essential for obtaining a comprehensive understanding on the existing processes experiencing wet particle agglomeration or for a successful development of new processes with high chances of wet particle agglomeration. This study presents a generalized agglomeration model on the basis of energy conservation before and after collision when colliding wet particles may differ in solid particle properties, liquid layer thicknesses, and impact speeds. The model was established based on the approximate values of energy losses that may happen during the collision. It incorporates body forces, solid-solid contacting, liquid capillary, and viscous contributions, as well as the liquid bridge volume effect. Predictions of the new model for collision outcomes of identical wet particles were like those from an analytical energy balance model developed recently by the group for identical wet particles. We also validated the new model by experimental data from literature. The results of a collision direction analysis indicated that the direction often has a minimal effect on the collision outcome in many practical scenarios. The results of Monte Carlo uncertainty analyses with the new model revealed that proper estimations of impact speed, under capillary limiting conditions, and thickness of coating layers and asperity heights, under viscous limiting conditions, are critical for the realistic prediction of collision outcomes at impact speeds close to critical impact speed, i.e., the minimum particle speed required for the particles to rebound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Thermoosmotic Transport in Nanochannels Grafted with pH-responsive Polyelectrolyte Brushes Modelled Using Augmented Strong Stretching Theory

In this paper, we develop a theory to establish that the thermoosmotic (TOS) effects, induced by the application of an axial temperature gradient, leads to a massive enhancement in liquid transport in nanochannels grafted with charged polyelectrolyte (PE) brushes. Here, we quantify the TOS transport by quantifying the induced electric field and the induced TOS flow field. The different components of the electric field, namely the ionic component, the thermal component, and the osmotic component, as well as the contributions of different ions to these components, are quantified. Furthermore, we express the TOS velocity as a combination of chemiosmotic (COS), thermal, and electroosmotic (EOS) components. The COS and the thermal components augment each other and the overall strength and direction of the TOS flow are dictated by the direction and the relative strength of the EOS component. Most importantly, we compare the cases of brush-grafted nanochannels with those of the brush-free nanochannels of identical surface charge densities: the TOS transport is massively augmented in the brush-grafted nanochannels attributed to the combination of the localization of the EDL (and hence any body force that depends on the EDL charge density) away from the nanochannel wall (i.e., the location of the maximum drag force) and the presence of a possible molecular slip (experienced by the liquid) along the brush surface.

42 ENGINEERING↗
DOE OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗
DOE OSTI2025

Spherical and sessile droplet dynamics by fluctuating hydrodynamics

We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.

Bell, John B↗
DOE OSTI2023

Toward precision Fermi-liquid theory in two dimensions

The ultracold and weakly coupled Fermi gas in two spatial dimensions is studied in an effective field theory framework. It has long been observed that universal corrections to the energy density to two orders in the interaction strength do not agree with Monte Carlo simulations in the weak-coupling regime. Here, universal corrections to three orders in the interaction strength are obtained, and are shown to provide agreement between theory and simulation. Special consideration is given to the scale ambiguity associated with the nontrivial renormalization of the singular contact interactions. Further, the isotropic superfluid gap is obtained to next-to-leading order, and nonuniversal contributions to the energy density due to effective range effects, p -wave interactions, and three-body forces are computed. Results are compared with precise Monte Carlo simulations of the energy density and the contact in the weakly coupled attractive and repulsive Fermi-liquid regimes. In addition, the known all-orders sum of ladder and ring diagrams is compared with Monte Carlo simulations at weak coupling and beyond.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗
DOE OSTI2025

Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

active matter↗
DOE OSTI2025

Deviations from the Isobaric Multiplet Mass Equation due to Threshold States

Recent studies have completed the A=16 isospin quintets for states with J π =0 + and 2 + . The dependence of their masses as a function of isospin projection shows evidence for deviations from quadratic behavior indicating isospin violation beyond the expectation from two-body forces. The deviation is most pronounced for the 2 + states. Predictions from the shell model embedded in the continuum (SMEC) allow us to explain that this isospin violation is associated with a modification of the nuclear structure due to the open-quantum-system nature of the proton-rich members of the quintet. In particular, the 0 + and 2 + states in 16 Ne and the 2 + state in 16 F are threshold resonances located just above a proton-decay threshold where s-wave coupling to the continuum is expected. The measured deviations of these threshold states from the quadratic behavior of the remaining members of the multiplets makes it possible to obtain information on the magnitude and the energy dependence of the continuum-coupling energy correction. Finally, continuum coupling is also indicated for the ground state of 8 C, but this time through p-wave coupling

coulomb energies & analogue states↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2024

A Solution Method for the Filtered Lifting Line Theory

The filtered lifting line theory presents a continuous form of the inviscid momentum equations of flow over a lifting device, such as a wing or rotor blade, using body forces without mathematical singularities. This theory is also consistent with an actuator line representation of a lifting device. In this work, we present a reformulation of the equations in terms of the local flow angle along the line, which allows solving the stand-alone equations using multivariate root-finding algorithms. This approach can be used to obtain a fast, computationally inexpensive solution of the loading distribution along a wing without the need to perform computational fluid dynamic simulations. We study the requirements in terms of resolution in the spanwise direction and establish the criteria for spacing and minimum amount of points required along the blade to obtain converged solutions. The solutions are compared to results from large-eddy simulations, and we observed excellent agreement with less than a percent difference in quantities along the blade between the methods.

17 WIND ENERGY↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗